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FIRE  AND  SMOKE  SIMULATOR  (FSSIM)  VERSION  1  -  USER’S  GUIDE 


1.0  INTRODUCTION 

1.1  Motivation  for  Developing  FSSIM 

In  both  peacetime  and  during  war,  fire  represents  a  significant  threat  to  any  ship.  A  fire,  whether 
started  by  a  mechanical  failure,  intentional  or  unintentional  human  actions,  or  damage  from  a 
weapon  hit,  threatens  the  ship  in  a  number  of  ways.  The  crew’s  health  and  ability  to  operate  the 
ship  is  affected  by  direct  exposure  to  the  fire  or  by  the  spread  of  smoke  and  toxic  gasses  through 
the  ship  by  either  natural  or  mechanical  ventilation.  Electrical  systems  can  be  degraded  by 
thermal  exposure,  exposure  to  acid  gasses  in  the  combustion  products,  or  by  electrical  failure 
resulting  from  soot  deposition,  which  might  include  hampered  cooling  or  dielectric  breakdown 
from  the  electrical  conductivity  of  the  soot.  Mechanical  systems  can  suffer  thermal  damage. 
Lastly,  fire  growth  and  spread  could  potentially  ignite  explosive  materials,  rocket  motors, 
aviation  fuel,  or  other  highly  flammable  substances,  which  could  possibly  result  in  temperatures 
or  overpressures  high  enough  to  affect  the  ship  structurally. 

There  currently  exist  a  number  of  computational  tools  for  examining  the  effects  of  a  fire  that  can 
be  applied  to  a  ship  and  its  crew.  One  could  use  hand  calculations  for  examining  simple 
scenarios  in  single  compartments.  Simple  rules  can  be  used  to  extend  this  approach  to  multiple 
compartments.  Zone  models  are  suitable  for  examining  somewhat  more  complex, 
time-dependent  scenarios  involving  multiple  compartments  and  levels,  but  stability  can  be  a 
problem  for  multi-level  scenarios,  scenarios  with  Heating,  Ventilation,  and  Air  Conditioning 
(HVAC)  systems,  and  for  post-flashover  conditions.  Computational  fluid  dynamics  (CFD) 
models  can  yield  detailed  information  about  temperatures,  heat  fluxes,  and  species 
concentrations;  however,  the  time  penalty  of  this  approach  currently  makes  using  CFD 
unfeasible  for  long  periods  of  real  time  or  large  computational  domains.  There  is  another  class 
of  models  which  have  traditionally  played  supporting  roles  in  fire  modeling.  This  class  is 
constituted  by  a  variety  of  network  models,  which  are  used  for  ventilation  systems  in  buildings 
or  fluid  flow  in  piping  networks.  These  models,  however,  lack  specific  physics  required  for  fire 
modeling. 

The  need  for  fire  modeling  occurs  through  the  “cradle-to-grave”  lifetime  of  a  ship,  which 
encompasses  the  phases  of  concept  design,  detailed  design  and  hazards  evaluation,  doctrine  and 
tactics  development,  ship  operation  and  maintenance,  and  crew  training.  In  the  concept  design 
phase,  fire  modeling  is  needed  to  evaluate  ship  designs  and  design  philosophies  in  order  to 
quickly  arrive  at  an  overall  concept  to  meet  required  performance  goals.  As  this  concept  is 
refined  into  a  detailed  design,  fire  modeling  is  continued  to  evaluate  ship  vulnerability  and 
recoverability  and  to  begin  the  process  of  defining  ship  operation.  Both  of  these  phases  could 
potentially  require  evaluating  the  fire  behavior  of  a  ship  and  its  system  for  hundreds  of  scenarios 
in  design  cycles  spanning  a  few  months  time.  Doctrine  and  tactics  planning  also  require 
evaluating  many  separate  fire  scenarios  in  combination  with  a  number  of  candidate  doctrines. 
During  the  operational  phase  of  a  ship’s  lifetime  fire  modeling  can  support  crew  training  through 
interactive  training  simulations.  Fire  modeling  can  also  aid  in  the  recovery  from  a  fire  event  by 
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providing  a  faster  than  real  time  predictive  capability  to  gauge  the  potential  effectiveness  of  a 
course  of  action. 

For  each  of  the  aforementioned  uses  of  fire  modeling,  speed  and  the  ability  to  model 
large  complexly  interconnected  spaces  with  ventilation,  detection,  and  suppression  with  a 
reasonable  degree  of  accuracy  are  performance  requirements  for  the  model.  Hand  calculations, 
while  fast,  have  limitations  in  applicability  and  large  uncertainties  in  their  results.  CFD 
computations  have  the  potential  to  be  accurate  to  any  desired  degree,  but  are  slow.  Zone  models 
balance  well  computational  time  and  accuracy;  however,  in  general  these  models  do  not  perform 
well  with  ventilation  systems  and  they  lack  the  ability  to  model  control  systems.  To  meet  the 
computational  speed  and  algorithm  requirements,  it  was  decided  to  develop  Fire  and  Smoke 
Simulator  (FSSIM)  as  a  network  fire  model. 

1.2  Network  Model  Concept 

A  network  model  represents  a  one-dimensional  abstraction  of  the  physical  world.  The 
prototype  being  modeled  is  represented  by  a  set  of  nodes  and  node  connections,  where  the  nodes 
contain  one  set  of  physical  variables  and  the  node  connections  represent  transfers  of  those 
variables  between  nodes.  Contrast  this  with  a  zone  model  which  uses  two  sets  variables  or  a  CFD 
model  which  uses  a  multitude.  A  network  representation  allows  for  maximum  physical  extent  of 
a  simulation  with  a  minimum  set  of  equations,  one  per  geometric  volume  of  interest.  Since  the 
number  of  locations  being  solved  for  is  kept  to  a  minimum,  a  network  model  also  has  the 
potential  for  the  faster  execution  speeds. 

There  are  currently  a  number  of  network  models  in  existence  that  are  used  for  safety 
related  simulations.  Some  of  these  include: 

•  TRAC  -  Transient  Reactor  Analysis  Code.  Used  for  simulating  normal  and  accident 
conditions  of  nuclear  power  plant  cooling  systems  [Spore,  et  al.,  2000]. 

•  RELAP  -  Reactor  Excursion  and  Leak  Analysis  Program.  Used  for  simulating 
normal  and  accident  conditions  of  nuclear  power  plant  cooling  systems  [USNRC, 
2001]. 

•  GOTHIC  -  Generation  of  Thermal  Hydraulic  Information  for  Containments.  Used 
for  simulating  normal  and  accident  conditions  of  nuclear  power  plant  containment 
buildings  and  reactor  buildings  [George,  et  al.,  2000a-b]. 

•  MELCOR  -  Used  for  simulating  normal  and  accident  conditions  of  nuclear  power 
plant  containment  buildings  and  reactor  buildings  including  aerosol  dispersion,  core 
melt/core  melt  interactions,  and  hydrogen  combustion  [Gauntt,  et  al.,  2000]. 

•  CONTAM  -  Used  for  simulating  the  dispersal  of  airborne  contaminants  through  a 
building  and  its  associated  HVAC  systems  [Dols,  Walton,  and  Denton,  2000]. 

•  COCOSYS  -  Containment  COdeSYStems  Used  for  simulating  accident  conditions 
in  light  water  nuclear  power  plant  containment  buildings  and  reactor  buildings 
[Allelein  et  al.,  1999;  Klein-Hepiing  et  ah,  2000]. 
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2.0  BACKGROUND 


FSSIM  is  a  network  fire  model  written  to  simulate  the  spread  of  fire  and  smoke  in  a  naval 
vessel.  However,  there  is  nothing  in  the  model  to  preclude  its  use  for  other  types  of 
compartmented  structures  such  as  a  building.  FSSIM  is  written  in  standard  Fortran  95  [Adams  et 
al.,  1997;  Lahey  Computer  Systems,  1998]  and  as  such  is  capable  of  being  compiled  on  any 
platform  for  which  a  Fortran  95  compiler  is  available.  In  FSSIM  each  compartment  in  a  structure 
is  represented  as  a  single  node  with  surfaces  (e.g.  bulkheads,  decks,  and  overheads)  and  vent 
openings  (e.g.  doors  and  hatches)  represented  as  node  connections.  There  is  no  practical  limit 
(2xl09  items  of  any  type)  to  the  size  of  an  FSSIM  simulation  other  than  available  computational 
resources  and  time. 

2.1  Documentation  Overview 

FSSIM  documentation  consists  of  three  volumes.  The  FSSIM  Theory  Manual  [Floyd  et  al. , 
2004a],  the  FSSIM  Users’  Manual  (this  document),  and  the  FSSIM  Validation  Manual  [Floyd  et 
al.,  2004b],  The  Theory  Manual  describes  the  equations  solved,  and  the  solution  algorithm  for 
the  heat  and  mass  transfer  along  with  the  equations  and  algorithms  for  FSSIM  sub-models.  The 
Users’  manual  discusses  creation  of  FSSIM  input  files  including  guidance  in  transforming  a 
prototype  structure  to  a  network  representation  and  explains  the  available  outputs.  Lastly,  the 
Validation  Manual  documents  FSSIM  performance  as  compared  to  experimental  data  and  other 
fire  models. 

There  is  a  companion  viewer  for  FSSIM  that  is  under  development  by  Mississippi  State 
University.  This  viewer  will  display  animated  results  of  FSSIM  calculations  in  a  three- 
dimensional  representation  of  the  ship  [Haupt  et  al.,  2004]. 

2.2  Features  of  FSSIM 

FSSIM  has  many  features,  including: 

•  One-dimensional  flow  model  including  friction  losses  and 
temperature-dependent  specific  heat. 

•  One-dimensional  multiple  layer,  temperature-dependent  heat  transfer. 

•  N-surface,  gray-gas  radiation  heat  transfer  including  radiation  streaming 
through  openings. 

•  Bidirectional  flow  through  horizontal  (hatches)  and  vertical  (doors)  flow 
connections. 

•  Combustion  product  species  tracking. 

•  Oxygen-  and  fuel-limited  combustion. 

•  Multiple  user-defined  fires  along  with  fire  spread  via 
compartment-to-compartment  heat  transfer. 
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HVAC  systems  including  ducts,  dampers,  and  fans  with  forward  and  reverse 
flow  losses  and  multiple  fan  models. 

Fire  detection  via  heat,  smoke,  and  fire  detection. 

Fire  spread  by  compartment-specific  criteria. 

Fire  suppression  via  sprinklers,  water  mist,  gaseous  agents,  aerosol  agents, 
and  foam. 

Fire  spread  prevention  via  boundary  cooling. 

Simple  control  systems  to  link  operation  of  equipment  to  sensors  or  times. 
Fast,  near-real-time  execution  speed. 

Can  be  run  in  a  standalone  configuration  or  as  a  callable  subprogram. 


2.3  Applications 

FSSIM  has  a  wide  range  of  applications,  including: 

•  Real-time  simulation  of  shipboard  conditions  during  a  fire. 

•  Damage  control. 

•  Event  tree  sequences. 

3.0  INSTALLING  AND  RUNNING  FSSIM 

3.1  Computer  Hardware  and  Software  Requirements 

FSSIM  requires  a  moderately  fast  central  processing  unit  (CPU)  with  a  sufficient  amount 
of  random  access  memory  (RAM)  to  accommodate  the  complexity  of  the  particular  simulation. 
FSSIM  will  alert  the  user  if  there  is  not  a  sufficient  amount  of  memory  for  an  attempted 
simulation.  For  a  Windows  based  personal  computer  (PC),  a  processor  that  is  as  fast  as  a  400 
Mhz  Pentium  II  should  be  minimally  sufficient  to  execute  basic  FSSIM  configurations.  More 
complex  simulations  should  be  executed  on  a  faster  system;  otherwise  real-time  simulations  may 
not  be  attained.  For  a  UNIX-based  station,  the  processor  should  be  at  least  as  fast  as  the  PC 
platform  specifications. 

3.2  Running  FSSIM 

3.2.1  Installing  FSSIM 

FSSIM  may  be  installed  from  a  compact  dick  (CD)  or  it  may  be  downloaded  from  a 
designated  site.  If  FSSIM  is  being  installed  from  a  CD,  the  following  applies: 

1 .  Create  a  sub-directory  where  FSSIM  will  be  executed. 

2.  Copy  all  files  from  the  CD  into  the  sub-directory  created.  Files  should  include: 
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•  fssimi  .  exe  -  the  executable  program. 

•  materials .  txt  -  material  properties  file  containing  information  for 
common  materials. 

•  usetypes .  txt  —  file  containing  fuel-loading  definitions  of  common 
compartment  classifications. 

•  fuels .  txt  -  fuel  properties  file  containing  information  for  common  fuels. 

•  samplel .  in,  sample2  .  in,  sample3 .  in,  and  sample4  .  in  —  sample  input 
files. 

3.  Run  the  sample  input  files  to  verify  proper  installation  (see  Section  3.2.3). 

If  FSSIM  is  being  downloaded  from  an  ftp  site  or  equivalent,  the  following  applies: 

1 .  Create  a  sub-directory  where  FSSIM  will  be  executed. 

2.  Download  the  compressed  file  fssimi  .  zip . 

3.  Extract  the  files  from  fssimi  .  zip  into  the  sub-directory  created  in  step  1. 

4.  Run  sample  input  files  to  verify  proper  installation  (see  Section  3.2.3). 


3.2.2  Creating  the  Input  File 

FSSIM  requires  an  input  file  containing  details  regarding  the  configuration  simulated. 

The  input  file  is  an  ASCII  file  that  may  be  edited  using  any  suitable  text  editor.  Section  4 
contains  details  on  preparing  the  input  file  and  Section  6  provides  several  input  sample  data  files. 


3.2.3  Starting  the  FSSIM  Calculation 

FSSIM  is  a  command  line  executable  program.  The  program  must  be  run  from  a  disk 
operating  system  (DOS)  or  UNIX  shell. 

•  Microsoft  Windows  based  PC 

Open  a  Command  Prompt  window,  change  directories  to  where  fssim.exe  and  the 
necessary  input  files  are  located,  then  run  the  model  as  follows: 

fssim  inputfile.in 

where  input  file .  in  is  any  valid  input  file.  The  input  file  will  be  read  and  the 
simulation  begun.  The  current  time  and  other  information  will  be  printed  to  the  screen 
so  that  the  user  may  monitor  the  progress  of  a  simulation.  If  the  input  file  is  not  valid, 
the  program  will  be  aborted  and  an  error  will  be  generated.  Note  that  the  executable 
and  the  input  files  are  not  case  sensitive  in  a  DOS  shell. 
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The  information  printed  to  the  screen  during  a  simulation  may  instead  be  saved  in  a 
file  by  running  FSSIM  as  follows: 

fssim  inputfile.in  >  output. out  (overwrites  output . out) 

or 

fssim  inputfile.in  »  output . out  (appends  to  output , out) 
where  output .  out  is  any  output  file  name. 

•  UNIX  platform 

Open  a  UNIX  shell,  change  directories  to  where  fssim.exe  and  the  necessary  input 
files  are  located,  then  run  the  model  as  follows: 

fssim  <  inputfile.in 

where  inputf  ile .  in  is  any  valid  input  file.  The  input  file  will  be  read  and  the 
simulation  begun.  The  current  time  and  other  information  will  be  printed  to  the  screen 
so  that  the  user  may  monitor  the  progress  of  a  simulation.  If  the  input  file  is  not  valid, 
the  program  will  be  aborted  and  an  error  will  be  generated.  Note  that  the  executable 
and  the  input  files  are  case  sensitive  in  a  UNIX  shell. 

The  information  printed  to  the  screen  during  a  simulation  may  instead  be  saved  in  a 
file  by  running  FSSIM  as  follows: 

fssim  <  inputfile.in  >  output. out  &  (overwrites  output . out) 
where  output .  out  is  any  output  file  name. 


3.2.4  Common  Runtime  Errors 

FSSIM  may  fail  to  run  for  a  variety  of  reasons.  Common  problems  are  summarized  in  the 
table  below. 
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Table  3-1.  Common  FSSIM  errors  and  Possible  User  Actions 


FSSIM  Error 

Possible  Cause 

User  Actions 

FSSIM  does  not  find  one  or  more 
compartments,  junctions,  surfaces, 
compositions,  materials,  ducts,  duct 
nodes,  fans,  suppressors,  and/or 
detectors. 

The  &CTRL  line  does  not  match 
number  of  desired  inputs. 

•  Make  sure  that  the  parameters 
on  the  &CTRL  line  match  the 
number  of  lines  that  should  be 
read  for  each  item. 

•  Make  sure  that  all  NAMELIST 
groups  are  correctly  spelled. 

•  Make  sure  that  that  NAMELIST 
groups  with  parameters  set  in 
auxiliary  input  files  are  properly 
set  up  for  file-read. 

Program  aborts  unexpectedly  or 
inexplicably. 

Input  file  format/setup  error. 

Examples  include  a  wrong  variable, 
a  NAMELIST  syntax  error,  and  an 
undefined  parameter  that  is  required. 

Verify  that  the  input  file  is  in  the 
correct  format  and  that  all  required 
parameters  are  included  on  each 
NAMELIST  line. 

Invalid  id  parameter  value.  An  id 
parameter  value  is  a  unique  integer 
that  identifies  various  FSSIM 
entities,  such  as  a  compartment,  an 
opening,  or  a  node.  These  values  are 
used  to  locate  openings,  nodes, 
ducts,  fans,  etc .  relative  to  each 
other.  FSSIM  does  not  check  for 
valid  id  parameter  values.  If  an 
invalid  id  is  used,  FSSIM  may  run, 
partially  run,  or  abort  immediately 
with  or  without  indicating  that  there 
is  an  invalid  id  parameter  value. 

Verify  that  all  id  parameters  are 
valid,  both  as  defined  and  as 
supplied  for  parameter  arguments. 

Improper  physical  geometry.  This 
condition  may  arise  when  an 
incorrect  id  parameter  value  is 
supplied  as  argument,  such  that  an 
opening  or  a  duct  is  connected  to  the 
wrong  compartment  or  node. 

Verify  the  input  configuration. 

FSSIM  Fails  to  converge. 

See  “Program  aborts  unexpectedly 
or  inexplicably”. 

See  “Program  aborts  unexpectedly 
or  inexplicably”. 

Convergence  parameters  must  be 
adjusted  for  configuration  modeled. 

Increase  the  number  of  inner  or 
|  outer  timestep  loops,  decrease 

1  minimum  timestep  size,  simplify 
mechanical  ventilation  system, 
and/or  add  natural  ventilation  to 
spaces  that  are  well  sealed. 

FSSIM  will  not  find  an  input  file. 

Incorrect  input  file  or  &FNAM 
specified  file  is  not  valid  or  not 
located  in  current  directory. 

Verify  that  the  input  file  is  in  current 
directory.  Verify  that  any  files  set 
using  the  &FNAM  NAMELIST 
group  and  used  as  input  vehicles  for 
die  particular  simulation  are  located 
in  the  current  directory. 
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FSSIM  Error 

Possible  Cause 

User  Actions 

The  results  do  not  make  physical 
sense. 

Incorrect  values  used.  An  example  is 
a  heat  release  rate  provided  in  W 
where  kW  is  required  for  input. 

Verify  that  the  values  provided  as 
parameter  arguments  have  the 
proper  units.  Verify  that  the  physical 
values  and  quantities  supplied  are 
the  values  intended. 

Insufficient  memory. 

Configuration  is  too  large  for 
computer. 

Reduce  the  size  of  the  problem  or 
use  a  computer  with  a  greater 
amount  of  RAM. 

FSSIM  fails  to  run  and  none  of  the 
user  actions  above  correct  the 
problem. 

Computer  code  error. 

Contact  the  primary  author: 

Jason  Floyd 

3660  Commerce  Drive,  Suite  817 
Baltimore,  MD  21227-1652 

iflovdfShaifire.com 

4.0  FSSIM  INPUT  FILE  STRUCTURE 

FSSIM  is  ran  using  a  valid  ASCII  input  file  that  contains  a  description  of  the  input 
geometry,  various  initial  and  subsequent  conditions,  and  other  user-defined  information.  FSSIM 
has  some  error  trapping;  however  there  is  a  significant  reliance  on  the  user  to  adequately  prepare 
the  input  file  for  the  configuration  desired,  both  in  terms  of  physical  representation  and  input  file 
integrity.  This  section  summarizes  the  structure  of  the  FSSIM  input  files.  Specifically,  the 
following  are  addressed: 

•  Key  elements  of  FSSIM  input  file  as  they  relate  to  physical  configurations. 

•  Input  parameters/variables  that  are  accessible  to  the  user  via  NAMELIST  groups 
and  are  used  to  build  a  simulation. 

•  The  type  of  input  variable  or  parameter  (integer,  real  number,  text,  or  logical). 

•  Appropriate  variable/parameter  input  ranges,  where  applicable. 

•  The  minimum  or  required  input  variables/parameters. 

•  The  default  conditions. 

NAMELIST  groups  allow  the  user  to  access  a  variety  of  model  parameters  or  variables, 
required  or  optional.  These  parameters  and  variables  are  identified  for  each  NAMELIST  group 
available  in  the  FSSIM  model  in  the  subsequent  sections.  The  input  format  is  noted  as  a  real 
number,  an  integer,  text,  or  a  logical  statement.  The  input  variable  or  parameter  maybe  a  single 
element  or  an  array  of  elements.  Note  that  real  numbers  maybe  used  as  single  precision  or 
double  precision  variables  within  the  model.  Because  the  input  file  format  is  transparent  to  this,  a 
distinction  between  the  two  variable  types  is  not  necessary  in  this  guide. 

Input  data  ranges  or  restrictions  that  are  intuitively  obvious  are  not  listed.  This  includes 
the  restrictions  that  time,  timesteps,  convergence  criteria,  area,  volume,  length,  width,  pressure. 
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density,  heat  capacity,  thermal  conductivity,  and  number  of  items  are  positive  values  greater  than 
zero  (unless  specifically  allowed  to  be  zero);  minimum  values  are  less  than  or  equal  to  maximum 
values;  identifications  for  unique  elements  are  unique;  and  related  restrictions. 

In  the  following  subsections  each  NAMELIST  group  will  be  described.  Not  all 
NAMELIST  groups  will  be  required  for  a  given  simulation.  In  addition,  for  those  NAMELIST 
groups  needed,  not  all  variables  in  the  group  will  need  to  be  specified.  Each  group  will  contain 
some  variables  that  the  user  must  set  for  a  successful  simulation.  NAMELIST  may  also  contain 
subsets  of  variables  of  which  only  one  may  be  set  (i.e.,  the  variables/parameters  are  mutually 
exclusive).  Many  parameters/variables  that  are  necessary  for  performing  a  simulation  have 
default  values,  and  thus,  do  not  have  to  be  set  by  the  user  unless  a  value  other  than  the  default  is 
desired.  Required  parameters,  default  values,  mutually  exclusive  parameters,  and  other  such 
conditions  are  identified  for  each  NAMELIST  group  in  the  subsequent  sections. 

4.1  Overview 

4.1.1  NAMELIST  Input  Format 

The  user  builds  a  simulation  model  in  FSSIM  using  NAMELIST  parameters  groups.  A 
NAMELIST  parameter  group  is  a  means  of  setting  variables  and  parameters  within  the  model 
that  are  associated  with  a  particular  group,  use,  or  function.  The  format  is  as  follows: 


&NLID  variablel=x,  variable2=.  FALSE. ,  variable3='  string' ,  ...  / 


where  snlid  is  the  NAMELIST  parameter  group  identifier  and  all  variables  that  are  set  are 
located  between  the  group  identifier  and  a  7’  symbol  and  are  separated  from  each  other  by  a 
comma.  A  NAMELIST  line  is  defined  as  all  of  the  information  located  between  the  c&’  and  */’ 
symbols  and  a  specific  NAMELIST  line  is  defined  as  all  of  the  information  between  the  group 
identifier  and  the  ‘  /  ’  symbol. 

FSSIM  uses  five  character  NAMELIST  group  identifiers.  The  group  identifier  must  be 
located  in  the  first  five  columns  of  the  input  file  and  variables  that  are  set  in  a  NAMELIST  group 
must  be  separated  from  each  other  by  a  comma.  Other  than  that,  there  are  no  formatting 
restrictions  on  the  input  file  and  NAMELIST  format  structure.  There  are  also  no  restrictions  on 
the  order  in  which  the  variables  appear  within  the  NAMELIST  group  or  the  order  in  which  the 
NAMELIST  groups  appear  within  an  input  file.  Thus, 


&NLID  variable2=.  FALSE . ,  variable3='  string'  ,  variablel=x/  ...  / 


is  equivalent  to  the  snlid  NAMELIST  parameter  group  example  provided  above  and  it  may 
appear  anywhere  within  the  input  file. 

Text  located  on  the  same  input  file  line  located  to  the  right  of  a  ‘ !  ’  symbol  is  treated  as  a 
user  comment.  Thus,  NAMELIST  lines  may  be  converted  to  user  comments  by  placing  the  ‘ !  ’ 
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symbol  in  column  1.  In  addition,  text  that  appears  outside  the  &nlid  /  construct  is  treated  as  a 
user  comment.  The  following  are  examples  of  user  input  file  comments: 


i&NLID  variable2=. FALSE, ,  variable3=' string' ,  variablel=x,  ...  /The  &NLID  line 
and  this  text  are  user  comments 

&NLID  variable2=. FALSE. ,  variable3='  string' ,  variablel=x,  ...  /Only  this  text 
is  a  user  comment 


The  input  is  not  case  sensitive;  however,  it  is  conventional  to  use  upper  case  for  the 
NAMELIST  parameter  identifier  and  whatever  case  or  combination  of  cases  are  used  by  the 
variable  or  parameter  within  the  NAMELIST  group.  The  only  exception  to  this  are  text  string 
variables,  which  are  case  sensitive.  This  is  of  particular  importance  when  naming  items  that  must 
be  retrieved  from  files.  The  names  must  match  exactly  for  a  successful  retrieval. 

There  may  be  one  or  more  instances  of  a  particular  NAMELIST  parameter  group 
contained  in  an  input  file.  In  such  case,  a  unique  identifier  is  normally  required  to  be  supplied  by 
the  user.  Multiple  instances  of  a  NAMELIST  parameter  group  are  necessary  when  there  is  more 
than  one  compartment;  multiple  surfaces  or  materials  within  a  compartment;  more  than  one  duct, 
fan,  or  duct  node;  multiple  fires;  multiple  openings;  etc.  This  is  described  in  greater  detail  in  the 
sections  that  follow.  If  there  are  more  instances  of  a  particular  NAMELIST  parameter  group 
defined  than  are  necessaiy  or  required,  only  those  up  to  the  required  number  are  actually  read  as 
input  Thus,  if  ten  instances  of  a  particular  NAMELIST  group  are  required  and  twelve  are 
provided,  the  first  ten  are  read  as  input  and  the  remaining  two  are  ignored. 


4. 1 .2  Variable  and  Parameter  Input  Format 

There  are  format  requirements  associated  with  each  type  of  variable  input.  They  are  as  follows: 

*  Real  numbers  are  entered  using  a  decimal: 
real_variable__example=l  ,75 


•  Integers  are  entered  without  using  a  decimal: 
integer_variable_example=5 


*  Logical  variables  are  entered  using  true  or  false  statements: 
logical_variable_examplel= . FALSE . 
logical_variable_example2= . TRUE . 

•  Text  strings  are  entered  using  single  quotes  and  should  not  contain  commas  to  avoid 
potential  output  formatting  errors: 
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text_string_examplel=' example  string  -  this  string  is  case  sensitive' 
text_string_example2=' Example  string  -  this  string  is  case  sensitive' 


Real  number  variables  that  are  entered  as  integers  will  be  converted  to  single  or  double  precision 
variables  in  the  model.  Conversely,  integer  variables  that  are  entered  as  real  numbers  are 
truncated  and  converted  to  an  integer  in  the  model.  Logical  states  may  be  expressed  as  either 
.true,  or  .false,  or  equivalently  as  .true,  or  .false.,  though  the  former  is  the 
conventional  input  style. 

An  integer  normally  has  a  value  between  -215  and  +215  (-32,767  and  +32,767). 

Exceptions  to  this  are  short  integers,  which  must  have  a  value  between  -27  and  +27  (-127  and 
+127),  and  long  integers,  which  must  have  a  value  between  -231  and  +231  (-2,147,483,647  and 
+2,147,483,647).  Short  and  long  integers  are  used  for  some  types  of  data  and  are  identified 
where  applicable.  Integer  values  entered  that  exceed  the  bounds  for  the  particular  integer  type 
will  produce  erratic  results.  Text  strings  that  exceed  the  maximum  length  are  truncated  to  the 
maximum  length.  Note  that  text  strings  are  case  sensitive,  such  that  the  two  string  examples  cited 
above  are  not  equal  because  the  ‘e’  is  capitalized  in  the  second  example. 

Some  NAMELIST  group  parameters  require  more  than  one  input  value.  These  are 
referred  to  in  this  document  as  parameter  arrays  and  they  are  input  on  a  NAMELIST  line  using 
comma  separators  between  the  array  elements.  The  following  are  examples  of  parameter  arrays: 


•  Integer  parameter  array  with  two  elements : 
example_integer_array=l, -1 

•  Real  number  parameter  array  with  six  elements: 
example_real_number_array=1.53,  2.66,  2.66,  1.0,  1.44,  1.354 


There  are  no  logical  or  text  string  variable  arrays  used  by  the  FSSIM  NAMELIST  groups.  The 
input  order  for  parameter  arrays  may  or  may  not  be  significant.  Descriptions  of  specific 
NAMELIST  group  parameter  arrays  will  indicate  whether  or  not  the  array  order  is  of 
significance. 


4.1.3  Model  Elements 

FSSIM  is  used  to  simulate  a  fire  and  determine  the  resulting  compartment  conditions 
within  one  or  more  compartments.  The  configuration  and  scenario  simulated  are  entered  using  up 
to  eighteen  different  NAMELIST  parameter  groups.  Each  parameter  group  sets  related 
parameters,  such  as  the  dimensions  of  a  compartment  or  the  location  and  size  of  openings  in  a 
compartment.  Figure  4-1  depicts  a  sample  compartment  and  identifies  the  parameter  groups  used 
to  set  the  parameters  needed  to  simulate  a  fire  scenario.  Table  4-1  summarizes  the  parameter 
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Compartment  Information  I  Genera!  Simulation  and  Data  Control 


groups  and  the  types  of  parameters  they  set.  Sections  4-3  through  4-5  summarize  the  syntax  of 
each  NAMELIST  parameter  group  shown  in  Table  4-1. 


Table  4-1.  FSSIM  NAMELIST  Parameter  Groups 


NAMELIST  Parameter 

Group 

Requirements/ 

Use 

Type  of  parameters  set  by  the  NAMELIST  Parameter 

Group 

General  Simulation  and  Data  Control 

&EXEC 

Required,  One  and  only 
one  &EXEC  line  must  be 
provided. 

*  Time  and  time-related  parameters. 

*  Convergence  control. 

*  Ambient  conditions. 

*  Output  format. 

&CNTR 

Required.  One  and  only 

one  &CNTR  line  must  be 
provided. 

•  The  number  of  specific  types  of  NAMELIST 
parameter  groups  to  be  input. 

&OUTP 

Optional.  If  used,  only 
one  &OUTP  line  is 
provided 

*  Specifies  outputs  to  be  written  if  the  variable 
outt ype  on  &EXEC  is  .  FALSE  * 

Optional,  If  used,  only  one 
&FNAM  line  is  provided. 

*  Names  of  files  used  to  provide  input  data  for 
specific  NAMELIST  parameter  groups. 

&CURV 

Optional.  One  or  more 
&CURV  lines  must  be 
provided  as  referenced  by 
other  NAMELIST 
parameter  groups. 

*  Mechanism  to  enter  tabular  data,  such  as  material 
properties  of  heat  release  rate. 

&CTRL 

Optional.  One  or  more 
&CTRL  lines  must  be 
provided  if  referenced  by 
another  NAMELIST 
group. 

•  Mechanism  to  turn  objects  on  or  off,  such  as  a 

fan  or  a  suppression  system.  Also  be  used  to 
open  or  close  a  door  or  a  fire  damper. 

Compartment  Information 

&COMP 

mm 

•  Compartment  dimensions  and  volume. 

•  Initial  compartment  conditions. 

•  Type  of  combustible  fuel  load  in  a  compartment. 

&  JUNC 

Optional.  One  &  JUNG  line 
is  provided  for  each 
opening,  if  any. 

*  Opening  dimensions,  location,  and  orientation. 

*  Bidirectional  or  one-way  flow. 

*  Type  of  opening  control. 

*  Opening  flow  loss  coefficient  and  inertial  length. 

&SURF 

Required.  Minimum  of  six 
for  each  unique  &COMP 
line. 

*  Type,  location,  and  orientation  of  boundary 
surfaces. 

*  Area  and  composition  of  boundary  surfaces. 

&CMPN 

Required,  unless  all 
surfaces  are  fake.  One 
&CMPN  line  is  required  for 
each  unique  composition 
identified  by  the  &SURF 
lines. 

•  The  thickness  and  type  of  material  composing  a 
solid. 

*  The  number,  thickness,  and  type  of  material 
composing  each  individual  layer  in  composite 
solid. 

&MTRL 

Required,  unless  all 
surfaces  are  fake.  One 
&MTRL  line  is  required  for 
each  unique  material 
identified  by  the  &CMPN 
parameter  groups. 

•  The  thermal  conductivity,  heat  capacity,  density, 
and  emissivisity  of  a  material. 
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Detection/Fire  Suppression  I  Source  Fire  and  Fuel  Load  I  Mechanical  Ventilation 


NAMELIST  Parameter 
Grou 


Requirements/ 

Use 


Optional.  One  &RDCT  line 
is  provided  for  each  duct 
segment  modeled,  if  any. 

Optional.  One  &RNOD  line 
is  provided  for  each  duct 
node  location  modeled,  if 

any. _ 

Optional.  One  &RFAN  line 
is  provided  for  each  fan 
modeled,  if  any. _ 


Optional.  One  &FIRE  line 
is  provided  per 
compartment  if  a  fire  is 
directly  specified  within 
one  or  more 
compartments. 

Optional.  One  &USES  line 
is  provided  per 
compartment  not 
containing  a  fire  (&FIRE) 
if  an  initially  un-ignited 
fuel  load  is  simulated. 


Optional.  &FUEL  lines 
must  be  provided  as 
referenced  by  the  &USES 
Iine(s). 


Optional.  One  &DTCT  line 
for  each  type  of  detection 
system  in  each 
compartment  where 
detection  is  simulated,  if 

any. _ _ 

Optional.  One  &SUPR  line 
for  each  type  of  detection 
system  in  each 
compartment  where 
suppression  is  simulated,  if 
any. 


Type  of  parameters  set  by  the  NAMELIST  Parameter 
Grour 


•  The  length  and  cross  section  of  duct  segments. 

•  Flow  pressure  loss  characteristics  through  ducts. 

•  Duct  damper  control. 

•  Type  of  fan  in  duct. 

•  Constant  mass/or  volume  flow  rates  through 

_ ducts. _ 

•  Location  and  elevation  of  node. 

•  Ducts  connected  to  node. 

•  Flow  pressure  loss  characteristics  at  node. 

•  Type  of  fan. 

•  Fan  performance  characteristics. 

•  Flow  pressure  loss  characteristics  through  fan. 


•  Fire  heat  release  rate  profile. 

•  Fire  location. 

•  Radiant  fraction  of  fire  heat  release  rate. 

•  Heat  of  combustion  and  heat  of  vaporization  of 
burning  material. 

•  Yield  of  various  gaseous  and  solid  combustion 
products  of  burning  material. 


Location  of  fuel  load. 
Type  of  fuel  load. 


Growth  rate  of  fire  when  fuel  material  ignited. 
Radiation  fraction  of  heat  release  rate  when  fuel 
material  ignited. 

Heat  of  combustion  and  heat  of  vaporization  of 
fuel  material. 

Fuel  material  ignition  criteria. 

Yield  of  various  gaseous  and  solid  combustion 
products  of  fuel  material. 


Type  of  detector. 

Location  of  detector. 
Detector  actuation  criteria. 


Type  of  suppression  system. 

Location  of  suppression  system. 

Suppression  system  performance  characteristics. 


4.2.1  Execution  Control  (&exec) 


The  &exec  group  allows  the  user  to  control  various  aspects  of  an  FSSIM  simulation 
including  initial  ambient  conditions,  the  solution  time,  convergence  criteria,  and  the  type  of 
output.  The  input  file  should  contain  only  one  &exec  parameter  group.  There  is  one  parameter 
that  must  be  set  on  the  &exec  line: 

REQUIRED  PARAMETER 


(Text  string)  Text  string  used  to  generate  output  files.  There  is  no  default 
and  the  string  must  be  40  characters  or  less. 


The  remaining  sexec  parameters  are  optional: 


OPTIONAL  PARAMETERS 


dt 

(Real  number)  The  first  timestep  (s).  Default  value  is  0.01  seconds. 

t_max 

(Real  number)  The  simulation  end  time  (s).  Default  value  is  100  seconds. 

dt_max 

(Real  number)  The  maximum  timestep  size  (s).  Default  value  is  1.0 
second. 

dt  min 

(Real  number)  The  minimum  timestep  size  (s).  Default  value  is  0.0001 
seconds. 

dt__dump 

(Real  number)  The  minimum  time  increment  that  output  is  generated  (s). 
The  default  value  is  0.0  seconds,  which  corresponds  to  every  timestep. 

humidity 

(Real  number)  The  relative  humidity  in  %.  The  default  is  70. 

ninnerjnax 

(Short  integer)  The  maximum  number  of  inner  loop  (velocity  solution) 
iterations  in  the  network  solver.  Default  is  10. 

nout  max 

(Short  integer)  The  maximum  number  of  outer  loop  (pressure  solution) 
iterations  in  the  network  solver.  Default  is  5. 

temp_criteria 

(Real  number)  The  temperature  convergence  criteria.  Default  is  0.1. 

pres__criteria 

(Real  number)  The  pressure  convergence  criteria.  Default  is  0.01. 

courant_criteria 

(Real  number)  The  Courant  convergence  criteria.  Default  is  0.5. 

tamb 

(Real  number)  The  ambient  temperature  (K).  Default  is  298. 1 5K. 

pa  mb 

(Real  number)  The  absolute  ambient  pressure  (Pa).  Default  value  is 

101,325  Pa. 

o2amb 

(Real  number)  The  ambient  oxygen  mass  fraction.  Default  is  0.23. 

o211 

(Real  number)  The  lower  limit  oxygen  mass  fraction  that  will  support 
combustion  (kg  02/kg  gas).  The  default  value  is  0.1  kg  02/kg  air. 
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outtype 

(Logical)  .  false  .  for  ASCII  text  output  and  .  true  .  for  Mississippi 

State  University  (MSU)  graphical  user  interface  (GUI)  output  format. 

The  default  state  is  .  true  . 

The  default  convergence  criteria  and  timestep  parameters  should  be  adequate  for  many 
simulations.  If  there  is  concern  regarding  solution  convergence,  the  temperature,  pressure,  and 
Courant  convergence  criteria  or  the  initial  and  maximum  timestep  could  be  reduced  and  the 
results  compared  to  verify  convergence.  Likewise,  the  time  required  to  obtain  a  solution  may  be 
decreased  by  increasing  the  tolerances  or  the  initial  timestep  size.  There  may  be  a  loss  of 
accuracy  if  either  is  decreased.  If  the  model  fails  to  converge  for  a  particular  simulation, 
increasing  the  number  of  inner  (velocity  and  fire  size)  and  outer  (pressure  and  courant  stability) 
iterations  allowed  or  decreasing  the  convergence  criteria  and  initial  timestep  may  improve  the 
model  stability  and/or  the  prospects  for  a  converged  solution.  The  user  is  referred  to  the  Theory 
Manual  for  a  discussion  of  the  iteration  loops. 

An  example  of  a  sexec  NAMELIST  parameter  group  line  is  as  follows: 


&EXEC  t_max=600 . 0,  dt_dump=1.0,  tamb=302.5,  outtype=. FALSE. ,  fileid=' sample' , 
pamb=9.8E4  / 


In  this  case,  a  solution  time  of  600  seconds  is  declared,  the  ambient  temperature  and 
pressure  are  set  to  302.5  K  and  98,000  Pa,  respectively,  the  output  files  will  begin  with  the  text 
string  ‘sample5,  and  the  output  will  be  in  ASCII  text  format.  Output  will  be  generated  at 
intervals  no  less  than  1  second.  The  default  convergence  criteria  and  timestep  parameters  are 
used. 


4.2.2  NAMELIST  Group  Input  Counters  (&cntr) 

The  &cntr  group  is  used  to  specify  the  number  of  specific  NAMELIST  groups  to  be  read 
where  the  number  may  vary  (i.e.,  compartment  and  compartment  related  information,  openings, 
ventilation  and  suppression  groups,  fire  and  fire  related  information).  There  is  one  parameter  that 
must  be  set  on  the  &cntr  line; 


REQUIRED  PARAMETERS 

(Integer)  Number  of  compartment  (&comp)  NAMELIST  groups  to  be 
read.  At  least  one  compartment  must  be  specified.  The  default  is  zero  and 
if  used  will  generate  an  error. 


The  remaining  &cntr  group  parameters  are  optional.  However,  they  may  be  required  to 
be  set  if  any  of  the  NAMELIST  groups  linked  to  the  scntr  group  variable  are  used.  The 
following  is  a  complete  list  of  these  parameters: 
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OPTIONAL  PARAMETERS 


nf  ire 

(Integer)  Number  of  initial  source  fires  (&fire)  NAMELIST  groups  to  be 
read.  At  least  one  initial  fire  must  be  specified,  though  it  may  have  a  zero 
heat  release  rate  or  power.  The  default  is  zero. 

nj  unc 

(Integer)  Number  of  junctions/openings  (& junc)  NAMELIST  groups  to 
be  read.  The  default  is  zero. 

nnotbidir 

(Integer)  Number  of  junctions/openings  (&  junc)  NAMELIST  groups  that 
do  not  have  bidirectional  flow.  The  value  must  be  equal  to  or  less  than 
njunc.  The  default  is  zero. 

nfuel 

(Integer)  Number  of  fuel  (sfuel)  NAMELIST  groups  to  be  read.  The 
default  value  is  zero. 

ntype 

(Integer)  Number  of  compartment  use  type  (suses)  NAMELIST  groups 
to  be  read.  The  default  value  is  zero. 

nsurf 

(Long  integer)  Number  of  duct  node  (&surf)  NAMELIST  groups  to  be 
read.  The  default  value  is  zero. 

ncmpn 

(Integer)  Number  of  composition  (&cmpn)  NAMELIST  groups  to  be  read. 
The  default  is  zero. 

nmat 

(Integer)  Number  of  material  (smtrl)  NAMELIST  groups  to  be  read. 

The  default  is  zero. 

nduct 

(Integer)  Number  of  duct  (&rdct)  NAMELIST  groups  to  be  read.  The 
default  is  zero. 

nnode 

(Integer)  Number  of  duct  (srnod)  NAMELIST  groups  to  be  read.  The 
default  is  zero. 

nf  an 

(Integer)  Number  of  fan  (&rfan)  NAMELIST  groups  to  be  read.  The 
default  is  zero. 

ncontrol 

(Integer)  Number  of  control  (&ctrl)  NAMELIST  groups  to  be  read.  The 
default  is  zero. 

ndetect 

(Integer)  Number  of  detector  (&dtct)  NAMELIST  groups  to  be  read.  The 
default  is  zero. 

nsuppress 

(Integer)  Number  of  suppression  (&supr)  NAMELIST  groups  to  be  read. 
The  default  is  zero. 

17 


The  number  of  openings,  or  junctions,  that  are  not  bidirectional  (i.e.,  the  solver  will  allow 
for  flow  in  each  direction  across  the  same  opening)  must  be  specified  separately  from  the  total 
number  of  openings.  The  model  uses  the  concept  of  zones  to  subdivide  a  network  into 
independent  volumes.  Spaces  that  are  connected  via  junctions  or  ductwork  share  zones.  The  zone 
concept  is  used  to  increase  the  efficiency  of  the  solver  and  the  model  currently  calculates  the 
number  of  zones  in  a  network.  Zoning  is  established  automatically  during  the  initialization 
process. 

Note  that  the  number  of  compartments,  openings,  detectors,  etc.  within  a  network  model 
are  limited  by  the  size  of  an  integer  (215),  a  sufficiently  large  number  for  nearly  any  practical 
system.  An  exception  to  this  is  the  number  of  surfaces,  which  is  limited  by  the  size  of  a  long 
integer  (231). 

An  example  of  a  scntr  NAMELIST  parameter  group  line  is  as  follows: 


SCNTR  ncomp=5,  nfire=l,  njunc=15,  nnotbidir=3,  nsurf=35,  ncmpn=2,  nmat=2, 
nduct=4,  nnode=4,  nfan=l,  ncontrol=l,  ndetect=l  / 


The  above  scntr  line  instructs  FSSIM  to  read  five  scomp,  one  spire,  fifteen  sjunc, 
thirty-five  ssurf,  two  scmpn,  two  &mtrl,  four  &rdct,  four  srnod,  one  srfan,  one  sctrl,  and 
one  &dtct  NAMELIST  groups.  In  addition,  twelve  of  the  openings  specified  by  the  &  junc 
NAMELIST  group  have  bidirectional  flow. 


4.2.3  Output  Data  Options  (soutp) 

The  soutp  group  Is  used  to  specify  the  output  files  written  when  outtype= .  true  .  on 
sexec.  All  parameters  associated  with  this  NAMELIST  group  are  optional.  The  output  file 
names  and  file  formats  are  discussed  in  Section  5.1.  The  following  is  a  complete  list  of 
parameters  set  by  the  soutp  group: 


OPTIONAL  PARAMETERS 


comp temp 

(Logical)  Write  compartment  temperature  output  file.  The  default  is 
.TRUE. 

comppres 

(Logical)  Write  compartment  pressure  output  file.  The  default  is  .TRUE . 

compq 

(Logical)  Write  compartment  radiation  absorption/emission  to  gas  and 
net  heat  flux  to  wall  output  file.  The  default  is  .FALSE . 

eompgas 

(Logical)  Write  compartment  gas  output  file.  The  default  is  .true  .  The 
constituent  gasses  are  indicated  at  the  end  of  this  table. 

fire 

(Logical)  Write  fire  theoretical  heat  release  rate  output  file  and 
compartment  actual  heat  release  rate  output  file.  The  default  is  .TRUE . 
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nodetemp 


(Logical)  Write  duct  node  temperature  output  file.  The  default  is  .TRUE . 


nodepres  (Logical)  Write  duct  node  pressure  output  file.  The  default  is  .TRUE . 


nodegas 

(Logical)  Write  duct  node  gas  mass  fraction  output  file.  The  default  is 
.true  .  The  constituent  gasses  are  indicated  at  the  end  of  this  table. 

juncvel 

(Logical)  Write  junction  velocity  output  file.  The  default  is  .false  . 

juncarea 

(Logical)  Write  junction  flow  area  output  file.  The  default  is  .FALSE . 

juncf low 

(Logical)  Write  junction  mass  flow  rate  output  file.  The  default  is  .TRUE . 

ductvel 

(Logical)  Write  duct  velocity  output  file.  The  default  is  .FALSE . 

ductarea 

(Logical)  Write  duct  flow  area  output  file.  The  default  is  .FALSE . 

duct flow 

(Logical)  Write  duct  mass  flow  rate  output  file.  The  default  is  .true  . 

surftemp 

(Logical)  Write  surface  temperature  output  file.  The  default  is  .true  . 

surf flux 

(Logical)  Write  surface  heat  flux  output  file.  The  default  is  .FALSE . 

detect 

(Logical)  Write  detector  alarm  status  output  file.  The  default  is  .TRUE . 

suppress 

(Logical)  Write  suppression  system  activation  status  output  file.  The 
default  is  .TRUE . 

gaso2 

(Logical)  Write  oxygen  mass  fraction  to  compartment/duct  node  gas 
mass  fraction  output  file.  The  default  is  .TRUE . 

gasco2 

(Logical)  Write  carbon  dioxide  mass  fraction  to  compartment/duct  node 
gas  mass  fraction  output  file.  The  default  is  .FALSE . 

gash2o 

(Logical)  Write  water  vapor  fraction  to  compartment/duct  node  gas  mass 
fraction  output  file.  The  default  is  .FALSE . 

gasco 

(Logical)  Write  carbon  monoxide  mass  fraction  to  compartment/duct 
node  gas  mass  fraction  output  file.  The  default  is  .FALSE . 

gasc 

(Logical)  Write  soot  mass  fraction  to  compartment/duct  node  gas  mass 
fraction  output  file.  The  default  is  .TRUE . 

gasfuel 

(Logical)  Write  unbumed  fuel  mass  fraction  to  compartment/duct  node 
gas  mass  fraction  output  file.  The  default  is  .FALSE . 

It  is  important  to  note  that  there  is  a  limit  to  the  number  of  concurrently  open  files  in 
Fortran.  This  limit  as  implemented  in  FSSIM  is  250  files.  As  will  be  discussed  in  Section  5.1, 
the  output  files  are  limited  to  255  columns  of  data  in  order  to  remain  useable  by  popular  COTS 
spreadsheet  software.  Thus,  given  a  sufficiently  large  problem  one  or  more  files  may  be  needed 
for  each  of  the  above  file  types.  A  fatal  error  will  occur  if  the  user  requests  more  files  to  be 
written  than  can  be  supported. 
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4.2.4  Input  Data  Files  (&fnam) 


The  &fnam  group  is  used  to  specify  names  of  input  files  to  be  read  for  specific 
NAMELIST  parameter  groups  data.  All  parameters  associated  with  this  NAMELIST  group  are 
optional,  as  is  using  additional  files  as  data  input  vehicles.  The  latter  is  controlled  by  other 
NAMELIST  groups.  The  following  is  a  complete  list  of  parameters  set  by  the  &fnam  group: 


OPTIONAL  PARAMETERS 


materialfile 

(Text  string)  Name  of  file  containing  smtrl  NAMELIST  input  data.  The 
default  is  ‘materials .  txt5  and  the  string  must  be  255  characters  or  less. 

fueifile 

(Text  string)  Name  of  file  containing  sfuel  NAMELIST  input  data.  The 
default  is  ‘fuels .  txt’  and  the  string  must  be  255  characters  or  less. 

usefile 

(Text  string)  Name  of  file  containing  suses  NAMELIST  input  data.  The 
default  is  ‘usetypes .  txt5  and  the  string  must  be  255  characters  or  less. 

compositionfile 

(Text  string)  Name  of  file  containing  scmpn  NAMELIST  input  data.  The 
default  is  ‘compositions .  txt5  and  the  string  must  be  255  characters  or 
less. 

fanfile 

(Text  string)  Name  of  file  containing  srfan  NAMELIST  input  data.  The 
default  is  ‘fans .  txt5  and  the  string  must  be  255  characters  or  less. 

detectorfile 

(Text  string)  Name  of  file  containing  &dtct  NAMELIST  input  data.  The 
default  is  ‘detectors .  txt5  and  the  string  must  be  255  characters  or  less. 

suppressorfile 

(Text  string)  Name  of  file  containing  &supr  NAMELIST  input  data.  The 
default  is  ‘suppressors .  txt’  and  the  string  must  be  255  characters  or 
less. 

An  example  of  a  sfnam  NAMELIST  parameter  group  line  is  as  follows: 


&FNAM  materialfile=' newmats .txt' ,  fuelfile=' AllFuels.dat' , 
fanfile=' fans.txt' ,  suppressorfile=' nozzles. dat'  / 


The  example  &fnam  line  changes  the  default  materiaif ile,  fueifile,  and  suppressorfile 
names  to  ‘newmats  .txt5,  ‘AllFuels .  dat5,  and  ‘nozzles .  dat’,  respectively.  The  fanfile  is 
also  set  to  the  default  ‘fans .  txt5  in  the  NAMELIST  group,  though  it  is  unnecessary  to  do  so. 


4.2.5  Tabular  Input  Data  (&curv) 

The  &curv  NAMELIST  group  is  used  by  other  NAMELIST  groups  as  a  means  of 
collecting  tabular  data.  Examples  of  tabular  data  that  may  be  entered  using  the  scurv 
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NAMELIST  group  parameters  include  thermal  conductivity  or  thermal  heat  capacity  that  vary 
with  temperature,  user-defined  fire  with  a  transient  heat  release  rate,  and  fan  performance  curves. 
The  &curv  parameters  are  read  when  directed  to  by  a  referring  NAMELIST  parameter  group. 

The  parameters  associated  with  this  NAMELIST  group  are  required;  however,  the  group 
itself  is  optional.  The  scurv  parameters  are  as  follows: 


REQUIRED  PARAMETER  WHEN  A  &CURV  NAMELIST  GROUP  IS  USED 


id 

(Text  string)  Name  of  &curv  data.  Used  to  associate  the  scurv  data  with 
a  referring  NAMELIST  group.  The  string  must  be  unique  for  each  set  of 
&curv  id  parameters  and  is  used  for  all  lines  used  to  enter  the  tabular 
data  for  a  single  function  or  rule.  There  is  no  default  and  the  maximum 
number  of  characters  in  the  id  string  is  40. 

The  following  parameter  is  optional;  however,  in  nearly  all  instances  the  user  will  supply 
the  appropriate  input  values: 

OPTIONAL  PARAMETER 

xy(2) 

(Real  number  array)  Array  of  two  real  numbers  separated  by  a  comma. 

The  two  numbers  are  the  independent  and  dependent  variables  of  the 
underlying  function  or  rule  used  to  generate  data,  respectively.  Note  that 
the  independent  variable  (x)  must  be  in  ascending  order.  The  default 

Value  is  0.0  for  both  elements  of  the  xy  array. 

A  table  is  created  by  using  multiple  scurv  lines  each  with  the  same  id.  The  number  of 
tabular  entries  is  limited  either  by  the  particular  NAMELIST  parameter  the  data  are  being 
assigned  to  or  the  size  of  an  integer  (215),  whichever  is  smaller.  Note  that  the  number  of  &curv 
lines  is  not  declared  in  the  scntr  NAMELIST  parameter.  Also,  multiple  NAMELIST  parameter 
groups  may  access  the  same  &curv  data  provided  that  they  reference  the  same  &curv  id. 
Typically,  the  scurv  data  used  to  generate  a  particular  table  are  located  beneath  the  NAMELIST 
parameter  group  that  references  the  scurv  data.  Refer  to  subsequent  sections  for  specific  details 
of  NAMELIST  groups  that  may  use  the  &curv  input  format. 

The  &curv  parameter  group  may  be  located  in  the  primary  input  file  or  in  any  of  the 
&fnam  auxiliary  input  files  associated  with  NAMELIST  parameters  that  use  tabulated  data; 
however,  it  must  be  located  in  the  same  file  as  the  referring  NAMELIST  group. 

An  example  of  a  &curv  NAMELIST  parameter  group  line  is  as  follows: 


&CURV  id='A36  steel  thermal  conductivity',  xy=0.0,  40.0  / 
&CURV  id='A36  steel  thermal  conductivity',  xy=100.0,  36.0  / 
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&CURV  id='A36  steel  thermal  conductivity',  xy=500.0,  30,0  / 
&CURV  ±d='A36  steel  thermal  conductivity',  xy=1000,0,  23,0  / 


which  describes  a  table  that  may  be  referenced  by  another  NAMELIST  group  using  the  text 
string ‘A36  steel  thermal  conductivity’  and  containing  four  pairs  of 
independent-dependent  variables. 


4.2.6  Input  Conditional  Control  (sctrl) 

The  sctrl  NAMELIST  group  is  used  to  control  the  status  of  various  types  of  devices, 
namely  doors,  hatches,  fans,  sprinklers,  nozzles,  and  other  active  suppression  components.  The 
control  functions  are  accessed  by  the  specific  NAMELIST  group  that  specifies  the  device  being 
controlled  using  the  sctrl  id.  There  must  be  a  sctrl  parameter  group  associated  with  devices 
set  by  the  srfan  and  ssupr  NAMELIST  groups.  Doors  and  hatches,  set  using  the  sjunc 
NAMELIST  parameter  group,  may  or  may  not  use  a  sctrl  NAMELIST  line.  Subsequent 
sections  describe  these  NAMELIST  parameter  groups  in  detail.  Currently,  the  sctrl  parameter 
group  only  turns  devices  on  and  off  or  open  and  shut.  It  is  expected  that  additional  functions  will 
be  available  in  this  group  that  may  be  used  to  reduce  or  increase  characteristic  functions  of 
various  devices,  such  as  flow  rates  and  concentrations. 

One  parameter  associated  with  the  sctrl  group  is  required  when  the  scurv  is  used: 


REQUIRED  PARAMETER  WHEN  A  &CTRL  NAMELIST  GROUP  IS  USED 


id 


(Long  integer)  Identification  number  for  an  individual  sctrl  group  line. 
This  number  is  referenced  by  other  NAMELIST  parameter  groups  to 
access  the  control  functions  set  by  the  sctrl  line.  The  id  must  be  unique 
among  other  sctrl  id  parameters  used  in  a  particular  FSSIM  simulation 
and  greater  than  zero. 


There  Is  no  requirement  for  the  sctrl  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  sctrl  id 
parameters  is  acceptable. 

The  remaining  parameters  are  optional;  however,  depending  on  the  particular  device 
being  controlled,  they  may  require  setting.  The  following  is  a  complete  list  of  the  optional 
parameters  that  maybe  set  on  the  sctrl  parameter  group  line: 
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OPTIONAL  PARAMETERS 


ctype 

. 

(Integer)  Type  of  control.  The  default  is  1 . 

1  -  On  or  off. 

2  -  Indicates  an  AND  function.  Used  with  inputvar  =  10  or  1 1. 

Other  types  not  yet  implemented. 

trip 

(Logical)  .  true  .  indicates  devices  changes  state  once;  .  false  .  indicates 
device  may  change  states  more  than  once.  The  default  state  is  .  true  . 

initial 

(Logical)  .  true  .  indicates  device  is  initially  on;  .  false .  indicates 
device  is  initially  off.  Used  when  ctype  is  set  to  1.  The  default  state  is 

. FALSE . 

inputvar 

(Integer)  Input  variable  used  to  change  the  state  of  the  device: 

1  -  time  (setpoint  units  are  seconds). 

2  -  compartment  temperature  (setpoint  units  are  K). 

3  -  duct  temperature  (setpoint  units  are  K). 

4  -  compartment  heat  release  rate  (setpoint  units  are  W). 

5  -  compartment  oxygen  mass  fraction  (setpoint  units  are 

kg/kg). 

7  -  Time  delay  from  trip  of  control  function  (s) 

8  -  compartment  surface  temperature  (setpoint  units  are  K) 

10  -  a  detector  (setpoint  is  not  used) 

1 1  -  dual  control  function 

If  the  inputvar  is  10,  the  condition  that  causes  a  state  change  is  detector 
activation.  The  default  is  1  (time). 

setpoint 

(Real  number)  The  condition  (e.g.  a  time,  temperature,  or  heat  release 
rate)  that  causes  the  device  to  change  state.  The  type  of  condition  is  set 
by  inputvar.  Default  is  1.0. 

location ( 2 ) 

(Long  integer)  The  location  (device,  duct,  or  compartment  identification 
number  [id])  of  the  inputvar  that  is  used  to  change  the  state  of  the 
device.  Not  used  if  inputvar  is  set  to  1.  Otherwise: 

•  &COMP  id  parameter  value  if  inputval  is  set  to  2, 4,  or  5  (location 
is  a  compartment).  Only  location  (l)  is  set. 

•  &rdct  id  parameter  value  if  inputvar  is  set  to  3  (location  is  a 
inside  a  duct).  Only  location  (l)  is  set. 

•  sdtct  id  parameter  if  inputvar  is  set  to  10  (location  is  a  specific 
detector).  Only  location  ( l )  is  set. 

•  &ctrl  id  parameter  if  inputvar  is  set  to  7  (location  is  a  control 
function).  Only  location  (l )  is  set. 

•  2  &ctrl  id  parameters  if  inputvar  is  set  to  1 1  (locations  are  two 
control  functions  to  be  ANDed).  Both  location  ( l )  and 
location  { 2 )  are  set. 

The  default  is  1 . 
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outputvar 

(Integer).  Output  parameter  used  when  ctype  is  not  set  to  1.  Not  yet 
implemented.  Default  is  1. 

maxbin 

(Integer).  Function  value  when  device  is  on.  For  use  when  ctype  is  not 
set  to  1.  Not  yet  implemented.  Default  is  1. 

minbin 

(Integer).  Function  value  when  device  is  off.  For  use  when  ctype  is  not 
set  to  1 .  Not  yet  implemented.  Default  is  0, 

As  noted,  sctrl  only  turns  devices  on  or  off  in  this  version  of  FSSIM;  however,  it  is 
expected  that  it  may  be  used  to  perform  other  types  of  function  modifications.  As  such, 
parameters  that  are  not  yet  implement  but  are  expected  to  be  used  for  this  purpose  are  listed 
above  for  completeness. 

Typically,  a  sprinkler  or  suppression  device  will  be  actuated  using  conditions  or  a 
detector  state  within  the  space  where  the  suppression  device  is  located.  Fans  are  typically  turned 
on  or  off  based  on  a  detection  state  in  a  compartment  or  duct  or  a  duct  temperature  condition. 
Note  that  a  device  may  initially  be  off  then  turned  on  and  vice  versa.  Also  note  that  when  the 
inputvar  is  set  to  1  (time),  the  trip  parameter  is  inherently  .  true  . ,  though  it  may  actually  be 
Set  to  .FALSE. 

Multiple  NAMELIST  parameter  groups  may  access  the  same  sctrl  parameter  group 
data  provided  they  reference  the  same  sctrl  id.  Typically,  each  sctrl  parameter  group  is 
located  beneath  the  NAMELIST  parameter  group  that  references  it  in  the  input  file,  though  this 
is  not  a  requirement. 

An  example  of  a  sctrl  NAMELIST  parameter  group  line  is  as  follows: 


SCTRL  id=3,  ctype=l,  inputvar=l,  initial=. FALSE.,  setpoint=120 . 0  / 


The  above  sctrl  line  has  an  id  reference  of  3,  the  control  is  on/off  with  the  device 
initially  being  off,  and  the  device  is  turned  on  at  120  seconds.  Other  NAMELIST  groups  thus 
access  this  control  using  the  integer  3.  The  sctrl  example  is  typical  of  an  exhaust  or  supply  fan. 

4.3  Compartments 

Compartments  are  defined  using  parameter  groups  that  set  the  dimensions  and  volume  of 
a  space,  the  exposed  materials  and  surface  area,  and  the  dimensions  of  openings  between  two 
spaces  and  between  a  space  and  the  surroundings  that  lie  beyond  the  modeled  domain  (ambient). 


4.3.1  Compartment  Definition  (scomp) 

The  scomp  group  allows  the  user  to  set  the  physical  dimensions  of  a  compartment.  One 
scomp  group  is  required  for  each  compartment  in  the  simulated  configuration.  There  must  be  at 
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least  one  compartment  for  FSSIM  to  successfully  run.  There  is  one  parameter  in  each  &comp 
group  line  that  must  be  set: 


REQUIRED  PARAMETER 

(Long  integer)  Identification  number  for  an  individual  scomp  parameter 
group  line.  The  id  must  be  unique  among  other  &comp  parameter  id 
values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 
is  no  default  value. 


There  is  no  requirement  for  the  &comp  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  scomp  id 
parameters  is  acceptable. 

The  remaining  parameters  are  optional,  though  in  most  instances  compartment 
dimensions  other  than  the  default  values  will  be  set.  The  following  is  a  complete  list  of  the 
optional  parameters  that  may  be  set  on  the  &comp  parameter  group  line: 


OPTIONAL  PARAMETERS 


compname 

(Text  string)  Name  of  compartment.  The  default  is  ‘compartment  id’ 
where  id  is  the  compartment  identification  number.  The  compartment 
name  is  used  primarily  for  labeling  output  data.  The  maximum  number  of 
characters  for  the  compartment  name  is  40. 

dx 

(Real  number)  Nominal  compartment  dimension  along  the  keel  or  in  the 
direction  parallel  to  the  x-axis  (m).  The  default  value  is  1 .0  m. 

dy 

(Real  number)  Nominal  compartment  dimension  along  the  beam  or  in  the 
direction  oriented  parallel  to  the  y-axis  (m).  The  default  value  is  1 .0  m. 

dz 

(Real  number)  Nominal  compartment  height  or  dimension  parallel  to  the 
z-axis  (m).  The  default  value  is  1.0  m. 

z 

(Real  number)  Nominal  elevation  of  floor  above  reference  datum  (m). 

The  default  is  0.0  m. 

vol 

(Real  number)  Free  air  volume  in  compartment  (nr3).  The  default  is  1.0 
m3. 

pres 

(Real  number)  Initial  pressure  in  compartment  (Pa).  The  default  is  pamb 
set  in  the  &exec  NAMELIST  parameter  group. 

rho 

(Real  number)  Initial  density  in  compartment  (K).  The  default  is 
calculated  from  tamb  set  in  the  &exec  NAMELIST  parameter  group. 

Note  do  not  specify  both  temp  and  rho.  If  both  are  set,  the  rho  will 
override  the  temp. 

25 


temp 

(Real  number)  Initial  temperature  in  compartment  (K).  The  default  is 
tamb  set  in  the  &exec  NAMELIST  parameter  group.  Note  do  not  specify 
both  temp  and  rho.  If  both  are  set,  the  rho  will  override  the  temp. 

usetype 

(Long  integer)  suses  NAMELIST  group  identification  number  (id)  of  a 
&uses  line  that  describes  the  combustible  fuel  load  and  the  type  of  fuel 
( &  fuel  id)  in  the  compartment.  A  value  of -1  indicates  that  there  is  no 
combustible  fuel  load  or  that  there  is  an  initial  fire  in  the  compartment. 

The  default  value  is  -1 . 

zone 

(Long  integer)  The  zone  number  where  the  compartment  is  located.  The 
default  is  calculated  by  FSSIM  and  the  value  entered  here  is  overridden 
regardless.  This  parameter  should  thus  not  be  used. 

The  free  air  volume  may  be  less  than  or  greater  than  the  compartment  volume  calculated 
from  the  nominal  dimensions  dx,  dy,  and  dz.  The  reference  datum  for  the  parameter  z  is  the 
elevation  that  the  ambient  pressure  (pamb)  and  temperature  (tamb)  set  in  the  &exec  group 
corresponds  to. 

Multiple  compartments  may  use  the  same  suses  line  (the  usetype  parameters  need  not 
be  unique  among  &comp  lines). 

An  example  of  a  scomp  NAMELIST  parameter  group  line  is  as  follows: 


&COMP  id=l,  z=l . 05,  dx=7 . 62,  dy=3.1,  dz=3.76,  vol=88.82,  temp=311.2, 
compname^' sample space  lr  / 


The  above  scomp  line  establishes  a  compartment  that  is  nominally  7.62  m  by  3.1  m  by 
3.76  m  high  and  has  a  free  air  volume  of  88.82  m3.  The  floor  of  this  compartment  is  located  1 .05 
m  above  the  location  where  the  ambient  pressure  and  ambient  temperature  are  defined  on  the 
&exec  NAMELIST  line.  The  initial  temperature  of  this  space  is  31 1.2  K,  possibly  different  from 
the  ambient  temperature,  and  the  compartment  name  is  ‘sampiespace  i\ 

4.3.2  Compartment  Opening  Definition  (& junc) 

The  &  junc  group  is  used  to  set  vent  opening  flow  paths  between  compartments  and  other 
compartments  and  between  compartments  and  the  surrounding  ambient  space.  Openings  set 
using  the  &  junc  NAMELIST  group  are  optional.  They  may  be  horizontal  as  in  a  deck/overhead 
boundary  or  vertical  as  in  a  bulkhead/wall  boundary.  Horizontal  openings  allow  vertical  flow  and 
vertical  vents  allow  horizontal  flow.  Flows  through  openings  set  using  the  &  junc  NAMELIST 
group  are  not  directly  forced.  Flows  that  are  forced  via  a  mechanical  system  or  systems  are  set 
using  the  &rdci,  srnod,  and  &rfan  NAMELIST  groups. 

If  an  opening  is  declared  using  a  &junc  parameter  group,  then  three  parameters  are 
required  to  be  set: 
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REQUIRED  PARAMETERS  WHEN  A  &  JUNC  NAMELIST  GROUP  IS  USED 


id 

(Long  integer)  Identification  number  for  an  individual  &  junc  parameter 
group  line.  The  id  must  be  unique  among  other  &  junc  parameter  id 
values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 
is  no  default  value. 

location (2) 

(Long  integer  array)  An  array  of  two  long  integers  separated  by  a 
comma.  The  array  references  the  compartment  identification  numbers  (id 
in  a  &comp  NAMELIST  group)  of  the  compartments  connected  by  the 
opening.  A  value  of-1  for  one  of  the  values  indicates  the  surrounding 
ambient  space.  At  least  one  compartment  id  parameter  used  on  an  &comp 
entry  must  be  specified,  and  a  compartment  must  not  be  connected  to 
itself  (i.e.,  parameter  array  elements  must  not  be  the  same).  There  are  no 
default  values. 

height 

(Real  number)  Opening  dimension  corresponding  to  the  height  for 
vertical  openings  (m).  The  default  value  is  set  0.0  m  to  avoid  causing 
memory  allocation  errors;  however,  for  a  vertical  opening  this  parameter 
must  be  set  to  a  positive  non-zero  value  to  correctly  simulate  flow 
through  the  opening.  For  a  horizontal  opening  this  value  should  remain 

0.0  m. 

There  is  no  requirement  for  the  &  junc  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  &  junc  id 
parameters  is  acceptable.  There  is  also  no  required  order  for  the  location  parameter  array 
elements.  If  an  opening  is  connected  to  ambient  (i i.e .,  the  location  parameter  array  element  is 
set  to  —1),  it  may  be  the  first  or  second  element  in  the  array. 

The  remaining  &  junc  parameters  are  optional: 


OPTIONAL  PARAMETERS 


juncname 

(Text  string)  Name  of  junction.  The  default  is  ‘junction  id’  where  id  is 
the  junction  identification  number.  The  junction  name  is  used  primarily 
for  labeling  output  data.  The  maximum  number  of  characters  for  the 
junction  name  is  40. 

elevation 

(Real  number)  Vertical  elevation  of  the  opening  using  the  same 
coordinate  system  as  for  the  compartment  elevations.  This  elevation 
corresponds  to  the  bottom  of  the  vent  for  a  vertical  opening  and  the 
height  of  the  plane  at  which  a  vent  is  located  for  a  horizontal  opening. 

The  default  value  is  set  0.0  m  to  avoid  causing  memory  allocation  errors; 
however,  this  parameter  must  be  set  to  correctly  simulate  flow  through 
the  opening. 

area 

(Real  number)  The  flow  area  of  the  vent  (m2).  The  default  value  is  1 .0 
m2. 

27 


door 

(Long  integer)  Control  function  &ctrl  id  used  to  open  or  close  a  vent.  A 
value  of — 1  indicates  that  the  vent  is  always  open.  The  default  value  is  -1. 

orientation 

(Short  integer)  The  orientation  of  the  opening  with  respect  to  the  first 
element  of  the  location  array; 

1  -  overhead  or  in  the  plane  with  maximum  compartment  z- 

dimension. 

2  -  deck/floor  or  in  the  plane  with  minimum  compartment  z- 

dimension. 

3  -  fore  or  in  the  plane  with  maximum  compartment  x-dimension. 

4  -  aft  or  in  the  plane  with  minimum  compartment  x-dimension. 

5  -  port/larboard  or  in  the  plane  with  maximum  compartment  y- 
dimension. 

6  -  starboard  or  in  the  plane  with  minimum  compartment  y- 

dimension. 

The  default  value  is  3  (fore  opening). 

bidirectional 

(Logical)  .  true  .  indicates  that  the  opening  has  bidirectional  flow  and 
.  false  .  indicates  that  the  flow  direction  in  an  opening  is  one-way.  The 
default  state  is  .  true  . 

vflow 

(Logical)  .  true  .  indicates  that  opening  allows  vertical  flow,  and 
.  false  .  indicates  that  the  opening  does  not  allow  vertical  flow.  The 
default  state  is  .  true  .  if  the  orientation  is  1  or  2. 

kloss 

(Real  number)  The  flow  loss  coefficient  for  discharge  through  an 
opening.  This  term  may  be  calculated  using  the  following  equation 
[Sabersky  et  ah,  1971;  Brater  et  ah,  1996]: 

Ld 

where  K  is  the  flow  loss  coefficient  and  cj  is  the  discharge  coefficient. 

The  discharge  coefficient  is  a  function  of  the  opening  geometry  and  the 
Reynolds  Number.  The  discharge  coefficient  is  tabulated  for  various 
types  of  openings  [e.g. ,  Idelchik,  1986;  Sabersky  et  ah,  1971;  Brater  et 
ah,  1996].  The  discharge  coefficient  for  typical  openings  in 
compartments  typically  between  0.6-0.7  [Klote,  1995].  The  default  value 
for  kloss  is  2.04,  which  corresponds  to  a  discharge  coefficient  of  about 
0.7. 

In  the  case  of  an  opening  that  spans  a  compartment  wall  (such  as  might 
occur  when  breaking  a  hallway  into  multiple  compartments),  kloss 
should  be  set  to  zero. 
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length 

(Real  number)  The  inertial  length  associated  with  an  opening  (m).  The 
inertial  length  corresponds  to  an  acceleration  distance  of  a  fluid  mass 
located  in  two  volumes  (presumably  with  different  pressure  profiles) 
connected  by  an  opening.  The  recommended  (and  default)  inertial  length 
is  the  horizontal  or  vertical  distance  between  the  center  of  two  connected 
volumes,  or  if  one  volume  is  ambient,  twice  the  horizontal  or  vertical 
distance  between  the  center  of  the  compartment  to  the  opening. 

zone 

(Long  integer)  The  zone  number  where  the  compartment  is  located.  The 
default  is  calculated  by  FSSIM  and  the  value  entered  here  is  overridden 
regardless.  This  parameter  should  thus  not  be  used. 

width 

(Real  number)  The  width  of  the  opening  (m).  The  default  is  4a  (m), 
where  A  is  the  flow  area  of  the  vent,  area.  This  parameter  is  currently  not 
used  by  FSSIM. 

If  an  opening  has  vertical  flow,  then  the  orientation  of  the  opening  must  either  be  in  the 
deck/floor  or  the  overhead  and  vf  low  must  be  set  to  .  true  .  Otherwise,  vertical  flow  will  not  be 
permitted  through  the  opening. 

Most  compartment  vents  should  be  modeled  as  bidirectional.  Some  cases  may  arise 
where  the  flow  direction  is  expected  to  be  in  one  direction  only,  such  as  spaces  that  are 
pressurized.  When  this  condition  arises,  setting  the  bidirectional  parameter  to  .  false  .  may 
reduce  the  number  of  calculations  performed  by  FSSIM  and  possibly  decrease  the  solution  time. 
This  parameter  should  not  be  set  to  .  false  .  if  it  is  not  known  whether  the  flow  is  bidirectional 
or  if  it  is  possible  that  bidirectional  flow  will  be  result. 

The  inertial  length  is  used  in  the  FSSIM  calculation  to  influence  the  acceleration  of  a 
fluid  as  it  moves  toward  and  through  an  opening.  The  default  inertial  length  is  based  on  other 
work  with  such  flow  calculations  [Guantt  et  al.,  2000]  and  should  be  adequate  for  most 
configurations.  Situations  that  may  require  modification  of  this  parameter  could  include  spaces 
connected  by  an  elongated  opening  such  that  the  acceleration  distance  is  greater  than  what  is 
determined  by  default  or  cases  where  the  acceleration  effects  on  the  solution  are  not  desired,  in 
which  case  the  inertial  length  should  be  set  to  a  small  positive  value. 

FSSIM  will  simulate  a  wide  variety  of  vent  conditions,  including  improperly  specified 
vent  height  and  elevations  (i.e.,  height  plus  elevation  exceed  vertical  dimension  of  compartment). 
As  such,  considerable  care  is  necessary  to  ensure  the  vent  specified  accurately  represents  the 
configuration  desired. 

An  example  of  a  &  junc  NAMELIST  parameter  group  line  is  as  follows: 

&JUNC  id=7/  location^lO, -1,  height=2. 13,  area=1.94,  elevation==0 . 0, 
orientation=3,  vflow= . FALSE . ,  juncname=' Ordinary  door'  / 


The  above  &  junc  line  sets  an  opening  between  a  compartment  with  a  &comp  id  parameter 
equal  to  10  and  the  surrounding  ambient  (-1  ).  The  opening  identification  id  is  7  and  is  referred 
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to  in  output  as  ‘ordinary  door’.  The  opening  is  located  in  the  forward  compartment  boundary 
and  is  thus  vertical.  It  is  2.13  m  tall,  the  bottom  is  located  at  elevation  of  0  m,  and  it  has  an  area 
of  1.94  m2.  The  flow  is  bidirectional  (default  condition)  and  not  vertical  (by  default  and  necessity 
[orientation  set  to  3,  a  vertical  boundary],  but  also  specified  with  vf  low  parameter  on  example 
sjunc  line).  The  vent  is  always  open  (the  default  door  parameter  state  of .  false,  is  used).  The 
default  flow  loss  coefficient  of  2.04  is  used  and  the  inertial  length  is  calculated  by  FSSIM. 


4.3.3  Compartment  Surface  Definition  (ssurf) 

The  &surf  group  is  used  to  define  the  boundary  surfaces  in  a  space  for  calculating 
thermal  flows.  A  minimum  of  six  surfaces  must  be  set  for  each  compartment  set  using  the  scomp 
parameter  group,  one  for  each  possible  orientation  (overhead/ceiling,  deck/floor,  fore,  aft, 
port/larboard,  and  starboard).  FSSIM  will  run  if  less  than  the  minimum  number  of  surfaces  is 
specified;  however,  the  results  will  not  be  correct.  The  user  should  ensure  proper  surface 
definition. 

Three  parameters  are  required  to  be  set  on  each  ssurf  line: 


REQUIRED  PARAMETERS 


id 

(Long  integer)  Identification  number  for  an  individual  ssurf  parameter 
group  line.  The  id  must  be  unique  among  other  ssurf  parameter  id 
values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 
is  no  default  value. 

stype 

(Short  integer)  The  type  of  surface: 

1  -  The  surface  is  between  two  compartments  or  a  compartment 

and  the  surrounding  ambient. 

2  -  The  surface  is  between  a  compartment  and  a  duct  space. 

There  is  no  default  value.  Note  that  compartment-duct  surfaces  are  not 
yet  implemented. 
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location 

(Long  integer  array)  An  array  of  two  long  integers  separated  by  a 
comma.  The  array  references  the  compartment  identification  numbers  (id 
in  a  &comp  NAMELIST  group)  if  stype  is  set  to  1  or  a  compartment 
identification  number,  and  a  duct  identification  number  (id  in  a  &rdct 
NAMELIST  group)  if  stype  is  set  to  2.  The  duct  id  value  must  be  the 
second  element  in  the  parameter  array  when  stype  is  set  to  2.  A 
compartment  id  value  of-1  indicates  the  surrounding  ambient  space.  If 
stype  is  set  to  1  at  least  one  compartment  id  value  must  be  positive  and 
the  parameter  array  elements  must  not  be  equal  (i.e.,  a  surface  within  a 
compartment).  There  are  no  default  values.  A  surface  shared  between  two 
compartments  should  only  be  specified  once.  The  order  of  the 
compartment  id  parameters  is  not  significant;  however  the  orientation 
parameter  must  be  properly  set.  Note  that  compartment-duct  surfaces  are 
not  yet  implemented. 

orientation 

(Short  integer)  The  orientation  of  the  surface  with  respect  to  the 
compartment  identified  by  the  first  element  of  the  location  array: 

1  -  overhead/ceiling  or  in  the  plane  with  maximum  compartment 

z-dimension. 

2  -  deck/floor  or  in  the  plane  with  minimum  compartment  z- 

dimension. 

3  -  fore  or  in  the  plane  with  maximum  compartment  x-dimension. 

4  -  aft  or  in  the  plane  with  minimum  compartment  x-dimension. 

5  -  port  or  in  the  plane  with  maximum  compartment  y-dimension. 

6  -  starboard  or  in  the  plane  with  minimum  compartment  y- 

dimension. 

There  is  no  default  value  for  this  parameter. 

There  is  no  requirement  for  the  ssurf  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  &surf  id 
parameters  is  acceptable.  There  is  also  no  required  order  for  the  location  parameter  array 
elements  when  stype  is  set  to  1.  If  the  surface  is  adjacent  to  ambient  (i.e.,  one  of  the  location 
parameter  array  elements  is  set  to  -1),  it  may  be  the  first  or  second  element  in  the  array.  If  stype 
is  set  to  2,  then  the  duct  identification  (srcdt  id)  must  be  the  second  element  in  the  parameter 
array.  The  first  element  may  be  either  a  compartment  identification  number  (scomp  id)  or  —1  (the 
ambient  surroundings).  Surfaces  must  separate  different  compartment  spaces  if  stype  is  set  to  1, 
such  that  a  surface  is  not  located  in  a  single  compartment  or  space. 

•  A  compartment  must  have  at  least  six  surfaces  defined,  one  for  each  orientation,  in  order 
for  the  radiation  heat  transfer  calculation  to  function  properly.  The  surfaces  may  be  shared  with 
another  compartment  or  they  may  be  between  the  compartment  and  ambient.  Note  that  if  the 
surface  is  shared,  and  the  compartment  id  is  the  second  element  in  the  location  array,  then  the 
surface  orientation  for  the  compartment  is  the  opposite  of  what  is  specified  by  the  orientation 
parameter. 

The  remaining  &surf  parameters  are  optional: 
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OPTIONAL  PARAMETERS 


surfnante 

(Text  string)  Name  of  surface.  The  default  is  ‘surface  c  location 
(1)  c  location  (2)  o  orientation’ where  location (1)  is  the  first 
compartment  identification  number  specified  in  the  location  parameter 
array,  location  (2)  is  the  second  compartment  identification  number 
or  the  first  duct  identification  number,  and  orientation  is  the 
orientation  parameter  value.  The  surface  name  is  used  primarily  for 
labeling  output  data.  The  maximum  number  of  characters  for  the 
surface  name  is  40. 

area 

(Real  number)  The  area  of  the  surface  (m2).  The  default  value  is  1 .0 
m2. 

composition 

(Long  integer)  scmpn  NAMELIST  group  identification  number  (id)  of 
a  &cmpn  line  that  describes  the  surface  composition  (materials).  The 
default  value  is  1. 

fake 

(Logical)  .  false,  indicates  that  the  surface  is  a  solid  boundary  and 
.  true  .  indicates  an  open  boundary  that  allows  thermal  radiation  heat 
flow  but  does  not  perform  the  one-dimensional  heat  transfer 
calculation.  An  example  of  the  latter  type  of  surface  would  be  the  open 
end  of  corridor  or  a  large  vent  opening  in  a  compartment.  The  default 
state  is  .false. 

reverse_composition 

(Logical)  .  true  .  indicates  that  the  composition  definition  for  the 
surface  is  in  the  opposite  direction  of  the  compartment  definition,  i.e. 
the  composition  layer  as  defined  using  the  first  element  of  both  the 
&cmpn  dx  and  materials  parameter  arrays  are  associated  with  the 
compartment  identified  by  second  element  of  the  location  array.  The 
default  state  is  .  false  . 

Fake  surfaces  are  used  to  calculate  the  radiant  heat  flow  through  openings.  This  may  be  a 
significant  energy  quantity  in  spaces  where  the  temperature  is  high  and  the  opening  area  is 
proportionately  large.  Note  that  there  is  no  practical  limit  to  the  number  of  surfaces  that  may  be 
set  in  a  compartment.  The  total  surface  area  should  nominally  be  about  equal  to  the  boundary 
surface  area. 

An  example  of  a  ssurf  NAMELIST  parameter  group  line  is  as  follows: 

&SURF  id=3,  stype=l,  loeation=3, 2,  area=1.94,  orientation=3,  compositional, 
surfname=# Wall  section*  / 


The  above  &surf  has  an  identification  (id)  value  of  3  and  is  also  identified  by  the  text 
string  ‘wall  section’.  The  surface  is  between  two  compartments  with  id  parameters  set  to  3 
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and  2.  The  surface  is  located  in  the  plane  with  the  maximum  x-dimension  (fore)  of  the 
compartment  with  an  id  parameter  set  to  3  and  has  a  surface  area  of  1 .94  m2.  The  heat  transfer 
calculations  will  be  performed  for  this  surface  (the  default  state  for  the  fake  parameter 
[ .  false  .  ]  is  used)  and  the  thickness  and  specific  materials  that  make  up  the  surface  boundary 
are  described  by  a  scmpn  line  with  an  id  parameter  of  1 .  If  the  fake  parameter  is  set  to  .  true  . , 
then  the  composition  parameter  is  not  used. 


4.3.4  Compartment  Surface  Composition  (scmpn) 

The  &cmpn  NAMELIST  group  is  used  to  set  the  type,  number,  and  thickness  of  materials 
that  make  up  a  surface  boundary.  The  &cmpn  NAMELIST  group  is  optional;  however,  if  there  is 
a  surface  that  is  not  fake,  then  at  least  one  scmpn  line  must  be  provided  with  an  identification 
number  that  matches  the  ssurf  composition  parameter  value. 

A  scmpn  NAMELIST  line  has  one  parameter  that  must  be  set: 


REQUIRED  PARAMETER  WHEN  A  SCMPN  NAMELIST  GROUP  IS  USED 


id 

(Long  integer)  Identification  number  for  an  individual  scmpn  parameter 
group  line.  The  id  must  be  unique  among  other  scmpn  parameter  id 
values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 
is  no  default  value. 

There  is  no  requirement  for  the  scmpn  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  scmpn  id 
parameters  is  acceptable. 

The  remaining  scmpn  parameters  are  optional: 

OPTIONAL  PARAMETERS 

cmpnname 

(Text  string)  Name  of  composition.  The  composition  name  is  used 
primarily  to  identify  the  appropriate  composition  input  data  to  read  from 
a  file  specified  by  sfnam  compositions  file  parameter.  The 
composition  name  is  also  used  for  labeling  output  data.  The  maximum 
number  of  characters  for  the  composition  name  is  40.  The  default  is 
‘composition  id’,  where  id  is  the  composition  identification  number. 
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fileread 

(Logical)  .  true  .  indicates  that  the  composition  data  should  be  read  from 
the  file  set  by  the  sfnam  compositionsfile  parameter  and  .  false  . 
indicates  that  the  composition  data  are  contained  within  the  scomp  id  line 
located  in  the  main  input  file.  The  default  state  is  .  true  .  If  the  cmpnname 
parameter  is  not  set  (t.e.,  the  default  is  used),  then  the  composition  data 
will  not  be  read  from  the  file  set  by  the  &fnam  compositionsfile 
parameter. 

layers 

(Short  integer)  The  number  of  material  layers  that  make  up  the 
composition  (and  referencing  surface  boundary).  The  default  value  is  1 
and  the  maximum  number  of  layers  permitted  is  50. 

materials 

(Long  integer  array)  An  array  of  long  integers  separated  by  commas 
consisting  of  the  material  identification  numbers  (&mtrl  id  parameter 
values)  that  make  up  each  layer  of  the  composition.  The  number  of 
elements  must  be  equal  to  the  layers  parameter  value  and  thus  not 
exceed  50.  The  first  material  identified  in  the  array  will  have  a  surface 
located  in  the  first  compartment  set  using  the  first  element  of  the  &surf 
location  parameter.  The  last  material  will  have  an  exposed  surface 
located  in  the  second  compartment  or  in  the  duct  set  using  the  second 
element  of  the  ssurf  location  parameter  array.  The  intervening 
materials  will  be  sequentially  adjacent;  that  is,  materials  array  element 

2  will  be  physically  located  between  materials  array  elements  1  and  3. 
The  default  value  for  this  parameter  is  1. 

dx 

(Real  number  array)  An  array  of  real  numbers  separated  by  commas 
consisting  of  the  thickness  of  each  layer  (m).  The  number  of  elements 
must  be  equal  to  the  layers  parameter  value  and  thus  not  exceed  50.  The 
first  value  corresponds  to  the  first  element  in  the  materais  array,  the 
second  value  corresponds  the  second  element,  and  so  forth.  The  default 
value  is  0.1  m. 

Input  associated  with  the  &cmpn  NAMELIST  group  maybe  located  in  a  separate  input 
file  with  a  name  specified  using  the  sfnam  compositionsfile  parameter.  If  a  different  input  file 
is  to  be  used,  then  the  following  must  be  set  on  the  scmpn  line  in  the  primary  input  file: 

•  f  ileread  must  be  set  to  .  true  .  (the  default  condition). 

•  cmpnname  must  be  set  to  a  string  that  can  be  found  in  the  file  set  by  the  &fnam 
compositionsfile  parameter. 

•  id  must  be  positive  and  unique  among  the  scmpn  id  parameters. 

The  input  file  specified  by  &fnam  compositionsfile  parameter  should  contain  a  scmpn  line 
with  the  cmpnname  equal  to  that  specified  in  the  primary  input  file.  The  line  should  contain  the 
additional  layers,  materials,  and  thickness  parameters  for  the  composition  material.  The  id 
parameter  must  not  be  reset  in  the  input  file  specified  by  sfnam  compositionsfile  parameter. 
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Two  examples  of  the  scmpn  NAMELIST  parameter  are  provided,  one  in  which  all 
information  is  entered  in  the  primary  input  file  and  one  in  which  a  secondary  input  file  is  used. 
The  first  example  is  as  follows  and  is  contained  entirely  within  the  primary  input  file: 

&CMPN  id=8,  fileread=. FALSE. ,  cmpnname=' Ins-steel-ins-paint' ,  layers=4, 
materials=l, 2,1,3,  dx=0.05,  0.01,  0.05,  0.001  / 


The  above  &cmpn  line  is  not  read  from  a  file.  It  has  an  identification  (id)  value  of  8  and  is 
also  identified  by  the  text  string  ‘ins-steel-ins-paint’.  The  composition  has  four  layers.  The 
first  and  third  layer  are  0.05  m  thick  and  consist  of  the  material  with  a  smtrl  id  of  1.  The  second 
layer  is  0.01  m  thick  and  uses  the  material  with  a  smtrl  id  of  2,  and  the  fourth  layer  is  0.001  m 
thick  and  uses  the  material  with  a  smtrl  id  of  3. 

The  second  example  uses  the  input  file  set  by  the  sfnam  compositionsf  ile  parameter  to 
enter  the  composition  layer  data.  The  following  line  appears  in  the  primary  input  file: 

&CMPN  id=8,  fileread=.TRUE. ,  cmpnname=f Ins-steel-ins-paint' / 


The  following  line  would  then  appear  in  the  file  set  by  the  &fnam  compositionsf  ile  parameter: 

&CMPN  cmpnname=' Ins-steel-ins-paint' ,  layers=4,  materials=l, 2, 1, 3,  dx=0.05, 
0.01,  0.05,  0.001  / 


The  second  scmpn  NAMELIST  example  is  identical  to  the  first  example  in  terms  of  the  layer 
information  entered. 


4.3.5  Surface  Composition  Material  Definition  (smtrl) 

The  &mtrl  NAMELIST  group  is  used  to  enter  the  thermal  properties  of  the  surface 
boundary  materials.  The  smtrl  NAMELIST  group  is  optional;  however,  if  there  are  one  or  more 
surfaces  that  are  not  fake,  then  there  must  be  at  least  one  scmpn  line  and  therefore  at  least  one 
&mtrl  line.  The  &mtrl  line  or  lines  must  have  an  identification  number(s)  that  corresponds  to  the 
&cmpn  materials  parameter  value(s). 

A  &mtrl  NAMELIST  line  has  one  parameter  that  must  be  set: 
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REQUIRED  PARAMETER  WHEN  A  &MTRL  NAMELIST  GROUP  IS  USED 

(Long  integer)  Identification  number  for  an  individual  &mtrl  parameter 
group  line.  The  id  must  be  unique  among  other  &mtrl  parameter  id 
values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 
is  no  default  value. 


There  is  no  requirement  for  the  smtrl  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  smtrl  id 
parameters  is  acceptable. 

The  remaining  smtrl  parameters  are  optional: 


OPTIONAL  PARAMETERS 


matname 

(Text  string)  Name  of  material.  The  material  name  is  used  primarily  to 
identify  the  appropriate  material  input  data  to  read  from  a  file  specified 
by  &fnam  materiaif  ile  parameter.  The  material  name  is  also  used  for 
labeling  various  output  data.  The  maximum  number  of  characters  for  the 
material  name  is  40.  The  default  is  ‘material  id’,  where  id  is  the 
material  identification  number. 

fileread 

(Logical)  .  true  .  indicates  that  the  material  data  should  be  read  from  the 
file  set  by  the  &  fnam  materiaifile  parameter  and  .false,  indicates 
that  the  material  data  are  contained  within  the  &mtrl  id  line  located  in 
the  main  input  file.  The  default  state  is  .  true  .  If  the  matname  parameter 
is  not  set  (i.e.,  the  default  is  used),  then  the  material  data  will  not  be  read 
from  the  file  set  by  the  &fnam  materiaifile  parameter. 

k 

(Real  number)  The  thermal  conductivity  of  the  material  (W/m-K).  The 
default  value  is  40.0  W/m-K. 

cp 

(Real  number)  The  heat  capacity  of  the  material  (J/kg-K).  The  default 
value  is  486.0  J/kg-K. 

rho 

(Real  number)  The  density  of  the  material  (kg/m3).  The  default  value  is 
7,753  kg/m3. 

e 

(Real  number)  The  surface  emissivity  of  the  material.  The  default  value 
is  0.9. 
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ktabie 

(Logical)  .  true  .  indicates  that  the  thermal  conductivity  data  are  input  as 
a  tabular  function  of  temperature  using  one  or  more  (up  to  500)  scurv 
lines.  .  false  .  indicates  that  the  thermal  conductivity  is  constant  and  is 
equal  to  the  value  set  by  the  k  parameter. 

The  thermal  conductivity  is  also  assumed  constant  for  temperatures  that 
are  above  or  below  the  temperature  range  set  by  the  scurv  lines.  If  the 
temperature  is  less  than  the  minimum  temperature  value,  then  the  thermal 
conductivity  is  assumed  equal  to  the  thermal  conductivity  specified  for 
the  minimum  temperature.  The  reverse  is  true  if  the  temperature  is 
greater  than  the  maximum  temperature  value. 

The  default  state  is  .  false  . 

cptable 

(Logical)  .  true  .  indicates  that  the  heat  capacity  data  are  input  as  tabular 
function  of  temperature  using  one  or  more  (up  to  500)  scurv  lines. 

.  false  .  indicates  that  the  heat  capacity  is  constant  and  is  equal  to  the 
value  set  by  the  cp  parameter. 

The  heat  capacity  is  also  assumed  constant  for  temperatures  that  are 
above  or  below  the  temperature  range  set  by  the  scurv  lines.  If  the 
temperature  is  less  than  the  minimum  temperature  value,  then  the  heat 
capacity  is  assumed  equal  to  the  heat  capacity  specified  for  the  minimum 
temperature.  The  reverse  is  true  if  the  temperature  is  greater  than  the 
maximum  temperature  value. 

The  default  state  is  .  false  . 

kname 

(Text  string)  The  name  of  the  scurv  line(s)  that  contain  the  thermal 
conductivity  data.  There  is  no  default  and  the  maximum  number  of 
characters  is  40.  This  parameter  is  not  used  if  ktabie  is  set  to  .  false. 

cpname 

(Text  string)  The  name  of  the  scurv  line(s)  that  contain  the  heat  capacity 
data.  There  is  no  default  and  the  maximum  number  of  characters  is  40. 

This  parameter  is  not  used  if  cptable  is  set  to  .  false  . 

The  smtrl  parameters  may  be  read  from  a  separate  input  file  as  specified  by  the  sfnam 
materiaif  ile  parameter.  The  method  is  identical  to  that  described  for  the  scmpn  NAMELIST 
parameter  group  in  Section  4.3.4  and  is  not  repeated  here. 

The  thermal  conductivity  and  heat  capacity  of  a  material  entered  using  the  smtrl 
NAMELIST  group  may  be  constant  or  a  function  of  temperature.  If  constant,  then  ktabie  or 
cptable  must  be  set  to  .  false  .  If  either  is  a  function  of  temperature,  then  ktabie  and/or 
cptabie  must  be  set  to  .  true  .  and  one  or  more  scurv  lines  with  the  id  parameter  set  to  kname 
or  cpname  must  be  provided.  The  scurv  lines  must  also  be  located  in  the  same  file  in  which  the 
cptable  or  ktabie  parameters  are  set.  Refer  to  Section  4.2.4  for  a  description  of  the  scurv 
NAMELIST  group. 
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The  following  is  an  example  of  a  &mtrl  NAMELIST  group  and  the  associated  scurv 
lines  that  use  the  input  file  set  by  the  &fnam  materiaifile  parameter  to  enter  the  composition 
layer  data: 


&MTRL  id=8,  fileread=.TRUE. ,  matname=' Insulation' /  This  is  located  in  the 
primary  input  file 


The  following  lines  would  then  appear  in  the  file  set  by  the  &fnam  materiaifile  parameter: 


&MTRL  matname=' Insulation' ,  e=0.91,  rho=128.0,  cp=1157.0,  cptable=. FALSE. , 
ktable=.TRUE. ,  kname=' Ins_k'  / 

&CURV  id='Ins_k',  xy=  0.0,  0.03  / 

&CURV  id='Ins_k',  xy=  1000.0,  0.25  / 


The  example  above  indicates  that  the  material  property  data  identified  with  an  id  parameter  of  8 
and  a  matname  of  ‘insulation’  should  be  read  from  the  input  file  set  by  the  &fnam 
materiaifile  parameter.  The  emissivity  is  set  to  0.91,  the  density  is  128  kg/m3,  and  the  heat 
capacity  is  a  constant  1,157  J/kg-K.  The  thermal  conductivity  varies  with  temperature  and  the 
temperature-conductivity  data  are  entered  on  two  &curv  lines  with  the  id  parameter  set  to 
‘ins_k\  The  &curv  lines  are  located  in  the  file  set  by  the  &fnam  materiaifile  parameter. 


4.4  Mechanical  Ventilation 

Forced  ventilation  conditions  are  simulated  using  duct  segments,  duct  nodes,  and  fans. 
Ducts  connect  two  nodes  and  are  assumed  uniform  in  terms  of  cross  section  and  friction  loss 
between  nodes.  Fans  are  added  to  duct  segments  and  are  essentially  represented  as  a 
discontinuous  pressure  change  that  is  a  function  of  the  local  pressure.  Mechanical  ventilation 
may  be  used  to  simulate  mechanical  supply,  exhaust,  or  other  types  of  inter-compartment 
connections  involving  the  transfer  of  air  through  ductwork. 


4.4.1  Mechanical  Ventilation  Nodes  (&rnod) 

The  &rnod  NAMELIST  group  is  used  to  enter  the  node  locations  of  duct  segments.  The 
nodes  are  used  to  connect  ducts  to  each  other,  to  compartments,  and  to  the  ambient  surroundings. 
The  &rnod  NAMELIST  group  Is  optional;  however,  at  least  two  srnod  lines  are  required  if  there 
is  a  srdct  NAMELIST  group  line. 

Nodes  should  be  defined  at  the  following  locations  within  a  ductwork  system: 

•  At  junctions  of  three  or  more  ducts. 

•  Between  multiple  fans 

•  At  compartment  or  ambient  terminal  locations. 
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In  addition,  nodes  may  also  be  used  where  ducts  change  shape,  cross-sectional  area,  or 
where  the  flow  changes  direction;  however,  it  is  not  necessary  to  do  so  provided  the  flow  loss 
coefficients  for  the  duct  and  all  associated  fittings  is  known. 

If  nodes  are  set  using  the  srnod  parameter  group,  then  three  parameters  are  required  to  be 
set  on  each  srnod  line: 


REQUIRED  PARAMETERS  WHEN  A  &RNOD  NAMELIST  GROUP  IS  USED 


id 

(Long  integer)  Identification  number  for  an  individual  &rnod  parameter 
group  line.  The  id  must  be  unique  among  other  srnod  parameter  id 
values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 
is  no  default  value. 

nducts 

(Integer)  The  number  of  ducts  connected  at  the  node.  If  the  node  is  a 
terminus  (terminal  set  to  .  true  . ),  then  nducts  must  be  set  to  2.  The 
maximum  number  of  ducts  that  may  be  connected  to  a  single  node  is  30. 
There  is  no  default  value. 

ducts 

(Long  integer  array)  An  array  of  long  integers,  separated  by  commas, 
containing  the  identification  numbers  (srdct  id  parameter  values)  of  the 
ducts  connected  at  the  node.  If  the  node  is  a  terminus  (terminal  set  to 
.  true  . ),  then  the  first  element  of  the  array  contains  the  compartment 
identification  number  (&comp  id  parameter  value)  or  a  -1  (ambient 
surroundings).  There  are  no  default  values  for  the  array  elements. 

There  is  no  requirement  for  the  &rnod  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  srnod  id 
parameters  is  acceptable.  The  ducts  connected  at  the  node  must  be  different;  the  same  duct 
should  not  be  multiply  connected  at  a  node. 

The  remaining  &rnod  parameters  are  optional: 


OPTIONAL  PARAMETERS 


nodename 

(Text  string)  Name  of  duct  node.  The  duct  node  name  is  used  primarily 
for  labeling  various  output  data.  The  maximum  number  of  characters  for 
the  duct  node  name  is  40.  The  default  is  ‘ductnode  id’,  where  id  is  the 
ductnode  identification  number,  if  terminal  is  .  false.,  and 
‘ductterminal-c  ducts  ( l )  ’,  where  ducts  ( 1 )  is  the  first  element  of 
the  ducts  array,  if  terminal  is  .  true  . 

elevation  " 

(Real  number)  The  elevation  of  the  node  above  the  reference  datum  (m). 
The  default  value  is  0.0  m. 
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terminal 


(Logical)  .  true  .  indicates  that  the  node  is  a  duct  terminus  (compartment 
or  ambient  surroundings)  and  .  false  .  indicates  that  the  node  is  not  a 
terminus.  The  default  state  is  .  false  . 


loss 


nosolve 


zone 


(Real  number  array)  An  array  of  real  numbers,  separated  by  commas,  that 
correspond  to  the  pressure  loss  coefficients  for  branch  flow  along  all 
possible  flow  paths  at  the  node.  The  number  of  elements  in  the  array  is 
equal  to  the  number  of  branch  flow  permutations  possible  at  the  node 
junction.  The  number  of  permutations  is  calculated  using: 

nekm=N-{N- 1) 

where  nejem  is  the  number  of  elements  in  the  loss  array  and  N  is  the 
number  of  ducts  connected  to  the  node  (nducts).  The  maximum  number 
of  elements  in  the  array  is  870  because  the  nducts  parameter  is  limited  to 
30.  The  sequence  of  the  array  is  based  on  the  order  of  the  srdct  id 
parameters  that  constitute  the  ducts  parameter  array.  The  order  is  as 
follows: 


Duct  that  flow  issues  from 

Duct  that  flow  enters  into 

loss  array 
element 

ducts  array  element  1 

ducts  army  element  2 

1 

2 

0 

0 

ducts  array  element  nducts 

ducts  army  element  2 

ducts  array  element  1 

ducts  array  element  3 

ducts  array  element  4 

ducts  array  element  nducts 

0 

■■Hi 

ducts  array  element  nducts 

D 

ducts  array  element 
nducts  —  1 

The  default  value  is  0,0  for  all  array  elements. 

(Logical)  .  true  .  indicates  that  the  mass  or  volume  flow  rate  through  the 
junction  is  fixed  (via  srdct  mdot  or  vdot  parameters)  and  .  false. 
indicates  that  the  mass  flow  rate  should  be  determined.  The  default  state 
is  .  FALSE . 

(Long  integer)  The  zone  number  where  the  compartment  is  located.  The 
default  is  calculated  by  FSSIM  and  the  value  entered  here  is  overridden 
regardless.  This  parameter  should  thus  not  be  used. 
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The  reference  datum  is  the  elevation  that  corresponds  to  the  sexec  pamb  and  tamb 
parameter  values.  The  pressure  loss  coefficients  are  determined  from  published  data  on  the 
particular  type  of  junction  (‘Y’  tee,  couple,  reducer,  cross,  or  ‘Y’  cross)  [Idelchik,  1986]. 
Typically,  the  pressure  loss  coefficient  is  a  function  of  the  Reynolds  number.  However,  as  with 
other  types  of  duct  fittings,  it  is  only  weakly  dependent  on  the  Reynolds  number  and  will  thus 
not  vary  significantly  over  a  large  range  of  Reynolds  number  values. 

There  may  be  circumstances  where  the  mass  or  volume  flow  rate  through  a  duct  is 
specified,  as  described  in  Section  4.4.1.  In  such  case,  setting  the  nosoive  parameter  to  .true. 
may  reduce  the  number  of  calculations  performed  and  possibly  decrease  the  computation  time 
for  a  simulation.  Note  that  in  order  to  conserve  mass,  where  two  or  more  duct  segments  ducts 
join  at  a  node,  the  net  mass  flow  must  balance.  As  previously  noted,  if  one  flow  path  is  left 
unspecified,  FSSIM  will  determine  the  appropriate  mass  flow  in  the  duct  segment.  In  this  case, 
the  nosoive  parameter  should  be  set  to  .  false  . 


An  example  of  a  srnod  NAMELIST  parameter  group  line  is  as  follows: 


&RNOD  id=3,  nducts=2,  terminals  TRUE . ,  ducts=-l,3,  elevation=l . 73,  loss=0.93, 
0.93,  nodename=' Duct  terminus  at  ambient'  / 


The  above  &rnod  line  is  identified  by  the  id  parameter  set  to  3  and  the  nodename 
parameter  set  to  ‘Duct  terminus  at  ambient’.  The  node  is  a  terminus  for  the  duct  with  a 
&rdct  id  parameter  of  3  at  the  ambient  surroundings.  The  elevation  of  this  node  is  1.73  m  above 
the  reference  datum.  The  pressure  loss  coefficient  of  0.93  applies  to  flow  in  both  directions. 


4.4.2  Ducts  (&RDCT) 

The  &rdct  NAMELIST  group  is  used  to  enter  the  location  and  physical  characteristics  of 
duct  segments.  The  srdct  NAMELIST  group  is  optional;  however,  it  is  an  essential  NAMELIST 
group  if  forced  ventilation  is  being  modeled.  A  duct  segment  is  located  between  two  nodes, 
defined  using  the  srnod  NAMLIST  group. 

If  one  or  more  ducts  are  set  using  the  &rdct  parameter  group,  then  two  parameters  are 
required  to  be  set  on  each  srdct  line: 
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REQUIRED  PARAMETERS  WHEN  A  &HDCT  NAMELIST  GROUP  IS  USED 


id 

(Long  integer)  Identification  number  for  an  individual  srdct  parameter 
group  line.  The  id  must  be  unique  among  other  &rdct  parameter  id 
values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 
is  no  default  value. 

node 

(Long  integer  array)  Array  of  two  &rnod  id  parameters  defining  the 
endpoints  of  the  duct  segment.  There  are  no  default  values  for  the  array 
elements. 

There  is  no  requirement  for  the  srdct  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  srdct  id 
parameters  is  acceptable.  The  endpoints  of  a  duct  segment  defined  by  the  node  array  must  use 
different  &rnod  id  parameter  values  (viz.,  a  duct  must  not  begin  and  end  at  the  same  location). 

The  remaining  srdct  parameters  are  optional: 


OPTIONAL  PARAMETERS 


duct name 

(Text  string)  Name  of  duct.  The  duct  name  is  used  for  primarily  labeling 
various  output  data.  The  maximum  number  of  characters  for  the  duct 
name  is  40.  The  default  is  ‘duct  id’,  where  id  is  the  duct  identification 
number. 

area 

(Real  number)  The  nominal  flow  area  of  the  duct  (m2).  The  default  value 
is  1.0  m2. 

perimeter 

(Real  number)  The  wetted  perimeter  of  the  duct  segment  (m).  There  is  no 

default  value. 

hd 

(Real  number)  The  hydraulic  radius  of  the  duct  segment  (m).  If 
perimeter  parameter  is  undefined,  then  the  default  hydraulic  radius  is 
■Ja  ,  where  A  is  the  equal  to  the  area  parameter  value.  Otherwise,  the 
default  hydraulic  perimeter  is: 

4  A 

P 

where  A  is  equal  to  the  area  parameter  value  and  P  is  equal  to  the 
perimeter  value. 
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kloss 

(Real  number  array)  An  array  of  numbers,  separated  by  commas, 
containing  the  turbulent  pressure  loss  coefficients  for  the  duct  segment 
and  any  fittings  that  may  be  located  within  the  duct  segment.  The  first 
element  in  the  array  is  the  coefficient  corresponding  to  positive  flow  and 
the  second  element  corresponds  to  negative  (reverse)  flow.  Positive  flow 
is  defined  as  flow  from  the  first  element  in  the  node  array  to  the  second 
element  and  vice  versa  for  reverse  flow.  The  pressure  loss  coefficient 
includes  the  pressure  losses  associated  with  flow  through  bends  and 
fittings.  It  may  also  include  the  friction  losses  to  the  duct  boundaries.  If 
boundary  friction  losses  are  included  in  this  parameter,  then  the 
roughness  parameter  should  be  set  to  0.0.  The  default  value  is  0.0. 

length 

(Real  number)  The  length  of  the  duct  segment  (m).  The  default  value  is 

1.0  m. 

roughness 

(Real  number)  The  surface  roughness  of  the  duct  segment  (m).  The 
default  value  is  0.00091  m.  If  the  friction  flow  losses  for  the  duct  are 
included  in  the  kloss  parameter,  roughness  should  be  set  to  0.0  m. 

damper 

(Long  integer)  Control  function  sctrl  id  used  to  open  or  close  a  damper 
in  a  duct  segment.  A  value  of -1  indicates  that  there  is  no  damper  in  the 
duct.  The  default  value  is  -1. 

fan 

(Long  integer)  Used  to  specify  that  a  fan  is  located  between  the  duct 
segment  endpoints  as  defined  by  the  first  element  in  the  node  array  to  the 
second  element.  The  value  is  equal  to  a  srfan  id  value  for  a  &rfan  line 
containing  the  appropriate  fan  information.  A  value  of-1  indicates  that 
there  is  no  fan  associated  with  this  duct  segment.  The  default  value  is  -1 . 

mdot 

(Real  number)  Mass  flow  rate  through  duct  (kg/s).  Positive  flow  is  from 
the  first  element  in  the  node  array  to  the  second  element.  The  mass  flow 
is  a  forced  condition  and  should  only  be  used  when  there  is  an  absence  of 
information  regarding  the  ductwork.  The  mass  flow  rate  must  be  greater 
than  -999  kg/s.  A  value  less  than  or  equal  to  -999  kg/s  disables  this 
feature.  The  default  value  is  -999  kg/s. 

vdot 

(Real  number)  Volume  flow  rate  through  duct  (mVs).  Positive  flow  is 
from  the  first  element  in  the  node  array  to  the  second  element.  The 
volume  flow  is  a  forced  condition  and  should  only  be  used  when  there  is 
an  absence  of  information  regarding  the  ductwork.  The  volume  flow  rate 
must  be  greater  than  -999  m3/s.  A  value  less  than  or  equal  to  -999  m3/s 
disables  this  feature.  The  default  value  is  -999  m3/s. 

zone 

(Long  integer)  The  zone  number  where  the  compartment  is  located.  The 
default  is  calculated  by  FSSIM  and  the  value  entered  here  is  overridden 
regardless.  This  parameter  should  thus  not  be  used. 

round 

(Logical)  .true,  if  the  cross-section  of  the  duct  is  round  and  .false,  if 
the  cross  section  is  not  round.  The  default  state  is  .  true  .  This  parameter 
is  currently  not  used  by  the  mode. 
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The  perimeter  variable  is  only  used  to  calculate  the  hydraulic  diameter,  hd.  It  is  thus  not 
necessary  to  specify  this  parameter  if  the  value  of  hd  is  known. 

The  pressure  loss  coefficient  is  determined  by  summing  the  contribution  from  the  various 
fittings  that  are  between  duct  segment  end  points,  and  may  include  friction  pressure  losses  to  the 
duct  boundary.  Typically,  the  pressure  loss  coefficient  is  a  function  of  the  Reynolds  number  and 
is  tabulated  for  the  type  of  fitting  present  [e.g.,  Crane  Company,  1991;  Idelehik,  1986;  Brater, 
1996].  In  many  cases,  the  pressure  loss  coefficient  is  only  weakly  dependent  on  the  Reynolds 
number  and  will  thus  not  vary  significantly  over  a  large  range  of  Reynolds  numbers.  FSSIM  will 
automatically  calculate  the  friction  pressure  loss  as  a  function  of  the  Reynolds  number  if  a  duct 
roughness  is  provided.  The  roughness  is  material  specific  and  may  be  estimated  from  various 
sources  [Crane  Company,  1991;  Idelehik,  1986;  Brater,  1996]. 

Duct  segments  may  contain  a  single  fan  with  characteristics  defined  at  the  referenced 
&rfan  NAMELIST  line.  Fans  are  located  midway  between  duct  segment  endpoints  in  the  model 
and  are  essentially  a  discontinuous  pressure  change  within  the  ductwork. 

If  the  duct  contains  a  damper,  it  may  opened  or  closed  based  on  an  input  condition  (time, 
temperature,  etc.)  using  a  sctrl  function. 

There  may  be  circumstances  where  there  is  inadequate  information  regarding  the 
ductwork  (length,  friction  loss,  and  fitting  losses).  In  this  case,  it  is  possible  to  force  a  mass  or 
volume  flow  rate  through  a  duct  segment  using  the  mdot  or  vdot  parameters.  When  either  is  set, 
the  pressure  loss  calculations  are  not  performed  and  this  flow  rate  is  maintained  regardless  of  the 
duct  segment  endpoint  conditions.  Note  that  in  order  to  conserve  mass,  where  two  or  more  duct 
segments  join  at  a  node,  the  net  mass  flow  must  balance.  If  either  feature  is  used,  then  one  flow 
path  should  remain  unspecified  such  that  the  flow  is  calculated  within  FSSIM  to  avoid  creating  a 
situation  where  mass  is  not  balanced.  Both  the  mass  flow  and  the  volume  flow  must  not  be  set 
within  a  single  duct;  otherwise,  FSSIM  will  generate  an  error  and  terminate  the  simulation. 

An  example  of  a  srdct  NAMELIST  parameter  group  line  and  the  associated  sctrl  line 
is  as  follows: 


SRDCT  id=3,  node=l,2,  area=0.607,  length=3.37,  kloss=3.89,  hd=0.82, 
roughness=0. 00055,  damper=2,  ductname=' Duct  Segment  A'  / 

SCTRL  id=2,  ctype=l,  inputvar=3,  location=3,  initial=. FALSE. ,  setpoint=76. 0  / 

The  above  srdct  line  is  identified  by  the  id  parameter  set  to  3  and  the  ductname 
parameter  set  to  5  Duct  Segment  a’.  The  duct  segment  endpoints  are  defined  by  nodes  with  the 
srnod  id  parameters  set  to  1  and  2.  The  length  of  the  duct  is  3.37  m,  the  nominal  flow  area  is 
0,607  m2,  and  the  hydraulic  diameter  is  0,82  m.  The  pressure  loss  coefficient  of  3.89,  which 
could  include  a  valve  and  an  elbow,  but  does  not  include  friction  losses  at  the  duct  boundary. 
The  latter  are  calculated  by  FSSIM  using  a  duct  surface  roughness  of  0.00055  m.  There  is  a 
damper  with  actuation  criteria  input  using  a  sctrl  line  with  an  id  parameter  set  to  2.  The 
referenced  sctrl  line  forces  the  damper  to  close  when  the  temperature  within  the  duct  segment 
with  an  id  parameter  set  to  3  (the  referencing  duct)  reaches  76  °C. 
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4.4.3  Fans  (&rfan) 


The  &rfan  NAMELIST  group  is  used  to  specify  the  fan  characteristics  (performance 
curve,  control  functions,  and  pressure  loss  coefficients).  The  srfan  NAMELIST  group  is 
optional;  however,  it  is  an  essential  NAMELIST  group  if  forced  ventilation  is  being  modeled. 

Fans  are  simulated  within  a  duct  segment  defined  by  a  &rdct  line.  If  a  &rfan  line  is 
included  in  the  model,  then  two  parameters  are  required  to  be  set  on  each  srfan  line: 


REQUIRED  PARAMETERS  WHEN  A  &RFAN  NAMELIST  GROUP  IS  USED 


id 

(Long  integer)  Identification  number  for  an  individual  srfan  parameter 
group  line.  The  id  must  be  unique  among  other  &rfan  parameter  id 
values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 
is  no  default  value. 

control 

(Long  integer)  Control  function  sctrl  id  used  to  turn  a  fan  on  or  off. 
There  is  no  default  value. 

There  is  no  requirement  for  the  &rfan  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  &rfan  id 
parameters  is  acceptable. 

The  remaining  srfan  parameters  are  optional: 


OPTIONAL  PARAMETERS 


fanname 

(Text  string)  Name  of  the  fan.  The  fan  name  is  used  primarily  for 
labeling  various  output  data.  The  maximum  number  of  characters  for  the 
fan  name  is  40.  The  default  is  ‘fan  id’,  where  id  is  the  fan  identification 
number. 

fileread 

(Logical)  .  true  .  indicates  that  the  fan  data  should  be  read  from  the  file 
set  by  the  &fnam  f  anf  ile  parameter  and  .  false  .  indicates  that  the  fan 
data  are  contained  within  the  srfan  id  line  located  in  the  main  input  file. 
The  default  state  is  .  true  .  If  the  fanname  parameter  is  not  set  (i.e.,  the 
default  is  used),  then  the  fan  data  will  not  be  read  from  the  file  set  by  the 
&fnam  fanfile  parameter. 
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ftype 

(Short  integer)  The  type  of  fan  in  terms  of  pressure  versus  flow 
characteristics: 

1  -  The  fan  produces  a  constant  pressure  increase  equal  to  delp 

(Pa)  regardless  of  the  flow. 

2  -  The  fan  produces  a  pressure  increase  that  is  a  quadratic 

function  of  the  flow.  The  flow  at  zero  pressure  (full)  and 
pressure  at  zero  flow  (shutoff)  are  provided. 

3  -  The  fan  flow  versus  pressure  data  are  provided  in  tabular 

format  using  scurv  lines  with  an  id  parameter  set  to  the  text 
string  table.  The  maximum  number  of  pressure  versus  flow 
data  points  is  32,000. 

The  default  value  is  1. 

delp 

(Real  number)  The  pressure  increase  (Pa)  at  the  fan  when  ftype  is  set  to 

1.  The  default  value  is  0.0  Pa. 

full 

(Real  number)  The  volumetric  flow  (m3/s)  at  zero  pressure  when  ftype 
is  set  to  2.  The  default  value  is  1 .0  m3/s. 

shutoff 

(Real  number)  The  pressure  (Pa)  at  zero  flow  when  ftype  is  set  to  2.  The 
default  value  is  1,000  Pa. 

table 

(Text  string)  The  name  of  the  scurv  line(s)  that  contain  the  flow  vs.  fan 
pressure  tabular  data.  There  is  no  default  and  the  maximum  number  of 
characters  is  40.  This  parameter  is  not  used  if  ftype  is  set  to  1  or  2. 

reverse 

(Logical)  .  true  .  indicates  that  the  fan  pressure  versus  flow  is  positive 
in  the  opposite  direction  of  positive  flow  and  .  false  indicates  that  the 
fan  pressure  versus  flow  is  positive  in  the  direction  of  positive  flow. 
Positive  flow  is  defined  as  flow  from  the  first  element  in  the  &rdct  node 
array  to  the  second  element  in  the  srdct  node  array  of  the  &rdct  line 
that  references  this  fan.  The  default  state  is  .  false  . 

kloss 

(Real  number)  The  pressure  loss  coefficient  associated  with  the  fan  when 
the  fan  is  not  operating.  The  default  is  1 .0. 

The  &rfan  variables  may  be  read  from  a  separate  input  file  as  specified  by  the  &fnam 
f  anf  ile  parameter.  The  method  is  identical  to  that  described  for  the  scmpn  NAMELIST 
parameter  group  in  Section  4.3.4  and  is  not  repeated  here. 

If  the  velocity  versus  pressure  information  is  entered  using  &curv  lines  (viz.,  f  type  is  set 
to  3),  then  the  profile  must  vary  with  pressure  and  it  should  resemble  a  typical  fan  curve  for  a 
centrifugal  or  axial  fan  (See  Figure  4-1).  The  pressure  is  assumed  constant  for  flow  rates  that  are 
above  or  below  the  flow  range  set  by  the  scurv  lines.  If  the  flow  is  less  than  the  minimum  flow 
value,  then  the  pressure  is  assumed  equal  to  the  pressure  specified  for  the  minimum  flow.  If  the 
flow  exceeds  the  maximum  flow,  then  the  pressure  is  assumed  equal  to  zero. 

The  pressure  loss  coefficient  for  the  fan  is  applicable  only  when  the  fan  is  not  operating 
but  there  is  flow  through  the  fan.  In  this  case,  the  fan  is  treated  as  a  fitting.  The  pressure  loss 
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coefficient  is  a  function  of  the  particular  fan  type  and  the  Reynolds  number.  The  coefficient  may 
be  available  directly  from  the  manufacturer.  Estimates  could  be  made  from  information  in  some 
references  [ e.g Crane  Company,  1991;  Idelchik,  1986;  Brater,  1996]. 

An  example  of  a  srfan  NAMELIST  parameter  group  line  and  the  associated  &curv  and 
&ctrl  lines  is  as  follows: 


&RFAN  id=6,  fileread=. FALSE . ,  control=l,  ftype=3,  table=' Supply  Fan  Data', 
reverse=. FALSE,  kloss=1.13,  fanname=' Supply  fan'  / 

&CTRL  id=l,  ctype=l,  inputvar=l,  initials  FALSE. ,  setpoint=120 . 0  / 

&CURV  id=' Supply  Fan  Data',  xy=0.0,  96.5  / 

&CURV  id=' Supply  Fan  Data',  xy=1.41,  49.7  / 

&CURV  id=' Supply  Fan  Data',  xy=1.99,  62.2  / 

&CURV  id=' Supply  Fan  Data',  xy=2.3,  5.2  / 


The  above  &rfan  line  is  read  from  the  primary  input  file  and  is  identified  by  the  f  type 
parameter  set  to  3  and  the  text  string  ‘supply  fan’.  The  fan  is  initially  off  but  is  turned  on  after 
120  seconds  from  the  start  of  the  simulation.  The  pressure  loss  coefficient  is  1.13  when  the  fan  is 
not  operating  and  the  pressure  versus  flow  data  are  provided  on  &curv  lines  with  the  id 
parameter  set  to  ‘supply  Fan  Data’. 


Figure  4-1.  Typical  Centrifugal  and  Axial  Fan  Curves  [Hayashi  et  al.,  1985;  ASHRAE,  1991] 
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4.5  Source  Fires  and  Fuel  Load 


Compartment  fires  are  either  initially  specified  in  terms  of  a  heat  release  rate  (power) 
versus  time  profile  or  arise  in  a  given  fuel  load  when  specific  conditions  (temperature)  are  met. 
An  initial  fire  is  set  using  the  &fire  parameter  group.  The  fuel  load  in  a  compartment  fuel  load  is 
defined  using  the  suses  NAMELIST  group  and  the  ignition  criteria  for  this  fuel  load  is  set  using 
the  sfuel  parameter  group, 

A  compartment  does  not  require  a  fire  or  a  fuel  load  for  FSSIM  to  successfully  run.  There 
are  some  applications,  such  as  verifying  the  initial  ventilation  conditions,  where  it  is  useful  to  run 
FSSIM  without  including  a  source  fire. 


4.5.1  Initial  Source  Fire  (&fire) 

The  &fire  NAMELIST  group  is  used  to  specify  the  initial  heat  release  rate  (power)  and 
species  yield  fractions  of  a  compartment  fire.  A  particular  compartment  should  not  have  more 
than  one  initial  fire  specified.  A  simulation  maybe  performed  without  setting  a  sfire  parameter 
group.  Such  a  model  may  be  useful  when  verifying  the  initial  or  ambient  ventilation  conditions 
in  compartments. 

There  is  one  parameter  that  must  be  set  if  the  &fire  line  is  used; 


REQUIRED  PARAMETER  WHEN  A  fiFIRE  NAMELIST  GROUP  IS  USED 

(Long  integer)  Identification  number  for  an  individual  &fire  parameter 
group  line.  The  id  must  be  unique  among  other  &fire  parameter  id 
"  d  values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 

is  no  default  value. 


There  is  no  requirement  for  the  sfire  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  sfire  id 
parameters  is  acceptable. 

The  remaining  &fire  parameters  are  optional: 
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OPTIONAL  PARAMETERS 


ftype 

(Short  integer)  The  type  of  heat  release  rate  versus  time  profile  that 
defines  the  initial  source  fire: 

1  -  The  heat  release  rate  is  constant. 

2  -  The  heat  release  rate  is  proportional  to  time  squared  (t2)  up 

to  the  time  set  by  t2_time  or  the  heat  release  rate  set  by 
t2_hrr,  after  which  the  fire  is  constant. 

3  -  The  heat  release  rate  versus  time  is  provided  in  tabular 

format  using  scurv  lines  with  an  id  parameter  set  to  the 
text  string  table.  The  maximum  number  of  heat  release 
rate  versus  time  data  points  is  32,000. 

The  default  value  is  1 . 

f subtype 

(Short  integer)  Fire  sub-type  parameter.  If  ftype  is  set  to  1,  then 
f  subtype  is  the  heat  release  rate  selector: 

l-100kW. 

2  —  250  kW. 

3- 1,000  kW. 

4  -  constant  heat  release  rate  set  by  the  hrr  variable. 

If  ftype  is  set  to  2,  then  f  subtype  is  the  growth  rate  constant  selector: 

1  -  0.02778  kW/s2  (‘slow’  growth  rate). 

2  -  0.01  111  kW/s2  (‘medium’  growth  rate). 

3  -  0.0444  kW/s2  (‘fast’  growth  rate). 

4  -  0.1778  kW/s2  (‘ultra  fast’  growth  rate). 

5  -  growth  rate  set  by  the  growth  parameter. 

The  default  value  is  1  regardless  of  the  value  of  ftype.  This  parameter  is 
not  used  if  ftype  is  not  set  to  1  or  2. 

location 

(Long  integer)  The  compartment  identification  (scomp  id  parameter) 
where  the  fire  is  located.  The  default  value  is  1 . 

hrr 

(Real  number)  The  heat  release  rate  of  the  source  fire  (kW).  This 
parameter  is  only  used  if  ftype  is  set  to  1  and  f  subtype  is  set  to  4.  The 
default  value  is  100  kW. 

growth 

(Real  number)  The  growth  rate  constant  for  a  t2  source  fire  (kW/s2).  This 
parameter  is  only  used  if  ftype  is  set  to  2  and  f  subtype  is  set  to  5.  There 
is  no  default  value  unless  f  subtype  is  not  properly  set  (/. e. ,  has  a  value 
greater  than  5  or  less  than  1),  in  which  case  the  default  growth  rate  is 
0.01173  kW/s2. 

t2_time 

(Real  number)  The  maximum  growth  period  for  a  t2  source  fire  (s).  This 
parameter  is  only  used  if  ftype  is  set  to  2.  The  default  value  is  1  x  105  s. 
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t2_hrr 

(Real  number)  The  maximum  heat  release  rate  for  a  tl  source  fire  (kW). 
This  parameter  is  only  used  if  ftype  is  set  to  2.  The  default  value  is  1  x 

109  kW. 

t_start 

(Real  number)  The  time  from  the  simulation  start  that  the  initial  source 
fire  begins  (s).  The  default  value  is  0.0  s. 

t_end 

(Real  number)  The  time  from  the  simulation  start  that  the  initial  source 
fire  ends  (s).  Note:  the  fire  may  end  earlier  due  to  suppression  or  oxygen 
consumption.  The  default  value  is  1  x  105  s. 

table 

(Text  string)  The  name  of  the  &curv  line(s)  that  contain  the  fire  heat 
release  rate  (power)  versus  time  (measured  from  t_start)  tabular  data. 
Time  (s)  is  the  first  element  of  the  &curv  xy  array  and  heat  release  rate 
(kW)  is  the  second  element.  There  is  no  default  and  the  maximum 
number  of  characters  is  40.  This  parameter  is  not  used  if  ftype  is  set  to  1 
or  2. 

tpyro 

(Real  number)  The  pyrolysis  temperature  (K).  The  default  value  is  370.0 

K. 

deltah  f 

(Real  number)  The  heat  of  combustion  of  the  fuel  material  (J/kg).  The 
default  value  is  44.6  x  106  J/kg. 

deltahjp 

(Real  number)  The  heat  of  pyrolysis  of  the  fuel  material  (J/kg).  The 
default  value  is  3.65  x  105  J/kg. 

chi  r 

(Real  number)  The  fraction  of  source  fire  heat  release  rate  energy  that  is 
released  as  thermal  radiation.  The  default  value  is  0.35. 
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yco2 

(Real  number)  The  carbon  dioxide  yield  (kg  carbon  dioxide  produced  per 
kg  fuel  consumed).  The  default  value  is  calculated  using  a  chemical  yield 
balance  and  depends  on  whether  or  not  the  oxygen  yield  is  defined.  If  the 
oxygen  yield  is  not  defined,  then  the  default  carbon  dioxide  yield  is 
calculated  whether  or  not  a  value  is  provided.  The  value  is  determined  as 
follows: 

11  11  11  11 

y  co2  ^  ry  y co  ^  soot  27  ^ a lO 

where  yCo,  is  the  yield  fraction  of  carbon  dioxide  (yco2  [kg/kg]),  yCo  is 
the  yield  fraction  of  carbon  monoxide  (yco  [kg/kg]),  ysoot  is  the  yield 
fraction  of  soot  (ysoot  [kg/kg]),  and  yn,o  is  the  yield  fraction  of  water 
(yh2o  [kg/kg]).  If  the  oxygen  yield  is  defined,  then  the  default  carbon 
dioxide  yield  is  only  calculated  if  a  value  is  not  provided  using  the 
following: 

-11  _n  11 

yco2  g  g  y n |  ^  yco 

where  y0,  is  the  yield  fraction  of  oxygen  (yo2  [kg/kg]).  To  avoid  potential 
rounding  error  in  the  model,  it  is  recommended  that  this  parameter  be 
calculated  by  FSSIM. 

yh2o 

(Real  number)  The  water  yield  (kg  water  produced  per  kg  fuel 
consumed).  The  default  value  is  1.596  kg/kg. 

yco 

(Real  number)  The  carbon  monoxide  yield  (kg  carbon  monoxide 
produced  per  kg  fuel  consumed).  The  default  value  is  0.012  kg/kg. 

yo2 

(Real  number)  The  oxygen  yield  (kg  oxygen  consumed  per  kg  fuel 
burned).  The  default  value  (undefined  condition)  is  calculated  using  a 
chemical  yield  balance  assuming  the  fuel  does  not  contain  significant 
amounts  of  oxygen: 

8  8  4 

y°2  ~  |  j  yco2  +  g  yHlo + y  yco 

where  yo,  is  the  yield  fraction  of  oxygen  (yo2  [kg/kg]),  yco,  is  the  yield 
fraction  of  carbon  dioxide  (yco2  [kg/kg]),  yHp  is  the  yield  fraction  of 
water  (yh2o  [kg/kg]),  and  yco  is  the  yield  fraction  of  carbon  monoxide 
(yco  [kg/kg]). 

ysoot 

(Real  number)  The  soot  yield  (kg  soot  produced  per  kg  fuel  consumed). 
The  default  value  is  0.042  kg/kg. 

Nearly  any  type  of  source  fire  heat  release  rate  profile  may  be  specified  in  FSSIM  using 
the  available  ftype  and  f  subtype  options.  The  particular  fire  will  depend  on  the  fuel  package 
burning  properties,  orientation,  and  surface  area.  References  and  approaches  for  estimating  the 
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heat  release  rate  are  found  in  various  sources  including  Lee  [1985],  Babrauskas  [2002],  Sardqvist 
[1993],  Nowlen  [1986],  NFPA  72  [2002],  and  Madrzykowski  [1996]. 

The  heat  release  rate  is  assumed  constant  for  times  that  are  above  or  below  the  time  range 
set  by  the  &curv  lines.  If  the  time  is  less  than  the  minimum  time  value,  then  the  heat  release  rate 
is  assumed  equal  to  the  heat  release  rate  specified  for  the  minimum  time.  The  reverse  is  true  if 
the  time  is  greater  than  the  maximum  time  value. 

The  fuel  characteristics  (heat  of  combustion,  pyrolysis  temperature,  the  heat  of  pyrolysis, 
and  the  species  yields)  are  constant  over  the  course  of  the  simulation  and  are  material  specific. 
Considerable  care  is  necessary  when  entering  the  yield  data  to  prevent  violations  of  species 
conservation.  Of  particular  note  is  the  situation  when  FSSIM  is  used  to  calculate  the  oxygen 
yield  and/or  the  carbon  dioxide  yield.  In  both  cases,  it  is  possible  that  unrealistic  (negative) 
yields  are  calculated  if  the  yield  of  water,  carbon  monoxide,  and  soot  are  over-specified.  Sources 
for  estimating  the  value  of  the  yield  parameters  include  Tewarson  [2002],  Babrauskas  [2003], 
and  Sardqvist  [1993]. 

An  example  of  a  sfire  NAMELIST  parameter  group  line  is  as  follows: 

&FIRE  id=12,  location=2»  ftype=2,  fsubtype=3,  t2_hrr=1500 . 0,  t_start=0.0, 
deltah_f=4 4 100000 . 0,  yh20=1.0,  yco=0 .  03 ,  ysoot=0.07  / 

The  above  &fire  line  is  identified  with  the  id  parameter  set  to  12,  is  located  in  the 
compartment  with  the  &comp  id  parameter  of  2,  and  begins  when  the  simulation  starts.  The  fire 
has  a  growth  rate  that  is  proportional  to  the  square  of  time  and  has  a  ‘fast’  heat  release  rate 
growth  constant.  The  heat  release  rate  becomes  constant  when  it  reaches  1,500  kW.  The  fuel  has 
a  carbon  monoxide  yield  of  0.03  kg/kg,  a  soot  yield  of  0.07  kg/kg,  and  a  water  yield  of  1 .0  kg/kg. 
The  carbon  dioxide  and  oxygen  yields  are  determined  by  FSSIM.  Based  on  the  equations  cited 
above,  the  yields  are  about  2.96  kg  CCL/kg  fuel  and  3.08  kg  02/kg  fuel,  respectively. 


4.5.2  Fuel  Load  Parameters  (&uses) 

The  &uses  parameter  group  is  used  to  set  the  fuel  load  and  to  associate  a  compartment  or 
compartments  with  a  &fuel  parameter  group  via  the  suses  id  parameter.  This  parameter  group 
is  optional;  however,  if  it  is  used  there  are  two  parameters  that  must  be  set: 


REQUIRED  PARAMETERS  WHEN  A  &USES  NAMELIST  GROUP  IS  USED 


id 

(Long  integer)  Identification  number  for  an  individual  suses  parameter 
group  line.  The  id  must  be  unique  among  other  suses  parameter  id 
values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 
is  no  default  value. 

fuel 

(Long  integer)  &fuel  id  parameter  containing  the  fuel  properties  and 
ignition  criteria.  There  is  no  default  value. 
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There  is  no  requirement  for  the  &uses  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  &uses  id 
parameters  is  acceptable.  A  &fuel  line  may  be  used  by  multiple  &uses  lines  (there  is  no 
requirement  for  the  fuel  parameter  to  be  unique  among  suses  lines). 

The  remaining  &uses  parameters  are  optional.  Many  are  the  same  type  of  parameter  as 
were  specified  for  the  &uses  parameter  group.  The  following  is  a  complete  list  of  the  optional 
parameters  that  may  be  set  on  the  &uses  line: 


OPTIONAL  PARAMETERS 


usename 

(Text  string)  Name  of  the  fuel  load.  The  fuel  load  name  is  used  primarily 
for  labeling  various  output  data.  The  maximum  number  of  characters  for 
the  fuel  load  name  is  40.  The  default  is  ‘usetype  id’,  where  id  is  the 
fuel  load  identification  number. 

fileread 

(Logical)  .  true  .  indicates  that  the  fuel  load  parameters  should  be  read 
from  the  file  set  by  the  &fnam  usefile  parameter  and  .  false  .  indicates 
that  the  fuel  load  data  are  contained  within  the  suses  id  line  located  in 
the  main  input  file.  The  default  state  is  .true.  If  the  usename  parameter 
is  not  set  (i.e.,  the  default  is  used),  then  the  uses  data  will  not  be  read 
from  the  file  set  by  the  sfnam  usefile  parameter. 

loading 

(Real  number)  The  combustible  fuel  load  (kg/m2).  The  default  value  is 

1 .0  kg/m2. 

The  &uses  variables  may  be  read  from  a  separate  input  file  as  specified  by  the  &fnam 
usefile  parameter.  The  method  is  identical  to  that  described  for  the  scmpn  NAMELIST 
parameter  group  in  Section  4.3.4  and  is  not  repeated  here. 

An  example  of  a  suses  NAMELIST  parameter  group  line  is  as  follows: 

&USES  ID=2,  fuel=6,  fileread-. FALSE. ,  loading=3.6,  usename=' Class  A  Fuel 
Load'  / 


The  above  &uses  line  is  read  from  the  primary  input  file  and  is  identified  with  the  id 
parameter  set  to  2  and  the  text  string  ‘class  A  Fuel  Load’.  The  fuel  combustion  properties  and 
ignition  criteria  are  defined  by  the  &fuel  line  with  the  id  parameter  of  6.  The  fuel  load  is  set  to 
3.6  kg/m2. 


4.5.3  Fuel  Load  (&fuel) 

The  &fuel  NAMELIST  group  is  used  to  specify  the  composition  and  ignition  criteria  for 
a  fuel  load  within  a  space.  The  information  provided  in  this  NAMELIST  group  is  generic  insofar 
as  it  is  not  placed  in  any  specific  location.  The  fuel  is  assigned  to  one  or  more  compartments 
using  the  &uses  id  parameter. 
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This  parameter  group  is  optional;  however,  if  it  is  used  there  is  one  parameter  that  must 

be  set: 


REQUIRED  PARAMETER  WHEN  A  &FUEL  NAMELIST  GROUP  IS  USED 

(Long  integer)  Identification  number  for  an  individual  &fuel  parameter 
group  line.  The  id  must  be  unique  among  other  &  fuel  parameter  id 
xu  values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 

is  no  default  value. 


There  is  no  requirement  for  the  sfuel  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  &fuel  id 
parameters  is  acceptable. 

The  remaining  sfuel  parameters  are  optional.  Many  are  the  same  type  of  parameter  as 
were  specified  for  the  sfire  parameter  group.  The  following  is  a  complete  list  of  the  optional 
parameters  that  may  be  set  on  the  sfuel  line: 


OPTIONAL  PARAMETERS 


fuelname 

(Text  string)  Name  of  fuel.  The  fuel  name  is  used  primarily  for  labeling 
various  output  data.  The  maximum  number  of  characters  for  the  fuel 
name  is  40.  The  default  is  ‘fuel  id5,  where  id  is  the  fan  identification 
number. 

f ileread 

(Logical)  .  true  .  indicates  that  the  fuel  data  should  be  read  from  the  file 
set  by  the  &fnam  f  ueif  ile  parameter  and  .  false  .  indicates  that  the  fuel 
data  are  contained  within  the  &fuel  id  line  located  in  the  main  input  file. 
The  default  state  is  .  true  .  If  the  fuelname  parameter  is  not  set  (i.e.,  the 
default  is  used),  then  the  fuel  data  will  not  be  read  from  the  file  set  by  the 
&FNAM  fuelfile  parameter. 

growth 

(Real  number)  The  heat  release  rate  growth  rate  (kW/s2).  The  default 
value  is  0.01 1 1 1  kW/s2  (‘medium’ t2  growth  rate). 

pyro  max 

(Real  number)  Maximum  fuel  pyrolysis  rate  (kg/m2).  The  default  Value  is 
0.01  kg/m2. 

tpyro 

(Real  number)  The  pyrolysis  temperature  (K).  The  default  value  is  370.0 

K. 

deltah_f 

(Real  number)  The  heat  of  combustion  of  the  fuel  material  (J/kg).  The 
default  value  is  44.6  x  1 06  J/kg. 

deltah  p 

(Real  number)  The  heat  of  pyrolization  of  the  fuel  material  (J/kg).  The 
default  value  is  3.65  x  105  J/kg. 
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chi  r 

(Real  number)  The  fraction  of  source  fire  heat  release  rate  energy  that  is 
released  as  thermal  radiation.  The  default  value  is  0.35. 

(Real  number)  The  carbon  dioxide  yield  (kg  carbon  dioxide  produced  per 
kg  fuel  consumed).  The  default  value  is  calculated  using  a  chemical  mass 
balance  and  depends  on  whether  or  not  the  oxygen  yield  is  defined.  If  the 
oxygen  yield  is  not  defined,  then  the  default  carbon  dioxide  yield  is 
calculated  whether  or  not  a  value  is  provided.  The  value  is  determined  as 
follows: 

_11  11  11  11 
y co :  ~  ^  y^  ^  y s,w<  27  y n 2° 

yco2 

where  yco,  is  the  yield  fraction  of  carbon  dioxide  (yco2  [kg/kg]),  yCo  is 
the  yield  fraction  of  carbon  monoxide  (yco  [kg/kg]),  ysoot  is  the  yield 
fraction  of  soot  (ysoot  [kg/kg]),  and  ynp  is  the  yield  fraction  of  water 
(yh2o  [kg/kg]).  If  the  oxygen  yield  is  defined,  then  the  default  carbon 
dioxide  yield  is  only  calculated  if  a  value  is  not  provided  using  the 
following: 

11  11  11 

yco2  -  8  yo2  9  yH2o  14  yco  t 

where  y0i  is  the  yield  fraction  of  oxygen  (yo2  [kg/kg]).  To  avoid  potential 
rounding  error  in  the  model,  it  is  recommended  that  this  parameter  be 
calculated  by  FSSIM. 

yh2o 

(Real  number)  The  water  yield  (kg  water  produced  per  kg  fuel 
consumed).  The  default  value  is  1.596  kg/kg. 

yco 

(Real  number)  The  carbon  monoxide  yield  (kg  carbon  monoxide 
produced  per  kg  fuel  consumed).  The  default  value  is  0.012  kg/kg. 

(Real  number)  The  oxygen  yield  (kg  oxygen  consumed  per  kg  fuel 
burned).  The  default  value  (undefined  condition)  is  calculated  using  a 
chemical  yield  balance  assuming  the  fuel  does  not  contain  significant 
amounts  of  oxygen: 

yo2 

8  8  4 

y°2  ~  j  j  fco2  +  g  y h2o  +  ^  Tco 

where  yo;  is  the  yield  fraction  of  oxygen  (yo2  [kg/kg]),  yCo;  is  the  yield 
fraction  of  carbon  dioxide  (yco2  [kg/kg]),  yHp  is  the  yield  fraction  of 
water  (yh2o  [kg/kg]),  and  yco  is  the  yield  fraction  of  carbon  monoxide 
(yco  [kg/kg]). 


ysoot 

(Real  number)  The  soot  yield  (kg  soot  produced  per  kg  fuel  consumed). 
The  default  value  is  0.042  kg/kg. 

ignite_ST 

(Real  number)  Non-overhead  surface  temperature  ignition  criteria  for 
fuel  (K).  The  default  value  is  473.0  K. 
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ignite  OV 

(Real  number)  Overhead  surface  temperature  ignition  criteria  for  fuel 
(K).  The  default  value  is  650.0  K. 

ignite_GT 

(Real  number)  Compartment  gas  temperature  ignition  criteria  for  fuel 
(K).  The  default  value  is  473.0  K. 

ignite  VT 

(Real  number)  Inflowing  gas  temperature  ignition  criteria  for  fuel  (K). 

The  default  value  is  600.0  K. 

The  &fuel  variables  may  be  read  from  a  separate  input  file  as  specified  by  the  &fnam 
fuel  file  parameter.  The  method  is  identical  to  that  described  for  the  scmpn  NAMELIST 
parameter  group  in  Section  4.3.4  and  is  not  repeated  here. 

As  with  the  &fire  parameter  group,  the  fuel  characteristics  (heat  of  combustion, 
pyrolysis  temperature,  the  heat  of  pyrolization,  and  the  species  yields)  are  constant  over  the 
course  of  the  simulation  and  are  material  specific.  Considerable  care  is  necessary  when  entering 
the  yield  data  to  prevent  violations  of  species  conservation.  Of  particular  note  is  the  situation 
when  FSSIM  is  used  to  calculate  the  oxygen  yield  and/or  the  carbon  dioxide  yield.  In  both  cases, 
it  is  possible  for  that  unrealistic  (negative)  yields  are  calculated  if  the  yield  of  water,  carbon 
monoxide,  and  soot  are  over-specified.  Sources  for  estimating  the  value  of  the  yield  parameters 
include  Tewarson  [2002],  Babrauskas  [2003],  and  Sardqvist  [1993]. 

The  fuel  may  ignite  via  any  of  four  criteria:  a  non-overhead  surface  temperature,  an 
overhead  surface  temperature,  the  compartment  gas  temperature,  or  the  temperature  of 
combustion/gas  products  entering  a  compartment.  If  one  or  more  of  these  criteria  are  not 
implemented,  then  the  default  values  should  be  adjusted  accordingly.  Note  that  the  heat  release 
rate  of  the  fuel  is  proportional  to  t2.  The  maximum  heat  release  rate  is  a  function  of  the  available 
oxygen  in  the  space  because  there  is  no  specific  fuel  package  burning  (i.e.,  a  flashover-like 
condition  is  assumed). 

An  example  of  a  &fuel  NAMELIST  parameter  group  line  is  as  follows: 

&FUEL  ID=11,  growth=0. 00293,  deltah_f=14100000 . 0,  yh20=1.0,  yco=0.03, 
ysoot=0.07,  fileread=. FALSE. ,  ignition_ST=436. 8,  ignition_0V=5000. 0, 
ignition_GT=773.0,  ignition_VT=1273. 0,  fuelname=' Class  A  Combustibles'  / 

The  above  &fuel  line  is  read  from  the  primary  input  file  and  is  identified  with  the  id 
parameter  set  to  1 1  and  the  text  string  ‘class  a  Combustibles’.  The  fire  has  a  growth  rate 
constant  equal  to  0.00293  kW/s2,  which  is  nearly  equal  to  a  ‘slow’ t2  fire.  The  heat  of  combustion 
for  the  materials  is  14  x  106  J/kg,  typical  of  wood-based  products.  The  fuel  has  a  carbon 
monoxide  yield  of  0,03  kg/kg,  a  soot  yield  of  0.07  kg/kg,  and  a  water  yield  of  1.0  kg/kg.  The 
carbon  dioxide  and  oxygen  yields  are  determined  by  FSSIM.  Based  on  the  equations  cited  above, 
the  yields  are  about  2.96  kg  CC>2/kg  fuel  and  3.08  kg  02/kg  fuel,  respectively.  The  ignition 
criteria  is  either  a  boundary  surface  temperature  of 436.8  K  (the  ASTM  El  19-98  temperature 
criteria  for  a  fire  barrier  wall  [ASTM  El  19-98,  1999]),  a  compartment  gas  temperature  of  773.0 
K  (flashover),  or  an  inflowing  gas  temperature  of  1273.0  K  (flames).  The  ignition  criteria  for  an 
overhead  surface  is  set  to  5000.0  K  to  disable  the  feature. 
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4.6  Detection  and  Fire  Suppression 


Detectors  and  fire  suppression  actions  may  be  included  in  a  FSSIM  simulation.  Detection 
may  be  used  as  a  means  of  controlling  fans,  dampers,  and  fire  doors.  The  detection  time  for  a 
particular  type  of  scenario  may  also  be  of  interest  for  reasons  beyond  the  scope  of  this  model, 
such  as  evacuation  and  response  time  estimates.  Suppression  is  a  means  of  directly  interfering 
with  the  combustion  process  within  a  space.  Several  types  of  extinguishing  systems  are  available 
in  FSSIM,  including  water-based  systems  or  systems  using  gaseous  agents. 


4.6.1  Detection  Systems  (sdtct) 

The  sdtct  parameter  group  is  used  to  specify  a  detector  (location  and  type).  This 
parameter  group  is  optional;  however,  if  it  is  used  there  are  two  (out  of  three)  parameters  that 
must  be  set: 


REQUIRED  PARAMETERS  WHEN  A  &DTCT  NAMELIST  GROUP  IS  USED 


id 

(Long  integer)  Identification  number  for  an  individual  sdtct  parameter 
group  line.  The  id  must  be  unique  among  other  sdtct  parameter  id 
values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 
is  no  default  value. 

comp 

(Long  integer)  A  compartment  identification  number  (scomp  id 
parameter  value)  where  the  detector  is  located.  -1  indicates  that  the 
detector  is  not  in  a  compartment.  This  parameter  must  be  -1  if  the  hvac 
parameter  is  greater  than  0.  The  default  value  is  -1 . 

hvac 

(Long  integer)  A  duct  identification  number  (srdct  id  parameter  value) 
where  the  detector  is  located.  -1  indicates  that  the  detector  is  not  in  a 
duct.  This  parameter  must  be  -1  if  the  comp  parameter  is  greater  than  0. 
The  default  value  is  -1 . 

There  is  no  requirement  for  the  sdtct  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  sdtct  id 
parameters  is  acceptable.  The  comp  and  hvac  parameters  are  mutually  exclusive  and  a  fatal  error 
will  occur  if  both  are  specified.  The  detector  may  be  located  in  either  a  compartment  or  a  duct, 
but  not  both.  Note  that  there  is  no  restriction  on  the  number  of  detectors  within  a  compartment  or 
a  duct.  Multiple  sdtct  parameter  groups  may  reference  the  same  compartment  or  duct;  however, 
if  the  detectors  are  identical  they  will  actuate  (or  fail  to  actuate)  at  the  same  time. 

The  remaining  sdtct  parameters  are  optional: 


57 


OPTIONAL  PARAMETERS 


detect  or  name 

(Text  string)  Name  of  the  detector.  The  detector  name  is  used  primarily 
for  labeling  various  output  data.  The  maximum  number  of  characters  for 
the  detector  name  is  40.  The  default  is  ‘detector  id’,  where  id  is  the 
detector  identification  number. 

fileread 

(Logical)  .  true  .  indicates  that  the  detector  parameters  should  be  read 
from  the  file  set  by  the  &  fnam  detectorf  lie  parameter  and  .  false  , 
indicates  that  the  detector  parameters  are  contained  within  the  &dtct  id 
line  located  in  the  main  input  file.  The  default  state  is  .  true  .  If  the 
detectorname  parameter  is  not  set  (i.e.,  the  default  is  used),  then  the 
detector  data  will  not  be  read  from  the  file  set  by  the  &  fnam  detector 
parameter. 

dtype 

(Short  integer)  The  type  of  detection  system: 

1  -  ionization  smoke  detector 

2  -  photoelectric  smoke  detector 

4  -  infrared  heat  detector 

5  -  ultraviolet  electromagnetic  radiation  detector 

6  -  rate  of  rise/thermal  element  heat  detector 

8  -  fixed  temperature  detector 

9  -  carbon  monoxide  detector 

1 0  -  carbon  dioxide  detector 

1 1  -  oxygen  sensor 

The  default  value  is  1. 

operational 

(Logical)  .  true,  indicates  that  the  detector  is  functioning  normally  and 

.  false  .  indicates  that  the  detector  is  not  functioning.  The  default  state  Is 
.TRUE. 

rti 

(Real  number)  The  Response  Time  Index  (RTI)  for  the  detector  device 
(m°'5-s0,5).  This  parameter  is  only  used  by  detectors  that  are  thermally 
actuated  (dtype  set  to  6).  The  default  value  is  10  m°'5-so  s). 

temperature 

(Real  number)  Actuation  temperature  (K)  for  detectors  that  are  thermally 
actuated  (dtype  set  to  6  or  8).  The  default  value  is  0.0  K. 

det_spaeing 

(Real  number)  Average  spacing  between  detectors  (m).  If  this  parameter 
is  not  set  or  a  negative  number  is  provided,  then  the  detector  spacing  will 
be  equal  to: 

where  L  is  the  detector  spacing  (m),  x  is  the  dimension  of  the 
compartment  containing  the  detector  specified  by  the  &comp  dx 
parameter,  andy  is  the  compartment  dimension  set  by  the  &comp  dy 
parameter.  This  parameter  is  not  used  if  the  detector  is  located  in  a  duct. 
The  default  value  is  (scomp  dx  *  scomp  dy)1/2. 
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min  fire 

(Real  number)  The  minimum  detectable  fire  size  (kW)  for  infrared  or 
ultraviolet  detectors  (dtype  set  to  4  or  5).  The  default  value  is  1.0  kW. 

visibitily 

(Real  number)  The  obscuration  density  (m'1)  at  which  actuation  occurs 
for  ionization  or  photoelectric  detectors  (dtype  set  to  1  or  2).  The  default 
value  is  0.0581  m'1. 

massfrac 

(Real  number)  The  mass  fraction  of  carbon  monoxide,  carbon  dioxide,  or 
oxygen  at  which  detection  occurs  for  a  carbon  monoxide,  carbon  dioxide, 
or  oxygen  detector  (dtype  set  to  9,  10,  or  1 1,  respectively).  The  default 
value  is  0.01  kg  CO,  C02,  or  02  per  kg  air. 

The  sdtct  variables  may  be  read  from  a  separate  input  file  as  specified  by  the  &fnam 
detectorf  ile  parameter.  The  method  is  identical  to  that  described  for  the  &cmpn  NAMELIST 
parameter  group  in  Section  4.3.4  and  is  not  repeated  here. 

Detector  actuation  information  is  specific  to  the  type  of  detector  installed.  Information 
regarding  detector  set  points  may  be  available  from  the  manufacturer  of  the  particular  device. 
Other  source  of  information  include  NFPA  72  [2002],  Geiman  and  Gottuk  [2001],  Schifiliti, 
Meachum,  and  Custer  [2002]  and  Beyler  and  DiNenno  [1991]. 

An  example  application  for  the  sdtct  NAMELIST  parameter  group  is  as  follows: 

&DTCT  ID=4,  comp=6,  fileread=. FALSE. ,  dtype=2,  visibility=0 . 18, 
detectorname=' Photoelectric  smoke  detector'  / 

&DTCT  ID=5,  comp=6,  fileread=. FALSE. ,  dtype=8,  temperature=343. 15, 
detectorname=' Fixed  temperature  detector  -  70C  '  / 


The  above  sdtct  lines  are  read  from  the  primary  input  file.  One  is  identified  with  the  id 
parameter  set  to  4  and  the  text  string  ‘photoelectric  smoke  detector’  and  the  other  is 
identified  with  the  id  parameter  set  to  5  and  the  text  string ‘Fixed  temperature  detector  - 
70c’.  Both  detectors  are  located  in  the  compartment  with  the  scomp  id  parameter  set  to  6  and  are 
operational.  The  first  detector  is  a  photoelectric  smoke  detector  that  will  actuate  when  the  smoke 
obscuration  density  is  greater  than  or  equal  to  0.18  m'1.  The  second  detector  is  a  fixed 
temperature  heat  detector  that  will  actuate  when  the  ambient  room  temperature  is  70°C 
(343.15  K)  or  greater. 

4.6.2  Suppression  Systems  (ssupr) 

The  &supr  parameter  group  is  used  to  specify  a  suppression  system  (location  and  type). 
This  parameter  group  is  optional;  however,  if  it  is  used  there  are  four  (out  of  six)  parameters  that 
must  be  set: 
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REQUIRED  PARAMETERS  WHEN  A  &SUPR  NAMELIST  GROUP  IS  USED 


id 

(Long  integer)  Identification  number  for  an  individual  &supr  parameter 
group  line.  The  id  must  be  unique  among  other  &supr  parameter  id 
values  used  in  a  particular  FSSIM  simulation  and  greater  than  zero.  There 
is  no  default  value. 

comp 

(Long  integer)  A  compartment  identification  number  (&comp  id 
parameter  value)  where  the  detector  is  located.  —1  indicates  that  the 
suppression  is  not  in  a  compartment.  This  parameter  must  be  -1  if  the 
hvac  parameter  is  greater  than  0.  Because  HVAC  suppression  is  not 
currently  implemented,  the  comp  parameter  must  be  positive.  The  default 
value  is  -1. 

hvac 

(Long  integer)  A  duct  identification  number  (srdct  id  parameter  value) 
where  the  suppression  is  located.  -1  indicates  that  the  suppression  is  not 
in  a  duct.  This  parameter  must  be  -1  if  the  comp  parameter  is  greater  than 

0.  Because  HVAC  suppression  is  not  currently  implemented,  this 
parameter  must  be  set  to  -1.  The  default  value  is  -1. 

diam 

(Real  number)  The  average  drop  diameter  (m)  for  a  mist  or  sprinkler 
suppression  system  (stype  set  to  1  or  2).  This  is  a  required  parameter 
because  a  mist  suppression  system  is  the  default  and  there  is  no  default 
value  for  the  droplet  diameter.  If  this  parameter  is  not  set,  an  error  would 
thus  be  generated. 

flowrate 

(Real  number  array)  An  array  or  real  numbers,  separated  by  a  comma, 
containing  the  mass  flow  rate  (kg/s)  of  the  suppressing  agents.  The  mass 
flow  rate  is  total  agent  flow  rate  into  a  particular  control  volume 
(compartment  or  duct)  from  all  discharge  devices  associated  with  the 
particular  suppression  system.  The  number  of  array  elements  is  1  if 
stype  is  set  to  1, 2, 4,  5, 6,  or  7.  If  stype  is  set  to  3,  the  number  of 
elements  is  equal  to  the  ngasses  parameter  and  the  elements  correspond 
to  the  individual  gas  components  as  specified  in  the  gas  parameter  array. 
The  order  of  the  array  elements  is  thus  significant.  There  is/are  no  default 
value(s). 

control 

(Long  integer)  Control  function  sctrl  id  parameter  value  used  to 
establish  actuation  criteria  for  suppression  system.  There  is  no  default 
value. 

There  is  no  requirement  for  the  &supr  id  parameter  to  be  sequential  or  in  ascending  or 
descending  order.  Any  valid  integer  greater  than  zero  and  unique  among  other  ssupr  id 


parameters  is  acceptable.  The  comp  and  hvac  parameters  are  mutually  exclusive.  The  suppression 
system  may  be  located  in  either  a  compartment  or  a  duct,  but  not  both.  However,  because  the 
HVAC  suppression  is  not  currently  implemented,  the  suppression  system  must  be  located  in  a 
compartment.  Note  that  there  is  no  restriction  on  the  number  of  suppression  systems  within  a 
compartment  (or  a  duct  when  implemented).  Multiple  &supr  parameter  groups  may  reference  the 
same  compartment  (or  duct). 
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The  remaining  &supr  parameters  are  optional: 


OPTIONAL  PARAMETERS 


suppressorname 

(Text  string)  Name  of  the  suppression  system.  The  suppression  system 
name  is  used  primarily  for  labeling  various  output  data.  The  maximum 
number  of  characters  for  the  suppression  system  name  is  40.  The  default 
is  ‘suppressor  id’,  where  id  is  the  suppression  system  identification 
number. 

fileread 

(Logical)  .  true  .  indicates  that  the  suppression  system  parameters  should 
be  read  from  the  file  set  by  the  &fnam  suppressorf  ile  parameter  and 
.  false  .  indicates  that  the  suppression  system  parameters  are  contained 
within  the  ssupr  id  line  located  in  the  main  input  file.  The  default  state 
is  .  true  .  If  the  suppressorname  parameter  is  not  set  ( i.e .,  the  default  is 
used),  then  the  suppressor  data  will  not  be  read  from  the  file  set  by  the 
sfnam  suppressorfile  parameter. 

stype 

(Short  integer)  The  type  of  suppression  system: 

1  -  mist  suppression  system 

2  -  sprinkler  system 

3  -  gaseous  agent  suppression  system 

4  -  hydrofluorocarbon  (Halon)  suppression  system 

5  -  handline  (not  implemented) 

6  -  foam 

7  -  dry  chemical 

8  -  boundary  cooling  for  a  single  surface 

9  -  boundary  cooling  for  a  compartment  (excludes  the  overhead) 

The  default  value  is  1. 

efficiency 

(Real  number  from  0  to  1)  Denotes  the  evaporation  efficiency  for  a 
water  mist  system.  This  parameter  will  reduce  the  mass  evaporation  by 
multiplying  the  calculated  evaporation  by  the  efficiency.  The  default 
value  is  1. 

operational 

(Logical)  .  true  .  indicates  that  the  suppression  system  is  functioning 
normally  and  .  false  .  indicates  that  the  suppression  system  is  not 
functioning.  The  default  state  is  .  true  . 

tf  low 

(Real  number)  The  discharge  temperature  of  the  suppression  system 
agent  (K).  Parameter  is  used  only  if  stype  is  set  to  1, 2, 3, 4, 5, 6,  8,  or  9. 
The  default  value  is  the  &exec  tamb  parameter  value. 

massfrac 

(Real  number)  The  required  mass  fraction  of  suppression  agent  for 
extinguishment  (kg  agent/kg  air).  This  parameter  is  used  only  if  stype  is 
set  to  4  or  7.  The  default  value  is  1  x  10'6  kg  agent/kg  air. 
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agent  cp 

(Real  number)  The  specific  heat  of  the  suppression  system  agent  (J/kg- 
K).  This  parameter  is  required  if  stype  is  set  to  4;  otherwise,  it  is  not 
used.  There  is  no  default  value. 

discharge  t 

(Real  number)  The  discharge  time  for  the  suppression  system  (s).  The 
default  value  is  1  x  109  s. 

ngasses 

(Short  integer)  The  number  of  gasses  composing  the  gaseous  suppression 
system  agent.  This  parameter  is  required  if  stype  is  set  to  3;  otherwise,  it 
is  not  used.  The  maximum  number  of  gasses  is  6.  There  is  no  default 
value. 

gas 

(Short  integer  array)  An  array  of  short  integers,  from  the  following  list, 
separated  by  commas,  specifying  the  gas  components  of  the  suppression 
system.  The  number  of  elements  is  equal  to  the  ngasses  parameter  and 
thus  may  not  exceed  6.  This  parameter  is  required  if  stype  is  set  to  3; 
otherwise,  it  is  not  used.  The  gas  components  are  identified  as  follows: 

1  -  nitrogen  (N2) 

2- oxygen  (02) 

3  -  fuel 

4  -  carbon  dioxide  (CO2) 

5  —  water  (H2O) 

6  -  carbon  monoxide  (CO) 

The  order  of  the  array  elements  is  not  significant  relative  to  which  gases 
are  specified;  however,  the  gas  flow  rates  specified  using  the  flowrate 
parameter  must  correspond  to  the  gases.  There  are  no  default  values  for 
any  of  the  array  elements. 

surf 

(Long  integer)  A  surface  identification  number  (&surf  id  parameter 
value)  where  the  suppression  is  targeted.  -1  indicates  that  the  suppression 
system  is  not  targeted  at  a  surface  and  is  thus  not  a  boundary  cooling 
device.  This  parameter  is  required  if  stype  is  set  to  8;  otherwise,  it  is  not 
used.  The  default  value  is  -1. 

expansion 

(Long  integer)  The  expansion  ratio  of  the  foam  used  to  suppress  a  fire. 

This  parameter  is  required  if  stype  is  set  to  6;  otherwise,  it  is  not  used. 

The  default  value  is  10.0. 

thickness 

(Long  integer)  The  thickness  of  the  foam  layer  used  to  suppress  a  fire 
(m).  This  parameter  is  required  if  stype  is  set  to  6;  otherwise,  it  is  not 
used.  The  default  value  is  0.1  m. 

The  ssupr  variables  may  be  read  from  a  separate  input  file  as  specified  by  the  sfnam 
suppressorf  ile  parameter.  The  method  is  identical  to  that  described  for  the  scmpn  NAMELIST 
parameter  group  in  Section  4.3.4  and  is  not  repeated  here. 

Suppression  system  parameters  are  specific  to  the  devices  used.  Information  regarding 
suppression  systems  maybe  obtained  from  the  manufacturers.  Information  is  also  available  in 
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various  reference  material  including  Flemming  [2002],  Scheffey  and  Hickey  [2002],  Grant 
[2002],  DiNenno  [2002],  Mawhinney  and  Back  [2002],  NFPA  1 1  [2002],  NFPA  12  [2000], 
NFPA  13  [2002],  NFPA  15  [2001],  NFPA  16  [2003],  and  NFPA  17  [2002]. 

An  example  application  for  the  ssupr  NAMELIST  parameter  group  and  the  associated 
sctrl  lines  is  as  follows: 

&SUPR  ID=4,  comp=6,  fileread=. FALSE. ,  stype=2,  flowrate=8 . 8E-4 ,  diam=0 . 0007, 
control=3,  suppressorname=' Ordinary  sprinkler  system'  / 

&CTRL  id=3,  ctype=l,  trip=.TRUE.,  initial* . FALSE . ,  inputvar=2, 
initial*. FALSE. ,  setpoint=338 . 15  / 

&SUPR  ID=5,  comp=6,  fileread*. FALSE. ,  stype=3,  flowrate=l . 6E-7,  4.3E-6, 
3.9E-6,  control=4 ,  ngasses=3,  gas=2,  1,  4,  discharge_t=35 . 0, 
suppressorname*' Gaseous  agent'  / 

&CTRL  id=4 ,  ctype=l,  trip*. TRUE.,  initial*. FALSE . ,  inputvar=10, 
initial* . FALSE. ,  location=2  / 

The  above  &supr  lines  are  read  from  the  primary  input  file.  One  is  identified  with  the  id 
parameter  set  to  4  and  the  text  string  ‘ordinary  Sprinkler  System’  and  the  other  is  identified 
with  the  id  parameter  set  to  5  and  the  text  string  ‘Gaseous  agent’.  Both  suppression  systems 
are  located  in  the  compartment  with  the  scomp  id  parameter  set  to  6  and  are  operational  but 
initially  not  actuated.  The  first  suppression  system  is  an  ordinary  sprinkler  system  that  actuates 
when  the  compartment  temperature  reaches  65°C  (338.15  K),  has  an  average  drop  diameter  of 
0.7  mm,  and  a  flow  rate  of  8.8  x  10-4  kg/s  (about  15  gpm).  The  second  suppression  system  is  a 
gaseous  agent  consisting  of  nitrogen,  oxygen,  and  carbon  dioxide.  The  mass  flow  rates  for  each 
constituent  are  1.6  x  10'7  kg  02/s,  4.6  x  10'6  kg  N2/s,  and  3.9  x  10'6  kg  C02/s.  The  gaseous 
suppression  system  actuates  by  detection  and  discharges  in  35  seconds.  The  detector 
characteristics  are  specified  by  a  &dtct  parameter  group  with  an  id  set  to  2. 


5.0  FSSIM  OUTPUT 

The  results  of  a  FSSIM  execution  are  interpreted  using  the  output  data.  Currently,  there 
are  two  output  methods  available:  one  that  generates  several  Comma  Separated  Value  (CSV) 
ASCH  files  that  may  be  used  by  a  graphing  or  spreadsheet  program  of  the  users  choice  (general 
output  format)  and  one  that  generates  a  single  Comma  Separated  Value  (CSV)  ASCII  file  that  is 
used  by  the  MSU  viewer  (MSU  viewer  output  format)  [Haupt  et  al,  2004].  As  noted  in  Section 
4.2.1,  the  particular  output  format  is  set  using  the  sexec  out  type  parameter. 

5.1  General  Output  Format 

When  the  sexec  outtype  parameter  is  set  to  .  false  . ,  FSSIM  will  generate  multiple 
CSV  output  files,  each  containing  data  for  related  variables.  The  number  of  files  that  are 
generated  will  depend  on  the  configuration  that  is  modeled.  The  output  files  are  generated  by 
appending  the  output  file  extension  to  the  f  ileid  parameter  set  by  the  &exec  group  (Section 
4.2.1)  along  with  an  integer  indicating  the  sequence  of  the  file  if  more  than  one  is  required  for 
the  output  variable.  For  example  if  there  were  256  compartments,  two  output  files  would  be 
required  and  they  would  be  named  [fiieid]_ctempsi .  csv  and  [f  iieid]_ctemps2 .  csv.  The 
specific  output  files  are  summarized  in  Table  5-1  and  are  described  below. 
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Table  5-1.  Summary  of  FSSIM  Output  Files  -  General  Output  Format 


Output  File 

Output  Data  Variables 

[f  ileid]__ctemps# .  csv 

Compartment  temperatures 

[f  ileid]_ntemps# .  cs? 

Duct  node  temperatures 

[f  ileid]_cpress#  *  csv 

Compartment  pressures 

[f  ileid]_npress#  *  cs? 

Duct  node  pressures 

[fileid]  chrri.esv 

Compartment  heat  release  rate 

[f  Held]  compq#  .  cs? 

Compartment  energy  balance  variables 

[fileid]  jvelst.esv 

Opening  flow  velocities 

[fileid]  dvelsf.csv 

Duct  flow  velocities 

[fileid]  j flows#. csv 

Opening  mass  flow  rates 

[f  ileid]_df  lows# .  csv 

Duct  mass  flow  rates 

[fileid]  jar ea#.cs? 

Opening  flow  areas 

[fileid]  darea# . csv 

Duct  flow  areas 

[fileid]  cgasi.esv 

Mass  fraction  of  gas  components  in  compartments 

[fileid]  ngasi.csv 

Mass  fraction  of  gas  components  in  the  ductwork 

[f ileid]_f  ire#  .  csv 

Fire  theoretical  maximum  heat  release  rate  (i.e.  non-02  limited) 

[f  ileid]__stemp# .  csv 

Surface  temperatures 

[f  ileid]_sf  lux# ,  csv 

Surface  heat  fluxes 

[fileid]  detect#. csv 

Actuation  status  of  detectors 

[fileid]  suppress#.csv 

Actuation  status  of  suppression  system 

Each  output  file  will  have  a  header  consisting  a  comma  separated  line  of  character 
strings.  The  format  for  all  headers  is: 


Time  (s),  item  name(1),  item  name(2)...,  item  name  (nitems) 


where  item  name  is  the  character  string  associated  with  the  output  quantity  (e.g.  for  a 
compartment  quantity  it  will  be  the  compartment  name,  for  a  surface  quantity  the  surface  name) 
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1.  Compartment  temperature  output  data  -  [fiieidl  ctemps# .  csv 


This  file  contains  the  compartment  temperature  data.  There  is  a  description  header 
followed  by  rows  of  output  data  in  the  following  format: 


time,  temp(cl),  temp(c2), ....  temp(ncomp) 


time  The  time  (s)  that  the  temperature  data  correspond  to. 

temp(cn)  The  gas  temperature  (K)  in  the  compartment  set  by  the  cnth  &comp  line, 

ncomp  The  number  of  &comp  lines  as  specified  on  the  &cntr  line. 

If  there  are  more  than  254  compartments  (254  temperatures  plus  the  time  gives  255 
columns),  than  additional  files  will  be  required  for  each  additional  block  of  254  compartments. 

2.  Duct  node  temperature  output  data  -  [f  iieid]_ntemps# .  csv 

This  file  contains  the  compartment  temperature  data.  There  is  a  description  header 
followed  by  rows  of  output  data  in  the  following  format: 


time,  temp(nl),  temp(n2), ....  temp(nnode) 


The  time  (s)  that  the  temperature  data  correspond  to. 

The  gas  temperature  (K)  at  the  duct  node  set  by  the  nnth  &rnod  line. 

The  number  of  &rnod  lines  as  specified  on  the  &cntr  line. 

If  there  are  no  duct  nodes,  then  this  file  will  not  be  created. 

3.  Compartment  pressure  output  data  -  (fiieid]_cpres# .  csv 

This  file  contains  the  compartment  pressure  data.  There  is  a  description  header  followed 
by  rows  of  output  data  in  the  following  format: 


time 

temp(nj) 

nnode 
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time,  pres(cl),  pres(c2),  pres(ncomp) 


The  time  (s)  that  the  pressure  data  correspond  to. 

The  pressure  (Pa)  in  the  compartment  set  by  the  cn,h  &comp  line. 

The  number  of  &comp  lines  as  specified  on  the  &cntr  line. 

4.  Duct  node  pressure  output  data  -  [f  iieid]_npres# .  esv 

This  file  contains  the  duct  node  pressure  data.  There  is  a  description  header  followed  by 
rows  of  output  data  in  the  following  format: 


time,  pres(nl),  pres(n2), ...,  pres(nnode) 


The  time  (s)  that  the  pressure  data  correspond  to. 

The  pressure  (Pa)  at  the  duct  node  set  by  the  nn*h  &rnod  line. 

The  number  of  &rnod  lines  as  specified  on  the  scntr  line. 

If  there  are  no  nodes,  then  this  file  will  be  created. 

5.  Compartment  heat  release  data  -  ff  ileid]_chrr# .  csv 

This  file  contains  the  compartment  heat  release  rate  data.  There  is  a  description  header 
followed  by  rows  of  output  data  in  the  following  format: 


time,  hrr(cl),  hrr(c2), ....  hrr(ncomps) 


The  time  (s)  that  the  heat  release  data  correspond  to. 

The  heat  release  rate  (kW)  at  the  compartment  set  by  the  en,h  &comp  line. 
The  number  of  scomp  lines  as  specified  on  the  &cntr  line. 


time 
hrr(cn) 
n comps 


time 

pres(nn) 

nnode 


time 

pres(cn) 

ncomp 
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If  there  are  no  nodes,  then  this  file  will  be  created. 


6.  Compartment  energy  balance  output  data  -  [f  ileid]_compq# .  csv 

This  file  contains  data  relating  to  the  energy  balance  in  a  compartment.  A  positive  value 
indicates  that  there  is  a  net  increase  in  the  energy  content  of  the  component.  There  is  a 
description  header  followed  by  rows  of  output  data  in  the  following  format: 


time,  qg(cl),  qw(cl),  qg(c2),  qw(c2), ....  qg(ncomp),  qw(ncomp) 


time  The  time  (s)  that  the  energy  balance  data  corresponds  to. 

qg(cn)  The  net  radiation  absorbed  by  the  gas  (W)  in  the  compartment  set  by  the 

cnth  &comp  line. 

qw(cn)  The  net  convection  heat  flux  (W)  to  the  boundaries  of  the  compartment 

set  by  the  cnth  scomp  line. 

ncomp  The  number  of  &comp  lines  as  specified  on  the  &cntr  line. 

Since  two  variables  are  being  written  per  compartment,  additional  files  are  required  if 
there  are  more  than  127  compartments. 

7.  Opening  flow  velocity  output  data  -  [ fileid|_i  vels# .  csv 


This  file  contains  flow  velocity  data  for  openings  specified  using  &  junc  NAMELIST 
group.  Positive  flow  for  &  junc  openings  is  defined  as  flow  from  the  first  compartment  specified 
by  the  first  element  of  &  junc  location  array  to  the  second  compartment  (or  ambient)  specified 
by  the  second  element. 

Bidirectional  &junc  type  openings  are  split  into  two  openings  in  FSSIM  and  the  output 
data  are  reported  accordingly.  The  lower  portion  of  a  vertical  opening  with  horizontal  flow 
retains  the  input  sequence  number.  The  upper  portion  of  a  vertical  opening  is  given  the  next 
available  sequence  number  that  is  greater  than  the  total  number  of  openings  as  set  by  the  &exec 
njunc  parameters.  Thus,  if  three  vertical  openings  were  defined  and  ordered  as  follows: 

&junc  id=6,  bidirectional=.  TRUE . ,  juncname-'  door  A'  ...  / 

&junc  id=3,  bidirectional^ .  FALSE . ,  ...  / 

&junc  id-12,  bidirectional . TRUE . ,  ...  / 

and  the  sexec  njunc  parameter  were  set  to  3  and  the  &exec  nnotbidir  were  set  to  1,  then 
FSSIM  would  create  five  openings  with  sequential  identification  numbers  as  follows: 
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Opening  1  -  lower  portion  of  opening  with  &  junc  id  set  to  6 
Opening  2  -  lower  portion  of  opening  with  &  junc  id  set  to  3 
Opening  3  -  lower  portion  of  opening  with  &  junc  id  set  to  12 
Opening  4  -  upper  portion  of  opening  with  &  junc  id  set  to  6 
Opening  5  -  upper  portion  of  opening  with  &  junc  id  set  to  12 

The  bidirectional  juncname  is  obtained  by  appending  ‘-2 5  to  the  name.  Thus  the 
juncname  for  opening  4  would  be  '  door  A-2'. 

The  flow  velocity  output  is  reported  sequentially  in  exactly  this  manner.  A  similar 
arrangement  applies  to  horizontal  openings,  except  the  notions  of  lower  and  upper  are  not 
applicable. 

The  output  file  contains  a  description  header  followed  by  rows  of  output  data  in  the 
following  format: 


time,  vf!ow(j1),  vflow(j2), ....  vflow(2*njunc-nnotbidir) 


The  time  (s)  that  the  opening  velocity  flow  data  corresponds  to. 

The  flow  velocity  (m/s)  through  the  jnth  sequential  opening. 

The  number  of  &  junc  lines  as  specified  on  the  &cntr  line. 

The  number  &  junc  openings  that  are  not  bidirectional  as  specified  on  the 
&cntr  line. 

If  there  are  no  openings,  then  this  file  will  not  be  created. 

8.  Duct  flow  velocity  OUtPUt  data  -  [f ileid]_dvels#  .  csv 


time 

vflow(jn) 

njunc 

nnotbidir 


This  file  contains  flow  velocity  in  ducts  specified  using  srdct  NAMELIST  group. 
Positive  flow  in  ducts  is  defined  as  flow  from  the  duct  node  specified  by  the  first  element  in  the 
srdct  node  array  to  the  second  duct  node  specified  by  the  second  element  of  the  array. 

The  output  file  contains  a  description  header  followed  by  rows  of  output  data  in  the 
following  format: 


time,  vflow(dl),  vfiow(d2), ....  vflow(nduct) 
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The  time  (s)  that  the  opening  velocity  flow  data  correspond  to. 

The  flow  velocity  (m/s)  through  the  duct  set  by  the  dnth  srdct  line. 
The  number  of  &rdct  lines  as  specified  on  the  &cntr  line. 

If  there  are  no  ducts,  then  this  file  will  not  be  created. 


time 

vflow(dn) 

nduct 


9.  Opening  mass  flow  output  data  -  [f  ileidj  j  flows# .  csv 


This  file  contains  mass  flow  data  for  openings  specified  using  &  junc  NAMELIST  group. 
Positive  opening  mass  flows  are  defined  in  the  same  manner  as  positive  opening  flow  velocities. 
The  mass  flow  through  openings  is  reported  in  the  same  sequential  manner  as  opening  flow 
velocities. 

The  output  file  contains  a  description  header  (that  indicates  when  a  bidirectional  opening 
is  split)  followed  by  rows  of  output  data  in  the  following  format: 


time,  mflow(j1),  mflow(j2) . mflow(2*njunc-nnotbidir) 


The  time  (s)  that  the  opening  mass  flow  data  correspond  to. 

The  mass  flow  (kg/s)  through  the  jnth  sequential  opening. 

The  number  of  &  junc  lines  as  specified  on  the  scntr  line. 

The  number  &  junc  openings  that  are  not  bidirectional  as  specified  on  the 
scntr  line. 

If  there  are  no  openings,  then  this  file  will  not  be  created. 

10.  Duct  mass  flow  output  data  -  |~f  iieid]_df  lows# .  csv 

This  file  contains  mass  flow  data  for  ducts  srdct  NAMELIST  group  and  for  ducts. 
Positive  duct  mass  flows  are  defined  in  the  same  manner  as  positive  duct  flow  velocities. 

The  output  file  contains  a  description  header  followed  by  rows  of  output  data  in  the 
following  format: 


time,  mflow(dl),  mflow(d2) . mflow(nduct) 


time 

mflow(jn) 

njunc 

nnotbidir 
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time,  mflow(d1),  mfiow(d2),  mflow(nduct) 


time  The  time  (s)  that  the  duct  mass  flow  data  correspond  to. 

mflow(dn)  The  mass  flow  (kg/s)  through  the  duct  set  by  the  dn,h  &rdct  line. 

nduct  The  number  of  &rdct  lines  as  specified  on  the  &cntr  line. 

If  there  are  no  ducts,  then  this  file  will  not  be  created. 

11.  Opening  flow  area  output  data  -  ffileid)_iarea# .  csv 

This  file  contains  opening  flow  area  data.  The  flow  area  of  an  opening  may  vary  with 
time  through  the  use  of  control  functions  set  using  a  sctrl  line  or  for  bidirectional  flows  where 
the  area  of  each  sub-opening  is  determined  as  part  of  the  mass  balance  across  the  opening.  The 
area  for  openings  set  using  the  &  junc  group  is  reported  in  the  same  manner  as  the  flow 
velocities. 

The  output  file  contains  a  description  header  (that  indicates  when  a  bidirectional  opening 
is  split)  followed  by  rows  of  output  data  in  the  following  format: 


time,  Area(j1),  Area(j2), ...,  Area(2*njunc-nnotbidir) 


The  time  (s)  that  the  opening  area  data  correspond  to. 

The  flow  area  (m  )  of  the  jn  sequential  opening. 

The  number  of  &  junc  lines  as  specified  on  the  &cntr  line. 

The  number  &  junc  openings  that  are  not  bidirectional  as  specified  on  the 
&CNTR  line. 


time 

Area(jn) 

njunc 

nnotbidir 


If  there  are  no  openings,  then  this  file  will  not  be  created. 

12.  Duct  flow  area  output  data  -  [f  iieid]_darea# .  esv 

This  file  contains  duct  flow  area  data.  The  duct  flow  area  may  vary  with  time  through  the 
use  of  control  functions  set  using  a  &ctrl  line. 
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The  output  file  contains  a  description  header  followed  by  rows  of  output  data  in  the 
following  format: 


time,  Area(dl),  Area(d2), Area(nduct) 


time  The  time  (s)  that  the  opening  area  data  correspond  to. 

Area(d1 )  The  flow  area  (m2)  of  the  duct  set  by  the  first  srdct  line, 

nduct  The  number  of  &rdct  lines  as  specified  on  the  scntr  line. 

If  there  are  no  ducts,  then  this  file  will  not  be  created. 


13.  Gas  mass  fraction  output  data  (compartments!  -  |~f  iieid]_caas# .  csv 


This  file  contains  the  mass  fraction  of  various  component  gases  in  each  compartment. 
There  is  a  description  header  followed  by  rows  of  output  data  in  the  following  format: 


time,  02(c1).  02(c2) . O2(r:comp),  CO(cl),  CO(c2) . CO(ncomp),  C02(c1),  C02(c2) . C02(ncomp),H20(c1),  H20(c2) . H20(ncomp), 

C(cl),  C(c2) . C(ncomp).  UHC(c1),UHC(c2) . UHC(ncomp) 


time  The  time  (s)  that  the  gas  concentration  data  correspond  to. 

02(cn)  The  oxygen  mass  fraction  (kg  02/kg  gas)  in  the  compartment  set  by  the 

cnth  &  comp  line.  Only  written  if  gaso2=.  TRUE,  on  the  &OUTP  line. 

CO(cn)  The  carbon  monoxide  mass  fraction  (kg  CO/kg  gas)  in  the  compartment 

set  by  the  cnth  scomp  line.  Only  written  if  gas  co= .  TRUE .  on  the 
&OUTP  line. 

C02(cn)  The  carbon  dioxide  mass  fraction  (kg  C02/kg  gas)  in  the  compartment  set 

by  the  cnth  &comp  line.  Only  written  if  gasco= .  TRUE .  on  the  &OUTP 
line. 

H20(cn)  The  water  vapor  mass  fraction  (kg  f^O/kg  gas)  in  the  compartment  set 

by  the  cnth  &comp  line.  Only  written  if  gash2o= .  TRUE .  on  the  &OUTP 
line. 

C(cn)  The  soot  mass  fraction  (kg  soot/kg  gas)  in  the  compartment  set  by  the 

cnth  scomp  line.  Only  written  if  gasc= .  TRUE .  on  the  &OUTP  line. 
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UHC(cn)  The  unbumed  hydrocarbon  mass  fraction  (kg  UHC/kg  gas)  in  the 

compartment  set  by  the  cnth  &comp  line.  Only  written  if 
gasfuel=.TRUE.  on  the  &  OUT  P  line. 

ncomp  The  number  of  &comp  lines  as  specified  on  the  &cntr  line. 


This  file  will  not  be  created  if  all  ga  s  ?  ?  ?  variables  are  set  to  .  FALSE .  on  the  &OUTP 
line.  Additional  files  will  be  required  if  there  are  more  than  254  /  (number  of  output  gases) 
compartments. 

14.  Gas  mass  fraction  output  data  (ducts')  -  ff ileidl  naas .  csv 


This  file  contains  the  mass  fraction  of  various  component  gases  at  each  duct  node.  There 
is  a  description  header  followed  by  rows  of  output  data  in  the  following  format: 


time,  02(n1),  02(n2) . 02{nnode),  CO(nl),  CO(n2) . CO(nnode),  C02(n1),  C02(n2), ....  CO2{nnode),H20(n1),  H20(n2) . H20{nnode), 

C(nl),  C(n2), ....  C(nnode),  UHC(n1).UHC(n2) . UHC(nnode) 


time 

02(nn) 

CO(nn) 


C02(nn) 

H20(nn) 


C(nn) 

UHC(nn) 


ncomp 


The  time  (s)  that  the  gas  concentration  data  correspond  to. 

The  oxygen  mass  fraction  (kg  02/kg  gas)  in  the  duct  node  set  by  the  nnth 
&rnod  line.  Only  written  if  gaso2=.  TRUE,  on  the  &  OUT  P  line. 

The  carbon  monoxide  mass  fraction  (kg  CO/kg  gas)  in  the  duct  node  set 
by  the  nnth  srnod  line.  Only  written  if  gasco= .  TRUE .  on  the  &OUTP 
line. 

The  carbon  dioxide  mass  fraction  (kg  C02/kg  gas)  in  the  duct  node  set  by 
the  nn,b  srnod  line.  Only  written  if  gasco= .  TRUE .  on  the  &OUTP  line. 

The  water  vapor  mass  fraction  (kg  f^O/kg  gas)  in  the  duct  node  set  by 
the  nnth  &rnod  line.  Only  written  if  gash2o= .  TRUE .  on  the  &OUTP 
line. 

The  soot  mass  fraction  (kg  soot/kg  gas)  in  the  duct  node  set  by  the  nn*h 
srnod  line.  Only  written  if  gas  c=.  TRUE,  on  the  &OUTP  line. 

The  unbumed  hydrocarbon  mass  fraction  (kg  UHC/kg  gas)  in  the  duct 
node  set  by  the  nn*  &rnod  line.  Only  written  if  gasfuel=.  TRUE .  on 
the  &OUTP  line. 

The  number  of  scomp  lines  as  specified  on  the  scntr  line. 
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If  there  are  no  ducts,  then  this  file  will  not  be  created.  Also  this  file  will  not  be  created  if 
all  gas???  variables  are  set  to  .FALSE,  on  the  &OUTP  line.  Additional  files  will  be 
required  if  there  are  more  than  254  /  (number  of  output  gases)  duct  nodes. 

15.  Energy  release  output  data  -  [f  iieid]_f  ire# .  csv 


This  file  gives  the  theoretical  maximum  energy  release  data  in  each  compartment.  The 
output  file  contains  a  description  header  followed  by  rows  of  output  data  in  the  following  format: 


time,  FireHRR(cl),  FireHRR(c2), ...,  FireHRR(ncomp) 


time 

FireHRR(cn) 

ncomp 


The  time  (s)  that  the  energy  release  data  correspond  to. 

The  heat  release  rate  of  the  fire  (kW)  as  specified  by  the  &fire  line  or  the 
suses  and  &fuel  lines  for  the  compartment  set  by  the  cnth  scomp  line. 

The  number  of  &comp  lines  as  specified  on  the  scntr  line. 


16.  Surface  temperature  output  data  -  [fiieid]_stemp# .  csv 


This  file  contains  the  surface  temperature  data  for  surfaces  set  using  a  &surf  line.  The 
surface  front  is  defined  as  the  face  of  the  surface  that  is  located  in  the  compartment  specified 
with  the  first  element  in  the  &surf  location  parameter  array.  The  surface  back  is  the  face  of  the 
surface  that  is  located  in  the  compartment  specified  by  the  second  element  in  the  &surf 
location  parameter  array. 

The  output  file  contains  a  description  header  followed  by  rows  of  output  data  in  the 
following  format: 


time,  stempF(sl),  stempB(sl),  stempF(s2),  stempB(s2), ....  stempF(nsurf),  stempB(nsurf) 


time 

stempF(sn) 
stempB  (sn) 
nsurf 


The  time  (s)  that  the  surface  temperature  data  correspond  to. 

The  front  temperature  (K)  of  the  surface  set  by  the  snth  &surf  line. 
The  back  temperature  (K)  of  the  surface  set  by  the  snth  &surf  line. 
The  number  of  &surf  lines  as  specified  on  the  scntr  line. 
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If  there  are  no  surfaces,  then  this  file  will  not  be  created.  Additional  output  files  will  be 
required  if  there  are  more  than  127  surfaces. 


17.  Surface  heat  flux  output  data  -  [fiieid]_sfiux# .  csv 


This  file  contains  the  surface  heat  flux  data  for  surfaces  set  using  a  ssurf  line.  The 
surface  front  and  back  are  as  described  for  the  Surface  temperature  data.  Positive  heat  flux  is 
from  the  compartment  to  the  surface.  The  output  file  contains  a  description  header  followed  by 
rows  of  output  data  in  the  following  format: 


time,  scFfsl),  srF(sl),  scB(sl),  srB(sl),  scF(s2),  srF(s2),  scB(s2),  srB(s2),  scF(nsurf),  srF(nsurf),  scB(nsurf),  srB(nsurf) 


time 

scF(sn) 

srF  (sn) 

scB(sn) 

srB  (sn) 

nsurf 


The  time  (s)  that  the  surface  temperature  data  correspond  to. 

The  front  net  convective  heat  flux  (W)  of  the  surface  set  by  the  snth 
ssurf  line. 

The  front  net  radiant  heat  flux  (W)  of  the  surface  set  by  the  snth  &surf 
line. 

The  back  net  convective  heat  flux  (W)  of  the  surface  set  by  the  sn* 
ssurf  line. 

The  back  net  radiant  heat  flux  (W)  of  the  surface  set  by  the  snth  ssurf 
line. 

The  number  of  &surf  lines  as  specified  on  the  scntr  line. 


If  there  are  no  surfaces,  then  this  file  will  not  be  created.  Additional  output  files  will  be 
required  if  there  are  more  than  63  surfaces. 


18.  Detector  actuation  status  data  -  |f iieid]_detect# .  csv 


This  file  contains  the  actuation  status  of  all  detectors  in  the  configuration  modeled.  A 
detector  is  either  actuated  (1  [ .  true  .  ])  or  not  actuated  (0  [ .  false  .  ]).  The  output  file  contains  a 
description  header  followed  by  rows  of  output  data  in  the  following  format: 


time,  status(dl),  status(d2), ....  status(ndetect) 
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time 

status(dl) 

status(d2) 

status(ndetect) 

ndetect 


The  time  (s)  that  the  detector  actuation  status  information  correspond  to. 

The  detector  actuation  status  of  the  detector  set  by  the  first  sdtct  line. 

The  detector  actuation  status  of  the  detector  set  by  the  second  &dtct  line 

The  detector  actuation  status  of  the  detector  set  by  the  ndetectth  &dtct 
line. 

The  number  of  sdtct  lines  as  specified  on  the  scntr  line. 


If  there  are  no  detectors,  then  this  file  will  not  be  created. 


19.  Suppression  system  actuation  status  data  -  ffileidl  suppress# .  csv 


This  file  contains  the  actuation  status  of  all  suppression  systems  in  the  configuration 
modeled.  A  suppression  system  is  either  actuated  ( 1  [ .  true  .  ])  or  not  actuated  (0  [ .  false  .  ]). 
This  is  solely  and  indication  of  the  control  function  that  operates  the  suppression  system  and  not 
whether  or  not  the  suppression  system  has  exceeded  its  discharge  t .  The  output  file  contains 
a  description  header  followed  by  rows  of  output  data  in  the  following  format: 


time,  status(sl),  status(s2),  ....  status(nsuppress) 


time  The  time  (s)  that  the  detector  actuation  status  information  correspond  to. 

status(sn)  The  suppression  system  actuation  status  of  the  suppression  system  set  by 

the  snth  &supr  line. 

nsuppress  The  number  of  &supr  lines  as  specified  on  the  scntr  line. 

If  there  are  no  suppression  systems,  then  this  file  will  not  be  created. 

5.2  MSU  Viewer  Output  Format 

When  the  &exec  out  type  parameter  is  set  to  .  true  . ,  FSSIM  will  generate  a  single  CSV 
output  file  containing  data  in  a  format  that  may  be  read  by  the  MSU  viewer  [Haupt,  et  al.,  2004]. 
The  output  file  is  named  [f  iieid]_out .  csv,  where  f  ileid  is  set  by  the  sexec  parameter  group 
line  (Section  4.2.1). 

The  MSU  output  file  contains  all  output  information  for  all  timesteps  within  a  single 
column  of  data.  The  output  file  consists  of  series  of  headers  followed  by  data.  Headers  typically 
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occupy  two  sequential  rows:  the  first  row  contains  text  that  identifies  the  variable  and  the  second 
row  lists  the  units  that  the  variable  is  reported  in.  Data  between  timestep  headers  correspond  to 
the  first  time  below  the  first  timestep  header.  Data  below  other  headers  corresponds  to  the 
variable  identified  by  the  header  the  data  are  below.  This  pattern  is  repeated  until  the  last 
timestep.  The  data,  in  the  row  order  in  which  it  is  contained,  are  as  follows: 


•  Time  (s)  -  single  row  below  timestep  header. 

•  Compartment  temperature  (K)  -  next  ncomp  rows  below  appropriate  header. 

•  Front  surface  temperature  (K)  -  next  nsurf  rows  below  appropriate  header. 

•  Back  surface  temperature  (K)  -  next  nsurf  rows  below  appropriate  header. 

•  Front  surface  convective  heat  flux  (kW/m2)  -  next  nsurf  rows  below  appropriate 
header. 

•  Back  surface  convective  heat  flux  (kW/m2)  —  next  nsurf  rows  below  appropriate 
header. 

•  Front  surface  radiant  heat  flux  (kW/m2)  -  next  nsurf  rows  below  appropriate 
header. 

•  Back  surface  radiant  heat  flux  (kW/m2)  -  next  nsurf  rows  below  appropriate 
header. 

•  Compartment  pressure  (Pa)  -  next  ncomp  rows  below  appropriate  header. 

•  Duct  node  pressure  (Pa)  -  next  nnode  rows  below  appropriate  header. 

•  Oxygen  mass  fraction  (kg  02/kg  gas)  -  next  ncomp  rows  below  appropriate 
header. 

•  Carbon  monoxide  mass  fraction  (kg  CO/kg  gas)  -  next  ncomp  rows  below 
appropriate  header. 

•  Soot  mass  fraction  (kg  soot/kg  gas)  -  next  ncomp  rows  below  appropriate  header. 

•  Compartment  heat  release  rate  (kW)  -  next  ncomp  rows  below  appropriate  header. 

•  Compartment  fire  size  (kW)  -  next  ncomp  rows  below  appropriate  header. 

•  Suppression  system  actuation  status  -  next  nsuppress  rows  below  appropriate 
header. 


where: 

ncomp  is  the  number  of  scomp  lines  as  specified  on  the  &cntr  line, 
nsurf  is  the  number  of  &surf  lines  as  specified  on  the  scntr  line, 
nnode  is  the  number  of  srnod  lines  as  specified  on  the  &cntr  line, 
nsuppress  is  the  number  of  &supr  lines  as  specified  on  the  scntr  line. 
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6.0  SAMPLE  CASES 


Four  samples  are  provided  to  illustrate  various  modeling  features  of  FSSIM.  The  samples 
are  as  follows: 

•  Simple  configuration  with  a  fire  in  a  space  that  is  connected  to  another  space  via  a 
door.  Each  space  is  vented  to  the  exterior  via  a  window.  Boundaries  are  steel. 

•  A  fire  in  a  single  space  with  an  irregular  floor  plan.  The  boundaries  are  insulated 
steel  and  there  is  a  fixed  water-based  suppression  system. 

•  A  fire  in  a  single  ‘U’-shaped  space  that  must  be  broken  into  three  compartments 
in  the  model.  The  boundaries  are  a  mix  of  steel  and  gypsum. 

•  Five  spaces  connected  by  a  single  corridor.  Three  of  the  spaces  are  mechanically 
vented. 


6.1  Sample  1:  Fire  in  Steel  Compartment  Connected  to  One  Adjacent  Compartment 

This  sample  consists  of  two  spaces  with  steel  boundaries  that  are  connected  via  a  single 
door  that  remains  open.  Each  space  is  vented  to  ambient  via  a  window.  A  fire  is  located  in  one  of 
the  spaces.  Refer  to  Figure  6-1  for  a  depiction  of  Sample  1.  The  following  details  apply: 

•  Each  space  is  5  m  x  5  m  by  3  m  tall.  The  boundaries  are  15.8  mm  thick  carbon 
steel. 

•  The  door  between  the  spaces  is  2. 1 3  m  by  0.9 1  m  wide.  The  door  is  located  in  the 
fore  of  the  fire  compartment  and  the  aft  of  the  adjacent  compartment. 

•  The  window  opening  in  each  space  measures  0.91  m  by  0.61  m  wide  and  is 
located  1.22  m  above  the  deck.  The  window  in  the  fire  compartment  is  located  on 
the  starboard  side  and  the  window  in  the  adjacent  compartment  is  located  on  the 
forward  side. 

•  The  source  fire  is  located  in  Room  1  and  is  a  ‘fast’ t2  fire  with  a  peak  heat  release 
rate  of  1,000  kW.  The  heat  of  combustion  of  the  burning  material  is  28,100  kJ/kg. 
The  default  species  yields  are  applicable  and  the  radiant  fraction  is  0.3. 

This  is  a  simple  example  that  illustrates  how  to  construct  an  input  file  and  set  various 
parameters  for  a  particular  simulation. 
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FORWARD 


Figure  6-1  Sample  1  (Not  to  Scale):  Two  adjacent  steel  compartments 


The  input  file  for  Sample  1  is  contained  in  Appendix  A-l.  The  simulation  is  1,200 
seconds  long.  The  material  properties  for  the  carbon  steel  are  contained  in  the  file 
materials .  txt.  The  portions  of  this  file  used  by  Sample  1  are  also  contained  in  Appendix  A-l, 
Note  that  fake  surfaces  are  applied  to  the  openings  to  allow  for  thermal  radiation  heat  transfer. 
The  calculated  temperature  and  pressure  in  the  two  compartments  are  shown  in  Figure  6-2. 
Figure  6-3  shows  the  concentration  of  various  gas  and  combustion  products  in  the  fire 
compartment  as  well  as  the  heat  release  rate  of  the  fire.  Other  outputs  available,  but  not  shown, 
include  the  boundary  temperature,  the  flux  to  the  boundaries,  and  vent  flow  information. 


- Temperature  -  fire  compartment 

—  —  —  Pressure  -  adjacent  compartment 

Figure  6-2.  Calculated  Compartment  Temperature  and  Pressure  as  a  Function  of  Time 
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Figure  6-3.  Calculated  Mass  Fraction  of  Various  Gases  in  Fire  Compartment  and  Fire  Heat  Release 

Rate 


6.2  Sample  2:  Fire  in  Insulated  Steel  Compartment  with  Irregular  Plan  and  Fixed 
Suppression 

This  sample  consists  of  a  single  space  with  an  irregular  floor  plan.  The  boundaries  are 
steel  with  insulation  backing.  There  is  one  door  and  one  vent  near  the  overhead  that  connects  the 
space  to  the  ambient  surroundings.  Refer  to  Figure  6-4  for  a  depiction  of  Sample  2.  The 
following  details  apply: 

•  The  space  has  four  walls  with  lengths  in  clockwise  order  beginning  with  the 
starboard  wall  of  4.0  m,  3.0  m,  7.0  m,  and  5.0  m.  The  overhead  height  is  3.0  m. 

•  There  is  a  2.13  m  tall  by  0.91  m  wide  door  in  the  aft  wall  that  connects  the  space 
to  ambient  conditions. 

•  There  is  a  0.3  m  by  1 .2  m  wide  vent  in  the  port  wall  located  0.5  m  below  the 
overhead  that  connects  the  space  to  ambient  conditions. 

•  The  walls  and  overhead  are  15.8  mm  thick  carbon  steel  with  25  mm  ceramic 
fiberboard  insulation  for  backing.  The  floor  is  25  mm  thick  carbon  steel  with  38 
mm  Fiberfrax  blanket  insulation  for  backing. 

•  The  source  fire  is  a  ‘medium’ t2  fire  with  a  peak  heat  release  rate  of  800  kW.  The 
heat  of  combustion  of  the  burning  material  is  44,000  kJ/kg.  The  default  species 
yields  are  applicable  and  the  radiant  fraction  is  0.3. 
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•  The  suppression  system  is  based  on  a  single  sprinkler  actuating  in  the  space.  The 
sprinkler  has  an  actuation  temperature  of  74°C  and  a  minimum  water  flow  rate  of 
1.36  kg/s.  The  average  drop  diameter  is  500  microns. 


FORWARD 


Figure  6-4  Sample  2  (Not  to  Scale):  Simple  space  with  irregular  floor  plan. 


Because  FSSIM  simulates  conditions  in  volumes  with  rectangular  floor  plans  and  flat 
decks  and  overheads,  some  modifications  are  necessary  to  simulate  a  fire  in  the  space  shown  in 
Figure  6-4.  The  FSSIM  compartment  definition  should  have  a  floor  plan  with  an  aspect  ratio 
similar  to  the  actual  floor  plan.  The  total  volume  of  the  compartment  defined  in  FSSIM  should 
not  be  less  than  the  actual  compartment  volume.  In  the  case  of  Sample  2,  the  dimension  in  the 
port-starboard  directions  is  3.0  m,  which  is  the  actual  plane  separation  of  the  port  and  starboard 
walls.  Thus,  dx  is  set  to  3.0  on  the  scomp  line.  Similarly,  the  dimension  in  the  fore-aft  separation 
range  from  a  minimum  of  4.0  m  to  a  maximum  of  7.0  m.  An  average  value  is  used  in  the 
simulation,  which  is  done  by  setting  dy  to  5.0  on  the  &comp  line.  In  this  case,  the  actual  volume 
of  the  compartment  equals  the  volume  determined  by  multiplying  dx,  dy,  and  dz.  It  would  also 
be  acceptable  to  set  dy  to  a  value  greater  than  5.0  but  less  than  or  equal  to  7.0  without  changing 
the  volume  (voi)  parameter.  This  would  increase  the  radiant  heat  transfer  between  perpendicular 
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surfaces  and  decrease  it  between  parallel  surfaces  because  the  radiation  configuration  factors 
would  change. 


7.0  m 


Compartment 

Definition 


Figure  6-5  Plan  View  of  Sample  2  (Not  to  Scale) 


The  input  file  for  Sample  2  is  contained  in  Appendix  A-2.  The  simulation  is  1,200 
seconds  long.  The  material  properties  for  the  carbon  steel,  ceramic  fiberboard,  and  Fiberffax 
blanket  insulation  are  contained  in  the  file  ‘materials .  txt\  The  portions  of  this  file  used  by 
Sample  2  are  also  contained  in  Appendix  A-2.  Note  that  fake  surfaces  are  applied  to  the  openings 
to  allow  for  thermal  radiation  heat  transfer.  The  calculated  temperature  and  pressure  in  the  two 
compartments  are  shown  in  Figure  6-6.  Figure  6-7  shows  the  concentration  of  various  gas  and 
combustion  products  in  the  fire  compartment  as  well  as  the  heat  release  rate  of  the  fire.  Other 
outputs  available,  but  not  shown,  include  the  boundary  temperature,  the  flux  to  the  boundaries, 
and  vent  flow  information.  Figure  6-7  indicates  that  the  sprinkler  actuates  and  begins  suppressing 
the  source  fire  around  4  minutes  after  ignition. 
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Figure  6-6.  Calculated  Compartment  Temperature  and  Pressure  as  a  Function  of  Time 
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Figure  6-7.  Calculated  Mass  Fraction  of  Various  Gases  in  Fire  Compartment  and  Fire  Heat  Release  Rate 


6.3  Sample  3:  Fire  in  ‘U ’-Shaped  Compartment  with  Gypsum  and  Steel  Boundaries 

This  sample  consists  of  a  single  space  with  a  ‘U’-shape  floor  plan  where  the  ’U’  opens  to 
the  aft  direction.  The  boundaries  are  a  mix  of  steel  and  gypsum  wallboard.  There  is  one  door  in 
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the  port  section  of  the  space  and  two  vents,  one  on  the  starboard  boundary  and  one  on  the  port 
boundary.  Sample  3  is  shown  in  Figure  6-8.  The  following  details  apply: 

•  The  fire  is  located  in  the  starboard  section  of  the  compartment.  The  default 
species  yields  are  applicable  and  the  radiant  fraction  is  0.3. 

•  The  starboard  and  port  walls  are  1 1  m  long;  the  forward  wall  is  12  m  long;  the 
starboard  and  port  sections  are  3  m  wide  and  the  forward  section  is  4  m  wide. 

•  The  overhead,  aft  and  port  walls  of  the  starboard  section  and  the  aft  and  starboard 
walls  of  the  port  section  are  15.8  mm  thick  gypsum.  The  remaining  surfaces  are 
12.7  mm  thick  carbon  steel. 

•  The  overhead  height  is  3  m  in  all  three  sections. 

•  The  source  fire  has  a  linear  growth  rate  from  0  kW  to  1,200  kW  for  120  seconds 
after  ignition.  The  heat  of  combustion  of  the  burning  material  is  44,000  kJ/kg. 


Because  FSSIM  simulates  conditions  in  volumes  with  rectangular  floor  plans  and  flat 
decks  and  overheads,  some  modifications  are  necessary  to  simulate  a  fire  in  the  space  shown  in 
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Figure  6-8.  The  space  was  divided  into  three  sub-compartments:  the  port  section,  the  forward 
section,  and  the  aft  section.  The  floor  plan  is  shown  in  Figure  6-9.  The  three  sub-compartments 
are  connected  using  an  opening  with  a  cross  sectional  area  equal  to  the  width  of  the  forward 
compartment  times  the  height  of  the  overhead.  A  fake  boundary  was  specified  at  the  connections 
to  allow  for  radiant  heat  transfer.  The  loss  coefficient  was  reduced  to  1.0,  which  corresponds  to  a 
discharge  coefficient  of  1.0. 
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ACTUAL  CONFIGURATION  CONFIGURATION  MODELED 

Figure  6-9  Plan  View  of  Sample  3  Configuration  Modeled  using  FSSIM  (Not  to  Scale) 


The  input  file  for  Sample  3  is  contained  in  Appendix  A-3.  The  simulation  is  1,200 
seconds  long.  The  material  properties  for  the  carbon  steel  and  gypsum  are  contained  in  the  file 
‘materials .  txt\  The  portions  of  this  file  used  by  Sample  3  are  also  contained  in  Appendix  A- 
3.  Note  that  fake  surfaces  are  applied  to  the  openings  to  allow  for  thermal  radiation  heat  transfer. 
The  calculated  temperature  and  pressure  in  the  two  compartments  are  shown  in  Figure  6-10. 
Figure  6-11  shows  the  concentration  of  various  gas  and  combustion  products  in  the  fire 
compartment  as  well  as  the  heat  release  rate  of  the  fire.  Other  outputs  available,  but  not  shown, 
include  the  boundary  temperature,  the  flux  to  the  boundaries,  and  vent  flow  information. 
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-  Temperature  -  starboard  segment  (fire  area) 

.  Temperature  -  forward  segment 

- Temperature  -  starboard  segment 


■  . .  Pressure  *  port  segment  (fire  area) 

. .  Pressure  -  forward  segment 

— — —  Pressure  -  starboard  segment 


Figure  6-10.  Calculated  Compartment  Temperature  and  Pressure  as  a  Function  of  Time 
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Figure  6-11.  Calculated  Mass  Fraction  of  Various  Gases  in  Fire  Segment  and  Fire  Heat  Release  Rate 


6.4  Sample  4:  Multiple  Compartment  Arrangement  Connected  by  a  Corridor 

This  sample  consists  of  five  spaces  connected  by  a  corridor.  The  boundaries  are  15.8  mm 
thick  gypsum.  Each  space  is  connected  to  the  corridor  via  a  single  door.  The  corridor  is 
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connected  to  ambient  at  one  end  via  a  standard  double  door  that  remains  open  for  the  duration  of 
the  scenario.  Three  of  the  spaces  are  connected  to  a  mechanical  ventilation  system.  The  two 
spaces  that  do  not  have  mechanical  ventilation  have  a  small  vent  to  ambient  located  near  the 
ceiling.  A  plan  view  of  Sample  4  is  shown  in  Figure  6-12.  The  following  details  apply: 

•  The  fire  is  in  Room  3  as  shown  in  Figure  6-12.  The  default  species  yields  are 
applicable  and  the  radiant  fraction  is  0.3. 

•  The  plan  dimensions  are  as  shown  in  Figure  6-12.  The  overhead  elevation  is  3  m. 

•  Single  doors  are  0.91  m  wide  by  2.13  m  tall;  double  doors  are  1.82  m  wide  by 
2.13  m  tall. 

•  The  source  fire  has  a  medium  growth  rate  with  a  maximum  heat  release  rate  of 
2,500  kW.  The  heat  of  combustion  of  the  burning  material  is  14,000  kJ/kg. 

•  All  ductwork  associated  with  the  mechanical  ventilation  system  has  a  rectangular 
cross  section  measuring  0.61  m  by  0.61  m.  The  duct  centerline  elevation  is  2.5  m. 
The  duct  lengths  are  as  shown  in  Figure  6-12. 

•  The  duct  roughness  is  0.0005  m. 

•  The  fan  performance  curve  is  quadratic.  The  flow  at  zero  pressure  is  1 .86  m3/s 
and  the  pressure  at  zero  flow  is  180  Pa. 

•  The  pressure  loss  coefficients  associated  with  the  duct ‘T’  connections  are  as 
follows: 

•  Flow  involving  ‘T’  stem  (flow  changes  direction):  1.63 

•  ‘T5  crosses  (no  change  in  flow  direction):  0.75 

•  The  pressure  loss  coefficients  associated  with  an  elbow  is  1.1  in  either  flow 
direction  through  the  elbows. 

•  The  pressure  loss  coefficient  associated  with  inlets  and  outlets  is  0.0. 
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Figure  6-12  Sample  4  (Not  to  Scale) 


As  was  the  case  for  Samples  2  and  3,  some  modifications  to  the  actual  floor  plan  are 
made  to  simulate  a  fire  in  the  space  shown  in  Figure  6-12.  Each  individual  space  was  simulated 
using  the  dimensions  shown  in  Figure  6-12.  The  corridor  was  subdivided  into  three 
approximately  equally  volume  segments  for  the  same  reasons  described  for  Sample  3.  While  this 
is  not  necessary,  it  may  improve  mass  and  energy  transfers  to  the  remote  compartments.  The 
floor  plan  of  the  configuration  simulated  using  FSSIM  is  shown  in  Figure  6-13.  The  three  sub¬ 
compartments  are  connected  using  an  opening  with  a  cross  sectional  area  equal  to  the  width  of 
the  forward  compartment  times  the  height  of  the  overhead.  A  fake  boundary  was  specified  at  the 
connections  to  allow  for  radiant  heat  transfer.  The  loss  coefficient  was  reduced  to  1.0,  which 
corresponds  to  a  discharge  coefficient  of  1.0. 
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Figure  6-13  Plan  View  of  Sample  4  Configuration  Modeled  using  FSSIM  (Not  to  Scale) 


The  input  file  for  Sample  4  is  contained  in  Appendix  A-4.  The  simulation  is  1,200 
seconds  long.  The  material  properties  for  gypsum  are  contained  in  the  file  ‘materials .  txt\ 

The  portions  of  this  file  used  by  Sample  4  are  also  contained  in  Appendix  A-4.  Note  that  fake 
surfaces  are  applied  to  the  openings  to  allow  for  thermal  radiation  heat  transfer.  The  calculated 
temperature  in  the  five  compartments  and  the  three  corridor  segments  is  shown  in  Figure  6-14. 
Likewise,  the  calculated  compartment  pressures  are  shown  in  Figure  6-15.  Figure  6-16  shows  the 
concentration  of  various  gas  and  combustion  products  in  the  fire  compartment  as  well  as  the  heat 
release  rate  of  the  fire.  Other  outputs  available,  but  not  shown,  include  the  boundary  temperature, 
the  flux  to  the  boundaries,  and  vent  flow  information 
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- Temperature  -  room  5 

-  Temperature  -  corridor  segment  1 
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-  Temperature  -  corridor  segment  3 


Figure  6-14.  Calculated  Compartment  Temperatures  as  a  Function  of  Time 
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Figure  6-15.  Calculated  Compartment  Pressures  as  a  Function  of  Time 
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Figure  6-16,  Calculated  Mass  Fraction  of  Various  Gases  in  Fire  Segment  and  Fire  Heat  Release  Rate 
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8.0  NOMENCLATURE 

Roman 

A 
c 
L 
n 
N 
K 
P 
x 

y 

Greek 
Superscript 
Subscript 

CO  Carbon  monoxide 


Area  (m2) 
velocity  coefficient 
Detector  separation  (m) 
number  of  elements  in  an  array 
Number  of  ducts 
Flow  loss  coefficient 
Perimeter  (m) 

Compartment  x-dimension  (m) 

Yield  (kg/kg);  Compartment  y-dimension  (m) 
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CO 2  Carbon  dioxide 

H20  Water 

02  Oxygen 

Denotes  velocity 


v 


Appendix  A  -  Sample  input  files 


A-l  Sample  1 

The  input  file  for  Sample  A-l  is  as  follows: 

Sample  1:  Two  room  configuration  with  t2  fire,  1MW  peak 

-  General  Simulation  and  Data  Control  - 

&EXEC  t_max=1200.0,  f ileid- 1 samplel * ,  out type- . FALSE . ,  pamb-101325 . 0  / 
&CNTR  ncomp-2,  njunc=3,  nnotbidir-O,  nfire-1,  nsurf=14, 

ncmpn=l,  nmat=l,  nfan-0,  nnode-0,  nduct-0,  ncontrol-0  / 

&FNAM  materialf ile= ’materials . txt *  / 


-  Source  Fire  and  Fuel  Load  - 

&FIRE  id-1,  f type-2,  fsubtype=3,  location-1,  t2_hrr-1000 . 0,  chi_r-0.3, 
deltah  f=28100000.0  / 


-  Compartment  Information  - 

&COMP  id-1,  compname- 1  Room  1',  dx=5.0,  dy-5.0,  dz=3.0,  vol-75.0  / 

&COMP  id-2,  compname- 1 Room  2',  dx-5.0,  dy=5.0,  dz=3.0,  vol-75.0  / 

&JUNC  id-1,  location-1,2,  height-2.13,  elevation-0.0,  area-1.948, 

bidirectional-. TRUE. ,  vf low-. FALSE.,  kloss=2. 04,  orientation-3  / 

&JUNC  id-2,  location-1, -1,  height-0.91,  elevation-1.22,  area-0.55, 

bidirect ional-. TRUE. ,  vf low-. FALSE. ,  kloss=2.16,  orientation-6  / 

&JUNC  id-3,  location-2, -1,  height-0.91,  elevation-1.22,  area-0.55, 

bidirect ional-. TRUE. ,  vflow-. FALSE. ,  kloss-2.16,  orientation-3  / 

—  Room  1  — 

&SURF  id-1,  stype-1,  location-1, -1,  orientation-1,  area-25.0,  composition-1  / 

&SURF  id-2,  stype-1,  location-1 , -1,  orientation-2,  area-25.0,  composition-1  / 

&SURF  id-3,  stype-1,  location-1,2,  orientation-3,  area— 13.05,  composition-1  / 
&SURF  id-4,  stype—1,  location— 1, -1,  orientation— 4 ,  area— 15.0,  composition-1  / 

&SURF  id-5,  stype-1,  location-1,—1,  orientation-5,  area-15.0,  composition-1  / 

&SURF  id-6,  stype-1,  location-1, -1,  orientation-6,  area-14.45,  composition-1  / 

—  Room  2  — 

&SURF  id— 7 ,  stype-1,  location— 2, -1,  orientation— 1 ,  area-25.0,  composition— 1  / 

&SURF  id— 8,  stype-1,  location-2, -1,  orientation-2,  area-25.0,  composition— 1  / 

&SURF  id-9,  stype-1,  location-2 , -1 ,  orientation-3,  area-14.45,  composition-1  / 
&SURF  id-10,  stype-1,  location-2, -1,  orientation-5,  area-15.0,  composition-1  / 

&SURF  id-11,  stype-1,  location-2, -1,  orientation-6,  area-15.0,  composition-1  / 

—  Fake  surfaces  for  Radiation  heat  flow  through  openings  — 

&SURF  id-50,  stype-1,  location-1, -1,  orientation-6,  area-0.55,  fake-. TRUE.  / 

&SURF  id-51,  stype-1,  location-1,2,  orientation-3,  area-1.948,  fake-. TRUE.  / 

&SURF  id-52,  stype-1,  location-1, -1,  orientation-3,  area-0.55,  fake-. TRUE.  / 

&CMPN  id-1,  cmpnname-1 Carbon  Steel  Boundary1,  fileread-. FALSE . , 
layers-1,  materials-1,  dx=0.0158  / 

&MTRL  id-1,  matname-' carbon  steel*  / 

- Mechanical  Ventilation - . - - 

—  None  — 
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— - - - Suppression  and  Detection 

None  — 


The  following  carbon  steel  material  property  lines  are  contained  in  the  materials .  txt 

file: 


&MTRL  matname=* carbon  steel * , e=0 , 8, rho=7860 . ,  ktable= . TRUE . , cptable= . TRUE . , 
kname- * carbon  steel  k* , cpname=* carbon  steel  cp */ 


&CURV  id=* carbon  steel  k*,xy=  200.15,81.5/ 
&CURV  id=*carbon  steel  k%xy=  298.15,71.1/ 
&CORV  id^* carbon  steel  k* ,xy=1073. 15, 29, 6/ 
&CURV  id=* carbon  steel  k* ,xy=1573 . 15, 31 . 4/ 


&CURV  id=* carbon  steel 
&CURV  id=* carbon  steel 
&CURV  id=* carbon  steel 
&CURV  id=* carbon  steel 
&CURV  id=* carbon  steel 
&CURY  id=* carbon  steel 
&CURY  id=* carbon  steel 
&CURV  id=* carbon  steel 
&CURV  id=* carbon  steel 
&CURV  id=* carbon  steel 


cp*,xy=  273.15,  438.9/ 
cp  *  , xy-  348.15,  501.6/ 
cp* ,xy=  473.15,  564.3/ 
cp  * , xy=  673.15,  627.0/ 
cp  * , xy=  873.15,  710.6/ 
cp  * , xy=  973.15,  836.0/ 
cp* ,xy=1036. 15,  836.0/ 
cp* ,xy=1041. 15, 6888 . 6/ 
ep*,xy=1046.15,  702.2/ 
cp*,xy=1223.15,  668.8/ 


A-2  Sample  2 


The  input  file  for  Sample  A-2  is  as  follows: 

Sample  2:  One  room  configuration  with  t2  fire,  800  kW  peak 

— — — _  General  Simulation  and  Data  Control  - - — - 

&EXEC  t_max=1200.0,  f ileid=* sample2 f ,  out type- . FALSE . ,  pamb=101325. 0  / 

&CNTR  ncomp=l,  njunc^2,  nnotbidir=0,  nfire=l,  nsurf=8, 

ncmpn=2,  nmat=3,  nfan-0,  nnode=0,  nduct=0,  ncontrol=l, 
nsuppress=l  / 

&FNAM  material file=*materials. txt *  / 


- - — —  Source  Fire  and  Fuel  Load  — - — - — — — - - — — — - - — — 

&FIRE  id=l,  ftype=2,  fsubtype=2,  locational,  t2_hrr=800.0,  chi_r=0.3, 
deltah  f-44000000.0  / 


— - - -  Compartment  Information  — — — — — - - - — — ; — — ’ — - 

£COMP  id-1,  compname^ * Room  1* ,  dx=5.5,  dy=3.0,  dz=3.0,  vol=49.5  / 

&JUNC  id=l,  locational, -1,  height=2.13,  elevation=0 . 0,  area=1.948, 

bidirectional=.TRUE. ,  vf low— . FALSE . ,  kloss=2.04,  orientation=4  / 

&JUNC  id=2,  location-1, -1,  height=0.3,  elevation=2.2,  area=0.36, 

bidirectional^ TRUE. ,  vflow=. FALSE, ,  kloss=2.16,  orientation=5  / 

— Room  1  — 

&SDRF  idKL,  stype^l,  locational, -1,  orientation=l,  area=16.5,  compositional  / 

&SURF  id=2,  stype=l,  location-1, -1,  orientation=2,  area-16.5,  composition^  / 

&SURF  id=3,  stype-1,  location=l, -1,  orientation=3,  area=15.G,  composition^  / 
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&SURF  id=4,  stype=l,  locational, -1,  orientation=4,  area=7.052,  compositional  / 

&SURF  id=5,  stype=l,  location=l, -1,  orientation-5,  area=20.64,  compositional  / 

&SURF  id=6,  stype=l,  locational, -1,  orientation=6,  area=12.0,  compositional  / 

— Fake  surfaces  for  Radiation  heat  flow  through  openings-- 

&SURF  id=50,  stype=l,  locational, -1,  orientation=5,  area=0.36,  fake=.TRUE.  / 

&SURF  id=51,  stype=l,  location=l, -1,  orientation=4 ,  area=1.948,  fake=.TRUE.  / 

&CMPN  id=l,  cmpnname= 1  Steel -Ins  Boundary  A*,  f ileread=. FALSE . , 
layers=2,  materials=l , 2 ,  dx=0 . 0158, 0 . 025  / 

&CMPN  id=2,  cmpnname=* Steel-Ins  Boundary  B*,  fileread=. FALSE. , 
layers=2,  materials=l, 3,  dx=0 . 025, 0 . 038  / 

&MTRL  id-1,  matname-* carbon  steel*  / 

&MTRL  id=2,  matname=* ceramic  fiberboard*  / 

&MTRL  id=3,  matname= ’ fiberfrax  blanket*  / 

-  Mechanical  Ventilation  - 

—  None  — 

-  Suppression  and  Detection  - - - - — 

&SUPR  id=l,  stype=2,  operational^ . TRUE . ,  compel,  diam=0.0005, 
flowrate=1.36,  control=l  / 

—  sprinkler  control  function  — 

&CTRL  id=l,  ctype=l,  trip=. FALSE. ,  initial^ . FALSE . ,  inputvar=2,  location=l, 
setpoint=347 . 0  / 

The  following  carbon  steel,  fiberfrax  blanket  insulation,  and  ceramic  fiberboard  material 
property  lines  are  contained  in  the  materials .  txt  file: 

&MTRL  matname=' carbon  steel  * , e=0 . 8, rho=7860 . ,  ktable=. TRUE. ,cptable=. TRUE. , 
kname=* carbon  steel  k *, cpname= *  carbon  steel  cp*/ 

&CURV  id=* carbon  steel  k*,xy=  200.15,81.5/ 

&CURV  id=* carbon  steel  k*,xy=  298.15,71.1/ 

&CURV  id=' carbon  steel  k* ,xy=1073 . 15, 29 . 6/ 

&CURV  id=* carbon  steel  k * , xy=1573 . 15, 31 . 4/ 

&CURV  id=* carbon  steel  cp*,xy=  273.15,  438.9/ 

&CURV  id=* carbon  steel  cp*,xy=  348.15,  501.6/ 

&CURV  id=* carbon  steel  cp*,xy=  473.15,  564.3/ 

&CURV  id=* carbon  steel  cp*,xy=  673.15,  627.0/ 

&CURV  id=* carbon  steel  cp*,xy=  873.15,  710.6/ 

&CURV  id-* carbon  steel  cp’,xy=  973.15,  836.0/ 

&CURV  id=* carbon  steel  cp* , xy-1036 . 15,  836.0/ 

&CURV  id=* carbon  steel  cp* ,xy=1041 . 15, 6888 . 6/ 

&CURV  id=’carbon  steel  cp’ ,xy=1046. 15,  702.2/ 

&CURV  id=* carbon  steel  cp* ,xy=1223 . 15,  668.8/ 

! Properties  from  manufacturer 

&MTRL  matname=*  fiberfrax  blanket  * , e=0 . 9, rho=128 . , ktable=.TRUE. , kname= *  fiberfrax 
blanket  k*,cp=1130./ 

&CURV  id=* fiberfrax  blanket  k*,xy=  473.15,0.04/ 

&CURV  id= * fiberfrax  blanket  k*,xy=  673.15,0.07/ 

&CURV  id=* fiberfrax  blanket  k*,xy=  873.15,0.12/ 

&CURV  id= * fiberfrax  blanket  k* , xy=1073 . 15, 0 . 18/ 

&MTRL  matname=* ceramic  fiberboard’ , e=0 . 9, rho=256. , ktable= . TRUE . , kname=* ceramic 
fiberboard  k*,cp=1130./ 

&CURV  id=* ceramic  fiberboard  k',xy=  311.15,0.043/ 
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&CURV  id- 'ceramic  fiberboard  k*,xy=  533.15,0,066/ 
&CURV  id-* ceramic  fiberboard  k*,xy=  803.15,0.104/ 
&CURV  id-* ceramic  fiberboard  k* ,xy-1366. 15, 0. 205/ 


A-3  Sample  3 


The  input  file  for  Sample  A-3  is  as  follows: 

Sample  3:  U-shaped  room  configuration  with  t2  fire 

— - - General  Simulation  and  Data  Control  — - — - — — - — 

&EXEC  t_max=1200.0,  fileid- * sample! * ,  out type-. FALSE. ,  pamb=101325.0  / 

&CNTR  ncomp-3,  njunc=5,  nnotbidir-0,  nf ire-1,  nsurf=21, 

ncmpn=2,  nmat=2,  nfan=0,  nnode-Q,  nduct=0,  ncontrol-1, 
nsuppress=0  / 

&FNAM  materialfile-'materials. txt *  / 


— - Source  Fire  and  Fuel  Load - — - ■ - 

&FIRE  id-1,  ftype-3,  location-1,  chi_r=0.3,  deltah_f-4 4 000000. 0, 
table-* SOURCE  FIRE*,  t_s tart-0.0  / 

&CURY  id-* SOURCE  FIRE*,  xy-0. 0,0.0  / 

&CURV  id- *  SOURCE  FIRE*,  xy=120, 1200 . 0  / 

&CURY  id—  *  SOURCE  FIRE*,  xy=1200, 1200 . 0  / 


--- — - -  Compartment  Information  — — - - - - - - — — - 

&COMP  id-1,  compname-* Starboard  Room*,  dx-11.0,  dy=3.0,  dz-3.0,  vol-99.0  / 
&COMP  id-2,  compname-* Center  Room*,  dx=4.0,  dy=6.0,  dz-3.0,  vol-72.0  / 
&COMP  id-3,  compname-* Port  Room*,  dx-11.0,  dy-3.0,  dz-3 . 0,  vol-99.0  / 


“  connections  between  sub- spaces  — 

SJUNC  id-1,  location-1,2,  height-3. 0,  elevation-0.0,  area-12.0, 

bidirectional-. TRUE. ,  vf low-. FALSE. ,  kloss— 1.00,  orientation— 5  / 

&JUNC  id-2,  location-2,3,  height-3, 0,  elevation-0.0,  area-12.0, 

bidirectional-. TRUE. ,  vf low— , FALSE. ,  kloss— 1.00,  orientation— 5  / 

—  openings  to  ambient  — 

&JUNC  id-3,  location-1, -1,  height-0.6,  elevation-1.9,  area-0.54, 

bidirectional-, TRUE. ,  vf low- . FALSE . ,  kloss-2.16,  orientation-6  / 

&JUNC  id-4,  location-3, -1,  height— 0.6,  elevation-1,9,  area-0.54, 

bidirectional-. TRUE. ,  vflow-. FALSE. ,  kloss=2.16,  orientation-5, 
door-1  / 

&CTRL  id-1,  ctype=l,  trip-. FALSE. ,  initial-. TRUE. ,  inputvar-2, 
location-3,  setpoint-373,15  / 

—  door  — 

& JUNC  id-5,  location— 3, -1,  height-2.13,  elevation-0.0,  area-1.948, 

bidirectional-. TRUE. ,  vflow-. FALSE. ,  kloss=2.G4,  orientation-4  / 

—  Starboard  Room  — 

&SURF  id-1,  stype-1,  location-1, -1,  orientation-1,  area-33.0,  composition-1  / 

&SURF  id-2,  stype-1 ,  location-1, -*i,  orientation-2,  area-33.0,  composition-2  / 

&SURF  id-3,  stype-1,  location-1, -1,  orientation-3,  area-9,0,  composition-2  / 

&SURF  id-4,  stype-1,  location-1, -1,  orientation-4,  area-9.0,  composition-1  / 
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&SURF  id-5,  stype=l,  locational , -1,  orientation-5,  area-21.0,  composition- 1  / 

&SURF  id-6,  stype-1,  location- 1 , -1 ,  orientation=6/  area-32.46,  composition- 2  / 

—  Center  Room  — 

&SURF  id=7,  stype-1,  locational, -1,  orientation-1,  area-18. 0,  composition=l  / 

&SURF  id-8,  stype-1,  location-1 , -1 ,  orientation-2,  area-18.0,  composition-2  / 

&SURF  id- 9,  stype-1,  location-1, -1,  orientation-3,  area=18.0,  composition-2  / 


&SURF 

id-10, 

stype-1 , 

location-1, 

-i. 

orientation-4, 

area-18 . 0, 

composition-1 

/ 

—  Port  Room 

— 

&SURF 

id=ll. 

stype-1. 

location-1, 

-i, 

orientation-1 , 

area-33 . 0, 

composition-1 

/ 

&SURF 

id-12, 

stype-1. 

location-1, 

-i. 

orientation-2 , 

area-33 . 0, 

composition-2 

/ 

&SURF 

id-13, 

stype=l. 

location-1, 

-i, 

orientation-3, 

area-9 .0, 

composition-2  / 

&SURF 

id-14, 

stype-1 , 

location-1, 

-i. 

orientation-4 , 

area-7 . 054 

,  compos it ion- 1 

/ 

&SURF 

id-15, 

stype-1 , 

location-1, 

-i. 

orientation-5, 

area-32 .46 

,  compos it ion- 2 

/ 

&SURF 

id-16, 

stype-1. 

location-1, 

-i, 

orient at ion-6, 

area-21 . 0, 

composition-1 

/ 

--  Fake  surfaces  for  Radiation  heat  flow  — 

&SURF  id-50,  stype=l,  location-1,2,  orientation-5,  area-12.0,  fake-. TRUE.  / 
&SURF  id=51,  stype-1,  location-2,3,  orientation-5,  area=12.0,  fake-. TRUE.  / 
&SURF  id-52,  stype=l,  locational, -1 ,  orientation-6,  area=0.54,  fake-. TRUE.  / 
&SURF  id-53,  stype-1,  location=3, -1,  orientation=5,  area-0.54,  fake=.TRUE.  / 
&SURF  id=54,  stype-1,  location-3, -1,  orientation-4,  area=1.948,  fake-. TRUE.  / 

&CMPN  id-1,  cmpnname- ’Gypsum  Boundary  * ,  f ileread- . FALSE . , 
layers=l,  materials=2,  dx-0.0158  / 

&CMPN  id=2,  cmpnname-* Steel  Boundary1,  fileread=. FALSE. , 
layers=l,  materials=l,  dx-0.0127  / 

&MTRL  id-1,  matname= * carbon  steel*  / 

&MTRL  id-2,  matname- *  gypsum*  / 


-  Mechanical  Ventilation  - 

—  None  — 

-  Suppression  and  Detection  - 

—  None  — 

The  following  carbon  steel  and  gypsum  material  property  lines  are  contained  in  the 
materials.txt  file: 

&MTRL  matname- ’carbon  steel ' , e=0. 8, rho-7860. ,  ktable-. TRUE. , cptable-. TRUE. , 
kname=* carbon  steel  k* , cpname-’ carbon  steel  cp*/ 

&CURV  id-* carbon  steel  k*,xy=  200.15,81.5/ 

&CURV  id-* carbon  steel  k*,xy-  298.15,71.1/ 

&CURV  id- 'carbon  steel  k’ ,xy=1073 . 15, 29. 6/ 

&CURV  id=  *  carbon  steel  k' , xy-1573 . 15, 31 . 4/ 

&CURV  id=*carbon  steel  cp*,xy=  273.15,  438.9/ 

&CURV  id=* carbon  steel  cp*,xy-  348.15,  501.6/ 

&CURV  id=* carbon  steel  cp*,xy=  473.15,  564.3/ 

&CURV  id- 'carbon  steel  cp*,xy=  673.15,  627.0/ 

&CURV  id=  *  carbon  steel  cp*,xy-  873.15,  710.6/ 

&CURV  id- ’carbon  steel  cp*,xy=  973.15,  836.0/ 

&CURV  id-' carbon  steel  cp ' , xy=1036. 15,  836.0/ 

&CURV  id-* carbon  steel  cp* , xy=1041 . 15, 6888 . 6/ 

&CURV  id- 'carbon  steel  cp’ , xy-1046 . 15,  702.2/ 

&CURV  id-’ carbon  steel  cp* , xy=1223 . 15,  668.8/ 
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&MTRL  matnairte=T gypsum* ,  e=0 , 8,  rho=211 . ,  k=0 ,75,  cp-1087  .  / 

A-4  Sample  4 

The  input  file  for  Sample  A-4  is  as  follows: 


Sample  4 ;  five  compartments  connected  with  a  corridor  and  a  t2  fire 

- - - —  General  Simulation  and  Data  Control  — - — - — — — — - 

&EXEC  t_max=1200 , 0,  fileid^' sample! f ,  outtype^. FALSE. ,  pamb=101325 . 0  / 

&CNTR  ncomp=8,  njunc=10,  nnotbidir=0,  nfire=l,  nsurf=49, 

ncmpn=l,  nmat=l,  nfan=l,  nnode=8,  nduct=7,  ncontrol=l, 
nsuppress=0  / 

&FNAM  materialfile^'materials.txt’  / 


- - —  Source  Fire  and  Fuel  Load  — — — - - - - - - 

&FIRE  id=l,  ftype=2,  fsubtype=2,  location^,  chi_r=0.3,  t2_hrr=2500.0, 
deltah  f=14 000000. 0,  t_start=0.G  / 


- — - —  Compartment  Information  — - * - — 

&CGMP  id=l ,  compname^ *  Room  1 1 ,  dx=7 . 0 ,  dy=4 , 0 ,  dz-3 . 0 ,  vol=8 4.0  / 

&COMP  id=2 ,  compname= *  Room  2  * ,  dx=5 . 0,  dy=4 . 0,  dz-3 . 0,  vol=60 . 0  / 

&CGMP  id=3,  compname= f Room  3*,  dx=8.0,  dy=4,G,  dz=3.0,  vol=96.0  / 

&CGMP  id=4,  compname= * Room  4%  dx=10.0,  dy=5.0,  dz=3,0,  vol=150 . 0  / 

&CGMP  id=5,  compname= * Room  5*,  dx=10.0,  dy=5.0,  dz=3.0,  vol=150 . 0  / 

—  Corridor  segments  — 

&COMP  id=10,  compname= * Corridor  Segment  1*,  dx=6.5,  dy-3.0,  dz=3 . 0,  vol=58.5  / 
&COMP  id=llf  compname= * Cor ridor  Segment  2’,  dx=7.0,  dy=3,0,  dz=3.0,  vol=63.Q  / 
&COMP  id=12,  compname= f Cor ridor  Segment  3*,  dx=6.5,  dy=3.0,  dz=3.0,  vol=58.5  / 

—  Doors  — 

&JUNC  id=l,  location=l, 10,  height=2.13,  elevation=0. 0,  area=1.94, 

bidirectional—, TRUE. ,  vf low- . FALSE . ,  kloss=2.04,  orientation=6  / 

& JUNC  id=2,  locations, 11,  height-2.13,  elevation=0 . 0,  area=1.94, 

bidirectional=,TRUE. ,  vflow=, FALSE. ,  kloss=2.04,  orientation=6  / 

& JUNC  id=3,  location=3, 12,  height=2.13,  elevation^O. 0,  area=1.94, 

bidirectional^ TRUE. ,  vflow=. FALSE. ,  kloss=2,04,  orientation=6  / 

&JUNC  id=4,  location=5, 12,  height=2.13,  elevation=0. 0,  area=1.94, 

bidirectional=.TRUE. ,  vflow=. FALSE. ,  kloss=2.04,  orientation=5  / 

& JUNC  id=5,  location=4, 11,  height=2.13,  elevation=0.0,  area=1.94, 

bidirectional=.TRUE. ,  vflow=. FALSE. ,  kloss=2.04,  orientation=5  / 

&  JUNC  id=6,  location-10, -1,  heights.  13,  elevation^O.O,  area=3.88, 

bidirectional^ TRUE. ,  vflow=, FALSE, ,  kloss=2.04,  orientation=5  / 

--  openings  to  ambient  — 

£ JUNC  id=7,  location=4,-l,  height=0,8,  elevation=2 .0,  area-1.6, 

bidirectional. TRUE. ,  vflow=. FALSE. ,  kloss=2.04,  orientation=6  / 

&JUNC  id=8,  location=5, -1,  heights. 8,  elevation=2, 0,  area=1.6, 

bidirectional=.TRUE. ,  vflow=. FALSE. ,  kloss=2.04,  orientation=6  / 

—  openings  between  corridor  segments  — 

&JUNC  id=9,  Iocation=10, 11,  height=3.0,  elevation=0 ,0,  area=9.G, 
bidirectional^ . TRUE , ,  vflow=. FALSE. ,  kloss=1.0,  orientation=3  / 

S JUNC  id=10,  location=ll, 12,  height=3.0,  elevation=0.0,  area=9.0. 
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bidirectional^ TRUE. ,  vflow= .  FALSE . ,  klossS.O,  orientation=3  / 


--  Room  X  — 

&SURF  idS,  stypeS,  locations, -1,  orientations,  areaS8.0,  compositions  / 

&SURF  idS,  stypeS,  locational , -1,  orientation=2,  area=28.0,  compositions  / 

&SURF  id=3,  stypeS,  locations , 2,  orientations,  areaS2.0,  compositional  / 
&SURF  id=4,  stypeS,  locations , -1,  orientations ,  areaS2.0,  composition=l  / 

&SURF  idS,  stypeS,  locational , -1,  orientations,  area=21.0,  compositions  / 

&SURF  id=6,  stypeS,  locations,  10,  orientation=6,  area=17.56,  compositions  / 
&SURF  id=7,  stypeS,  locational,  11,  orientations,  areaS.  5,  compositions  / 


—  Room  2  — 


&SURF 

id=8. 

stypeS, 

locations,  -1, 

orientations. 

areaSO .  0, 

compos  it  ionS  / 

&SURF 

id=9. 

stypeS, 

locations,  -1, 

orientations , 

areaSO .  0, 

compositions  / 

&SURF 

id=10, 

stypeS, 

locations,  3, 

orientations. 

areaS2 . 0, 

compositions  / 

&SURF 

id=ll. 

stypeS , 

locations,  -1, 

orientations, 

areaSS .  0, 

compos  it  ionS 

/ 

&SURF 

idS  2, 

stypeS, 

locations,  11, 

orientation=6, 

areaS3 . 06 

,  compos  it  ionS 

/ 

—  Room  3  — 

&SURF 

idS3 , 

stypeS , 

locations,  -1, 

orientations, 

areaS2 . 0, 

compos  it  ionS 

/ 

&SURF 

id=14 , 

stypeS, 

locations,  -1, 

orientations, 

areaS2.0, 

compos  it  ionS 

/ 

&SURF 

id=15, 

stypeS, 

locations,  -1, 

orientations , 

areaS2 . 0, 

compos  it  ionS 

/ 

&SURF 

id=16. 

stypeS , 

locations,  -1, 

orientations , 

areaS4 . 0, 

compositions 

/ 

&SURF 

id=17 , 

stypeS, 

locations,  11, 

orientation=6, 

areaS  .5, 

compos  it  ionS  / 

&SURF 

id=18. 

stypeS , 

locations,  12, 

orientation=6, 

areaS7 . 56 

,  compositions 

/ 

—  Room  4  — 

&SURF 

idS  9, 

stypeS, 

location=4, -1, 

or  ientat  ionS , 

areaSO .  0, 

compos  it  ionS 

/ 

&SURF 

idSO, 

stypeS , 

location=4, -1, 

or  ientat  ionS , 

areaSO .  0, 

compos  it  ionS 

/ 

&SURF 

idSl, 

stypeS, 

location=4, 5, 

orientations. 

areaSS .  0, 

compositions  / 

&SURF 

idS  2, 

stypeS, 

locations,  -1, 

orientations, 

areaS5. 0, 

compos  it  ionS 

/ 

&SURF 

idS  3, 

stypeS , 

locations ,  10, 

or  ientationS , 

areaS9 . 5, 

compositions 

/ 

&SURF 

idS  4, 

stypeS, 

locations,  11, 

or  ientat  ionS, 

area=8 .56, 

compositions 

/ 

&SURF 

idS  5, 

stypeS, 

locations, -1, 

orientation=6, 

areaS  9 . 2, 

compositions 

/ 

—  Room  5  — 

&SURF 

idS  6, 

stypeS , 

locations,  -1, 

orientations, 

areaSO .  0, 

compositions 

/ 

&SURF 

idS  7, 

stypeS, 

locations,  - 1 , 

orientations, 

areaSO .  0, 

compositions 

/ 

&SURF 

idS  8, 

stypeS, 

locations, -1, 

orientations , 

areaS5 . 0, 

compositions 

/ 

&SURF 

idS  9, 

stypeS , 

locations,  11, 

or  ientat  ionS, 

areaSO .  5, 

compos  it  ionS 

/ 

&SURF 

idSO, 

stypeS, 

locations,  12, 

orientat  ionS , 

areaS7 . 56 

,  compositions 

/ 

&SURF 

idSl, 

stypeS, 

locations,  -1, 

orientation=6, 

area-29.2, 

compos  it  ionS 

/ 

—  Compartment  segment  1  — 

&SURF  idS2,  stypeS,  locationSO,  -1, 
&SURF  id=33,  stype=l,  location=10, -1, 
&SURF  id=34,  stype=l,  location=10, -1 , 

—  Compartment  segment  1  — 

&SURF  id=35,  stypeS ,  locationSl, -1 , 
&SURF  id=36,  stype=l,  location=ll, -1, 

—  Compartment  segment  1  — 

&SURF  id=37,  stype=l,  location=12, -1, 
&SURF  idS8,  stypeS ,  locationS2, -1, 
&SURF  id=39,  stype=l,  location=12, -1, 


orientations , 

areaS9 .5, 

compositions 

/ 

orientation=2 , 

areaS  9 . 5, 

compositions 

/ 

orientations , 

area=3 .54, 

compos  it  ionS 

/ 

orientations , 

areaS  9 . 5, 

compositions 

/ 

orientation=2 , 

areaS9 . 5, 

compositions 

/ 

orientations , 

area-19 . 5, 

compos  it  ionS 

/ 

orientations , 

areaS9 . 5, 

compos  it  ionS 

/ 

orientation-3 , 

area=9 . 0, 

compos  it  ionS 

/ 

—  Fake  surfaces  for  Radiation  heat  flow  — 

&SURF  id=50,  stypeS,  locations,  10,  orientation=6,  areas .  94,  fake=.TRUE.  / 
&SURF  id=51,  stype=l,  locations, 11,  orientation=6,  area=1.94,  fake=.TRUE.  / 


A-7 


&SURF  id*52,  stype=l,  location*3, 12,  orientation=6,  area*1.94,  fake*. TRUE,  / 
&SURF  id*53,  stype*l,  location=4, 11,  orientation=5,  area=1.94,  fake*. TRUE.  / 
&SURF  id=54,  stype=l,  location=5, 12,  orientation=5,  area=l,94,  fake*. TRUE.  / 
&SURF  id*55,  stype=l,  location=10, -1,  orientation*4,  area=3.88,  fake*. TRUE,  / 
&SURF  id=56,  stype=l,  location*4, -1,  orientatio n*6,  area*0.8,  fake*. TRUE.  / 
&SURF  id=57,  stype*l,  location*5, -1,  orientation=6,  area=0.8,  fake*. TRUE.  / 
&SURF  id=58,  stype=l,  location=10, 11,  orientation=3,  area=9.0,  fake*. TRUE.  / 

5 SURF  id*59,  stype=l,  location*!!, 12,  orientation=3,  area*9,0,  fake*. TRUE.  / 

&CMPN  id=l,  cmpnname* * Gypsum  Boundary*  ,  f ile read*. FALSE. , 
layers*!,  materials*2,  dx=0.0158  / 

&MTRL  id=2,  mat name* 1  gypsum1  / 


- -  Mechanical  Ventilation  - - - - — 

--  nodes  — 

&RNOD  id-1,  nducts=2,  ducts=l,l,  terminal*.TRUE. ,  elevation=2.5, 
loss*0. 0,0.0  / 

&RNOD  id=2,  nducts*2,  ducts=l,4,  terminal*. FALSE. ,  elevation=2 . 5, 
loss=l. 1,1.1  / 

&RN0D  id*3,  nducts-2,  ducts=2,2,  terminal*.TRUE. ,  elevation=2 . 5, 
loss*0, 0,0.0  / 

£RN0D  id=4,  nduets*3,  ducts*2,4,5,  terminals  FALSE. ,  elevation=2 . 5, 
loss— 1 . 63, 1 . 63, 1 . 63, 0 . 75, 1 . 63, 0 . 75  / 

&RNOD  id=5,  nducts*2,  ducts=3,3,  terminal*.TRUE. ,  elevation=2 . 5, 
loss=0. 0,0.0  / 

&RNOD  id*6,  nducts=2,  ducts*3,6,  terminals  FALSE. ,  elevation=2 . 5, 
loss=l .1,1.1  / 

&RNOD  id*7,  nducts*3,  ducts=5,6,7,  terminal— . FALSE . ,  elevation=2 . 5, 
loss=0. 75, 1.63, 0.75, 1.63, 1.63, 1.63  / 

&RNOD  id=8,  nducts=2,  ducts=~l,7,  terminal=.TRUE. ,  elevation-2.5, 
loss=0. 0,0.0  / 

—  ducts  — 

&RDCT  id*l,  node*l,2,  area=0.3721,  perimeter=2 . 44,  length=2.0, 
roughness=0. 0005  / 

&RDCT  id*2,  node*3,4,  area*0.3721,  perimeter=2. 44,  length=2.0, 
roughness=0. 0005  / 

&RDCT  id*3,  node=5, 6,  area=0.3721,  perimeter-2.44,  length-2.0, 
roughness=0.0005  / 

&RDCT  id*4,  node=2,4,  area*0.3721,  perimeter=2 . 44,  length*6.0, 
roughness=0.0005  / 

&RDCT  id*5,  node*4,7,  area=0.3721,  perimeter*2. 44,  length=3.0, 
roughness=0. 0005  / 

&RDCT  id=6,  node*6,7,  area*0.3721,  perimeter*2 . 44,  length=2.0, 
roughness*0 .0005  / 

£RDCT  id*7,  node=7,8,  area=0,3721,  perimeter=2 . 44,  length=6.0, 
roughness*0,G005,  fan*l  / 

—  fan  — 

&RFAN  id*!,  f ileread— . FALSE . ,  ftype*2,  full*1.83,  shutoff*180, 
reverse* . FALSE . ,  control*!  / 

&CTRL  id*l,  ctype*l,  trip* . FALSE . ,  initial*. TRUE . ,  inputvar*!, 
setpoint=10000.00  / 


- — - -  Suppression  and  Detection 

--  None  — 
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The  following  gypsum  material  property  line  is  contained  in  the  materials .  txt  file: 


&MTRL  matname= ' gypsum ' , e=0 . 8, rho=211 . , k=0 . 75, cp=1087 . / 
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